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1 6.  Abstract  (1MC  200  wort*) 


A  procedure  for  three-dimensional  nonlinear  material  modelling  of  fiber-reinforced  laminated  composites 
is  presented.  The  material  modelling  procedure  has  a  two-level  hierarchical  structure.  At  the  bottom  level, 
constitutive  information  about  the  fiber  and  the  matrix  phases  are  synthesized  using  a  micromechanical 
model  to  yield  the  effective  stress-strain  response  of  a  unidirectional  lamina.  At  the  top  level,  a  three-di¬ 
mensional  lamination  scheme  is  employed  which  assembles  the  laminae  within  a  sublaminate,  and  delivers 
the  effective  stress-strain  response  of  the  sublaminate.  Local  stresses  and  strains  in  a  lamina  or  in  fiber 
and  matrix  phases  can  be  recovered  from  the  effective  values  at  any  stage.  The  material  modelling  proce¬ 
dure  enables  the  use  of  standard  displacement-based  finite  elements.  The  matrix  material  is  characterized 
using  nonlinear-elastic  Ramberg-Osgood  relations.  Micromechanical  failure  criteria  are  used  for  determin¬ 
ing  various  modes  of  failure,  including  compression  kink-banding.  The  accuracy  of  the  micromechanical 
model  is  demonstrated  by  comparing  its  predictions  with  results  from  other  micromechanical  models  and 
experimental  data.  Examples  are  also  presented  for  laminated  structures;  the  results  are  in  good  agreement 
with  analytical  and  experimental  results  available  in  the  literature. 
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CHAPTER  1 


INTRODUCTION 


1.1  General 

A composite  maybe  defined  as  (Lee  [1])  a  multi-phase  material  formed  by  com¬ 
bining  materials  differing  in  composition  and  form,  in  a  manner  that  they  remain  bonded 
together  and  retain  their  identities  and  properties.  The  constituents  of  a  composite  act  in 
concert  to  provide  improved  specific  or  synergistic  properties  not  obtainable  by  one  of  the 
components  acting  alone. 


Composite  materials  are  ideally  suited  for  products  and  applications  where  high  stiff- 
ness-to- weight  and  strength-to-weight  ratios  are  desirable,  such  as  in  aircraft  and  spacecraft 
structures.  Composites  are  resistant  to  corrosion  and  fatigue,  and  typically  have  good  impact 
tolerance.  All  these  attributes  make  die  use  of  composites  an  attractive  option  for  a  wide 
range  of  applications,  and  consequently  there  has  been  a  steady  rise  in  die  use  of  composites 
in  the  aerospace  industry,  automobile  industry,  and  sporting  goods  industry.  A  detailed  ac¬ 
count  of  various  applications  of  composite  materials  may  be  found  in  [2]. 

Composites  may  be  broadly  classified  as  (Lee  [1]): 

•  Fibrous  (composed  of  fibers  in  a  matrix) 

•  Laminated  (layered,  of  different  materials) 

•  Particulate  (composed  of  particles  or  flakes  in  a  matrix) 

•  Hybrid  (combinations  of  any  of  the  above  types) 

The  present  study  focuses  on  developing  analysis  procedures  for  unidirectionally 
fiber-reinforced  laminated  composites.  Such  a  composite  consists  of  unidirectionally 
fiber-reinforced  laminae,  stacked  together,  with  the  fiber  direction  in  each  lamina  being 
possibly  different,  to  achieve  the  desired  properties  in  each  direction. 
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A  typical  unidirectional  fiber-reinforced  lamina  is  about  0.005-0.01  inches  in  thickness 
and  its  cross-section  may  contain  numerous  fine  fibers  arranged  in  a  random  order;  al¬ 
though  to  facilitate  analytical  modelling,  it  is  conveniently  assumed  that  the  fibers  are  ar¬ 
ranged  in  neat  packing-geometries,  such  as  triangular,  square,  or  hexagonal. 

When  structural  analysis  procedures  for  structures  made  of  traditional  materials  of 
construction  such  as  steel  are  compared  with  analysis  procedures  for  laminated  structures, 
two  important  differences  stand  out: 

•  Unlike  most  metals,  a  fiber-reinforced  lamina  tends  to  be  trans¬ 
versely  isotropic,  with  the  axis  of  isotropy  being  parallel  to  the  fiber 
direction.  Furthermore,  the  lamina  can  exhibit  nonlinear  response, 
particularly  when  subjected  to  shear  stress  in  the  plane  of  lamina. 

•  In  shell  and  plate  structures  made  of  metals,  there  is  rarely  a  varia¬ 
tion  of  material  properties  in  the  thickness  direction.  Whereas,  in 
laminated  plates  and  shells,  the  material  properties  vary  from  ply  to 
ply  because  of  the  different  orientation  angle  of  each  ply. 

These  two  important  differences  affect  the  way  in  which  the  analysis  of  laminated 
composites  is  accomplished:  the  first  bears  on  the  constitutive  modelling  of  a  lamina,  and 
the  second  on  the  structural  modelling  of  the  laminated  structure.  These  two  important 
considerations  are  described  in  greater  detail  in  the  following  sections  in  order  to  bring 
out  the  motivation  and  aims  of  the  present  study. 

1.1.1  Constitutive  Modelling  of  the  Lamina 

This  entails  a  description  of  the  material  behaviour,  including  a  suitable  strength 
theory  or  criteria  employed  to  define  the  onset  of  various  modes  of  failure  in  the  material. 
The  constitutive  theories  used  to  model  the  behaviour  of  composite  materials  may  be  clas¬ 
sified  as  Micromechanical  ox  Macromechanical  in  nature.  Micromechanical  theories  pre¬ 
dict  the  overall  behaviour  of  a  composite  lamina  from  the  known  properties  of  the  indi- 
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vidual  constituents,  i.e.  the  fiber  and  the  matrix,  and  their  detailed  interaction.  This 
usually  involves  considering  a  Representative  Volume  Element  of  the  composite  consisting 
of  the  fiber  phase  and  the  matrix  phase.  Principles  of  mechanics  are  applied  to  determine 
the  local  stress  and  strain  fields  within  the  representative  element  and  the  local  fields  are 
averaged  to  yield  the  effective  or  equivalent  homogeneous  properties  of  the  composite 
lamina.  Macromechanical  theories  on  the  other  hand  consider  the  lamina  to  be  a  homo¬ 
geneous  anisotropic  material  whose  properties  are  to  be  determined. 

One  advantage  of  using  a  micromechanical  theory  is  that  no  apriori  knowledge  of 
the  lamina  response  is  needed,  since  this  response  is  predicted  from  the  properties  of  the 
constituent  fiber  and  matrix  phases.  The  constituents  themselves  are  homogeneous,  and 
the  vast  amount  of  knowledge  accumulated  regarding  the  behaviour  of  homogeneous 
materials  can  be  directly  employed  to  construct  the  effective  response  of  a  composite 
lamina.  A  second  advantage  is  that  since  only  the  properties  of  the  constituents  are 
required,  the  effect  of  using  different  fiber  matrix  combinations  and  different  fiber 
volume  fractions  can  be  easily  studied.  Also,  the  interaction  of  the  fiber  and  the  matrix 
material  can  be  accurately  modelled. 

1.1.2  Structural  Modelling  of  Laminates 

The  laminated  nature  of  composite  structures  resulting  from  the  stacking  together 
of  laminae  poses  a  problem  because  of  the  discontinuities  in  material  properties  across 
the  interface  of  two  laminae  with  different  fiber  orientations.  To  address  this  problem, 
various  laminated  plate  theories  have  been  proposed.  The  simplest  of  these  theories  is  the 
Classical  Laminated  Plate  Theory  (CLPT),  which  is  an  extension  of  the  classical  Kirchhoff 
plate  theory  to  laminated  composites.  Similar  extensions  to  Mindlin  plate  theory  have 
been  proposed,  and  are  referred  to  as  First-order  Shear  Deformation  Theories  (FSDT) 
for  laminates.  CLPT  and  FSDT  predict  the  response  of  thin  to  moderately  thick  laminates 
quite  well.  The  perceived  need  for  further  refinement  and  accuracy  in  the  analysis  of  thick 
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laminates  led  to  the  development  of  Higher-order  Shear  Deformation  Plate  Theories 
(HSDT),  in  which  the  assumed  displacement-field  is  of  a  still  higher  polynomial  order  in 
the  thickness  direction  than  in  the  FSDT.  Aside  from  the  finite  elements  based  on  these 
various  laminated  plate  theories,  various  types  of  hybrid  elements  based  on  hybrid 
stress-strain  formulations  have  also  been  proposed.  It  is  easy  to  see  that  with  this 
approach,  in  order  to  extend  the  capabilities  of  the  existing  finite  element  analysis 
packages  to  include  laminated  composite  structures,  a  major  effort  would  be  required  to 
develop  a  separate  element  library,  containing  elements  specifically  for  the  analysis  of 
laminated  composite  structures. 

The  present  investigation  is  aimed  at  developing  efficient  and  accurate  procedures 
for  constitutive  modelling  of  fiber-reinforced  composites  which  enable  the  use  of  stan¬ 
dard  displacement-based  finite  elements  to  model  laminated  structures.  The  aims  and 
objectives  of  he  present  study  are  formally  outlined  in  the  next  section. 

1.2  Objectives 

The  aims  of  this  study  are  as  follows: 

•  To  model  the  response  of  laminated  composites  using  the  properties 
and  the  state  of  its  basic  constituents,  i.e.  fiber  and  matrix.  This 
basically  means  that  some  sort  of  micromechanical  constitutive 
theory  is  needed  to  describe  the  material  behaviour  of  the  compos¬ 
ite.  As  previously  stated,  micromechanical  theories  have  certain 
advantages  over  macromechanical  theories,  since  they  provide  de¬ 
tailed  information  about  the  state  of  the  basic  constituents,  and  offer 
greater  flexibility  in  modelling  various  types  of  nonlinearities  ob¬ 
served  in  composites. 

•  It  is  desired  that  the  micromechanical  theory  chosen  should  be  accu¬ 
rate,  and  at  the  same  time,  simple  enough  to  render  its  incorporation 
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in  a  general  structural  analysis  finite  element  package  to  be  practi¬ 
cal. 

•  The  micromechanical  model  chosen  should  have  sufficient  flexibil¬ 
ity  to  enable  the  modelling  of  different  types  of  nonlinearity  and 
damage  mechanisms. 

•  The  modelling,  or  the  analysis  procedure,  should  be  such  that  the 
equilibrium  requirements  on  tractions  at  the  laminae  interfaces  are 
not  violated. 

•  This  material  modelling  of  the  laminated  composite  is  to  be  done 
in  such  a  manner  that  the  standard  displacement-based  finite  ele¬ 
ments  can  be  used  for  analyzing  laminated  composite  structures. 

Having  outlined  the  aims  of  the  present  study,  a  review  of  the  constitutive  theories 
employed  to  describe  the  behaviour  of  composites  will  be  presented  next. 

1.3  Constitutive  Theories  for  Composites 

In  the  present  study  it  was  decided  to  employ  a  micromechanical  theory;  therefore 
in  the  following  review  of  constitutive  theories,  greater  emphasis  is  placed  on  microme¬ 
chanical  theories. 

1.3.1  Micromechanical  Theories 

The  field  of  micromechanics  has  been  an  area  of  much  research  activity.  The  bulk  of 
the  earlier  literature  is  devoted  to  micromechanical  theories  aimed  at  determination  of 
effective  elastic  moduli  of  a  composite  lamina  from  the  properties  of  constituent  phases.  The 
simplest  micromechanical  model  for  determination  of  effective  elastic  moduli  is  the 
Mechanics  of  Materials  model  (Fig.  1.1),  in  which  the  constituent  phases  are  assumed  to  be 
either  connected  in  a  series  (Reuss)  arrangement,  or  in  a  parallel  (Voigt)  arrangement  The 
latter  arrangement  of  constituents  is  also  commonly  referred  to  as  the  Rule  of  Mixtures 
approach.  A  description  of  the  Mechanics  of  Materials  approach  can  be  found  in  the  works 
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of  Jones  [3]  and  Harris  [4].  The  axial  stiffness  and  the  axial  Poisson’s  ratio  of  a  unidirectional 
composite  is  quite  accurately  predicted  by  the  Mechanics  of  Materials  approach.  To  deter¬ 
mine  the  transverse  direction  properties,  the  fibers  and  the  matrix  are  assumed  to  be 
connected  in  series;  the  resulting  estimates  of  transverse  modulus  are  considerably  lower 
than  the  experimentally  observed  values,  and  in  fact,  die  estimates  so  obtained  define  the 
Reuss  Lower  Bound  on  the  transverse  modulus  of  a  unidirectional  composite.  To  remedy  this 
situation  Combination  Type  models  have  been  proposed,  in  which  the  constituents  are 
arranged  in  such  a  manner  that  a  portion  of  them  is  acting  in  series  and  the  remainder  is 
acting  in  parallel  with  one  another.  An  example  of  this  approach  is  the  micromechanical 
model  proposed  by  Shaffer  [12].  Shaffer’s  model  resulted  in  some  improvement  over  the 
Reuss  Lower  Bound  for  transverse  modulus.  The  model  could  also  be  extended  to  consider 
the  response  of  the  composite  in  the  plastic  range,  i.e  after  the  onset  of  yielding  in  the  matrix. 

In  other  related  works  by  Evkall  [15]  and  Abolin’sh  [16],  the  increased  transverse 
direction  stiffness  of  unidirectional  composites  as  compared  to  the  prediction  by  the 
Mechanics  of  Materials  approach  was  investigated.  Abolin’sh  [16]  pointed  out  that 
when  a  unidirectional  composite  is  stressed  in  the  direction  transverse  to  the  fibers,  the 
matrix  is  restrained  in  the  fiber  direction,  because  it  is  attached  to  the  fiber  which  has 
high  axial  stiffness.  Thus  the  apparent  stiffness  of  the  matrix  is  greater  than  the  stiffness 
of  the  pure  matrix  in  bulk.  When  this  apparent  or  effective  stiffness  of  the  matrix  phase 
is  taken  into  account,  the  resulting  transverse  direction  stiffness  of  composite  is  greater 
than  that  given  by  the  Reuss  Lower  Bound.  This  approach  is  referred  to  as  the  Restrained 
Matrix  Concept. 

Other  available  schemes  include  the  Self-consistent  Field  method, which  includes 
the  Two-phase  model  [19,20,21]  (Fig.  1.2)  and  the  Three-phase  model  [17,18,22,23]  (Fig. 
1.3).  In  the  two-phase  model  version  of  the  Self-consistent  scheme,  the  effective 
properties  of  the  composite  are  determined  by  solving  two  problems.  First,  an  inclusion 


6 


of  the  fiber  phase  is  considered  to  be  imbedded  in  a  medium  whose  properties  are 
equivalent  to  the  unknown  properties  of  the  composite.  A  uniform  state  of  stress  or 
strain  is  then  prescribed  at  a  distance  sufficiently  far  from  the  inclusion.  The  average 
stresses  and  strains  in  the  inclusion  are  expressible  in  terms  of  the  elastic  moduli  of  the 
composite.  The  same  problem  is  then  repeated  for  the  inclusion  made  up  from  the 
matrix  phase  of  the  composite.  The  average  stresses  or  strains  in  the  two  phases  are  then 
combined  to  yield  the  effective  moduli  of  the  composite.  Among  the  short-comings  of 
the  Two-phase  model  cited  in  literature  are:  that  it  is  appropriate  only  for  low  fiber 
volume  fractions,  it  results  in  poor  representation  of  the  microstructure  of  the 
composite,  and  it  does  not  take  into  account  the  interaction  among  the  constituent 
phases. 

The  Three-phase  variant  of  the  Self-consistent  field  method,  as  applicable  to 
fiber-reinforced  composites,  is  attributed  to  Hill  [17,18]  in  which  solutions  for  four  of  the 
five  elastic  moduli  associated  with  transversely  isotropic  composite  lamina  were  given. 
The  correct  solution  for  the  fifth  property  (  shear  modulus  in  the  plane  of  isotropy)  was 
given  later  by  Christensen  and  Lo  [23].  The  model  consists  of  a  fiber  surrounded  by  a 
matrix  shell,  embedded  in  a  medium  whose  properties  are  equivalent  to  the  composite. 
This  model  is  a  more  reasonable  representation  of  the  microstructure  of  the  composite 
lamina,  as  compared  to  the  IWo-phase  model. 

Hashin  and  Rosen  [25]  presented  the  solution  for  effective  moduli  of  a  composite 
comprised  of  circular  fibers  of  varying  diameters  arranged  in  a  random  manner.  The 
resulting  composite  is  transversely  isotropic  in  an  overall  sense.  The  assumed  micro¬ 
structure  and  the  corresponding  representative  volume  element  are  shown  in  Fig.  1.4. 
The  model  is  referred  to  in  the  literature  as  the  Composite  Cylinder  Assemblage  (CCA) 
model.  Closed-form  expressions  for  E\  ,£23  ,Gi  and  vi  were  presented;  closed  form 
expressions  for<?23  *£2  and  V23  could  not  be  obtained,  and  these  quantities  could  only 
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be  bounded.  The  results  presented  by  Hashin  and  Rosen  [25]  contained  some  errors; 
these  were  subsequently  presented  in  correct  form  by  Dow  and  Rosen  [26],  and  by 
Hashin  [10,28].  An  interesting  relation  between  the  Three-phase  model  and  the  CCA 
model  was  observed  by  Hashin  [9]:  the  expressions  for  .^23  >G\  and  vi  obtained 
by  Hermans  [22]  for  the  Three-phase  model  are  exactly  the  same  as  given  by  the  CCA 
model.  Hermans’  derivation  of  the  expression  for  G23  contained  some  errors  which 
were  removed  later  by  Christensen  and  Lo  [23],  who  gave  the  result  for  G&  in  the  form 
of  a  root  of  a  quadratic.  It  was  concluded  by  Hashin  [27]  and  Christensen  and  Lo  [23] 
that  the  Three-phase  model  and  the  CCA  model  are  related  to  one  another,  although  a 
rigorous  proof  of  the  equivalence  of  the  two  methods  was  not  given. 

Various  bounds  on  the  effective  elastic  moduli  of  fiber-reinforced  composites 
have  been  proposed.  These  include  the  bounds  given  by  the  CCA  model  [25,26]  and  the 
general  bounds  established  by  Hill  [17]  and  Hashin  [24]. 

To  study  the  effect  of  fine  details  of  the  microstructure  of  a  composite,  such  as  the 
shape  of  the  fibers  and  the  fiber  packing  geometry,  upon  the  effective  elastic  moduli, 
various  numerical  studies  have  been  conducted,  e.g.  Foye  [29],  Adams  and  Doner  [3031] 
and  Chen  and  Cheng  [3233].  These  numerical  procedures  usually  employ  either  the 
finite  element  method  or  the  finite  difference  procedure  in  order  to  model  a  representa¬ 
tive  volume  element  of  the  composite  in  detail.  The  effective  moduli  of  the  composite  are 
determined  by  applying  appropriate  displacement  constraints  and  loadings  on  the  bound¬ 
aries  of  the  representative  volume  element  to  correspond  to  uniform  nominal  strain.  Such 
numerical  procedures  are  for  too  complicated  and  computationally  intensive  for  the 
purpose  of  structural  analysis  of  composite  structures.  However,  they  provide  useful  and 
accurate  results  which  can  be  used  to  assess  the  accuracy  of  other  proposed  micromechan¬ 
ical  models.  Also,  they  provide  useful  insights  into  the  behaviour  of  composites  as 
affected  by  the  interaction  of  the  fiber  and  the  matrix  at  the  fiber/matrix  interface. 
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This  completes  the  over-view  of  micromechanical  models  used  for  predicting  the 
linear-elastic  properties  of  fiber-reinforced  composites.  For  detailed  discussion  and 
description  of  various  models,  the  reader  is  referred  to  the  works  of  Chamis  and 
Sendekyj  [8]  and  comprehensive  reviews  by  Hashin  [9]  and  McCullough  [11].  Next,  an 
account  of  the  micromechanical  models  which  have  been  used  to  describe  the  nonlinear 
behaviour  of  fiber-reinforced  composites  is  given. 

The  earliest  work  to  describe  the  nonlinear  behaviour  of  fiber-reinforced  materi¬ 
als  in  the  context  of  micromechanics  is  by  Hill  [18].  Hill  employed  the  Self-consistent 
field  method  to  obtain  the  instantaneous  moduli  of  a  composite  cylinder  under  axisym- 
metric  loading  conditions;  bounds  on  the  instantaneous  moduli  were  also  obtained. 

A  number  of  finite  element  based  micromechanical  studies  have  been  conducted 
by  various  researchers  to  investigate  the  nonlinear  behaviour  of  fiber-reinforced  compos¬ 
ites,  e.g.  Foye  [35],  Adams  [36],  and  Adams  and  Crane[37].  Initial  yield  surfaces  for 
various  metal  matrix  composites  subjected  to  mechanical  and  thermal  loadings  were 
investigated  by  Dvorak  et  al.[39,40],  using  the  finite  element  method.  Their  results 
showed  that  matrix  yielding  in  fibrous  composites  usually  starts  at  the  fiber/matrix 
interface.  The  initial  yield  surface  was  found  to  be  an  irregular  ellipsoid,  with  its  longest 
axis  inclined  towards  the  hydrostatic  axis  in  the  stress  space.  The  important  factors 
affecting  the  onset  of  yielding  were  found  to  be  the  ratio  of  fiber  modulus  to  matrix 
modulus  and  the  fiber  volume  fraction.  It  was  found  that  relatively  small  temperature 
changes  can  cause  yielding.  An  interesting  find  reported  by  them  was  that  unlike  homoge¬ 
neous  metals,  the  application  of  hydrostatic  loading  results  in  yielding  of  the  composite. 

Dvorak  and  Bahei-El-Din  [41]  employed  the  Self-consistent  field  method  to 
study  the  elastic-plastic  behaviour  of  fibrous  composites.  The  Self-consistent  scheme 
was  modified  to  correct  some  of  the  problems  associated  with  the  TWo-phase  version  of 
this  method,  as  reported  by  Hutchinson  [38],  who  observed  high  estimates  of  initial  yield 
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stresses  and  low  estimates  of  plastic  strains.  The  unmodified  and  modified  Self-consistent 
schemes  gave  similar  results  for  axisymmetric  loading.  For  the  case  of  longitudinal  shear 
loading  the  modified  self-consistent  scheme  performed  well,  whereas  the  unmodified 
scheme  gave  values  of  initial  yield  stresses  substantially  greater  than  the  finite  element 
predictions.  The  extension  of  the  Self-consistent  schemes  to  handle  non-symmetric 
loadings  was  found  to  be  difficult. 

In  an  effort  to  devise  a  micromechanical  model  which  is  simple,  yet  applicable  to 
general  stress  states,  Dvorak  and  Bahei-El-Din  [41,42]  proposed  the  Vanishing  Fiber 
Diameter  Model ,  in  which  the  composite  is  modelled  as  a  continuum  reinforced  by  fibers 
of  vanishingly  small  diameter  without  changing  the  fiber  volume  fraction.  In  this  manner 
the  axial  strain  in  both  phases  was  constrained  to  be  the  same,  thus  incorporating  the 
restrained  matrix  approach.  The  transverse  direction  stresses  in  the  fiber  and  the  matrix 
were  assumed  to  be  the  same  —  note  that  this  is  tantamount  to  arranging  the  fiber  and 
matrix  phases  in  series  for  loading  in  transverse  direction,  hence  their  observation  [43] 
that  some  of  the  moduli  given  by  this  model  at  times  fell  below  the  Hashin-Rosen 
bounds[25]. 

Min  and  Crossman  [50]  presented  a  Mechanics  of  Materials  based  micromechani¬ 
cal  model  for  describing  the  nonlinear  behaviour  of  fiber-reinforced  composites  sub¬ 
jected  to  thermomechanical  loadings.  The  fiber  and  the  matrix  phases  were  assumed  to  be 
in  a  state  of  plane-stress.  The  predictions  from  the  model  showed  reasonable  agreement 
with  the  experimental  data  in  most  cases. 

Sun  and  Chen  [55]  proposed  a  simple  micromechanical  model  to  investigate 
elastic-plastic  behaviour  of  fibrous  composites.  The  fiber  was  assumed  to  be  linearly 
elastic,  and  the  matrix  was  modelled  as  elastic-plastic  material  following  the  J2  flow- 
rule.  A  macromechanical  orthotropic  plasticity  model  was  also  proposed  to  describe  the 
lamina  response.  The  authors  suggested  that  the  micromechanical  model  could  be  used 
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to  predict  the  responses  of  desired  fiber  matrix  combinations,  following  which  the  macro¬ 
mechanical  model  could  be  calibrated  and  tuned  to  match  the  micromechanical  predic¬ 
tions.  The  macromechanical  model  was  used  in  a  finite  element  based  on  First-order 
Shear  Deformation  Theory  to  predict  the  response  of  notched  and  tapered  laminates  [56]. 
In  the  development  of  the  micro-model,  the  constituent  phases  were  assumed  to  be  in  the 
state  of  plane-stress;  thus  the  three-dimensional  nature  of  the  stress  field  was  not 
accounted  for.  The  matrix  restraining  effect  of  the  fibers  was  also  not  fully  accounted  for. 

A  micromechanical  theory  based  on  the  analysis  of  a  typical  repeating  cell 
representing  a  unidirectional  composite  was  presented  by  Aboudi  [7,51,52,53,54].  The 
composite  is  assumed  to  be  comprised  of  square  fibers  arranged  in  a  doubly  periodic 
array  as  shown  in  Fig.  (1.5).  The  representative  cell  is  divided  into  four  subcells;  one 
fiber  subcell  and  three  matrix  subcells.  The  theory  proceeds  by  assuming  a  linear 
displacement  field  in  each  subcell,  and  the  continuity  of  displacements  and  tractions  at 
the  fiber/matrix  interface  is  satisfied  in  an  average  sense.  The  theory  has  been 
successfully  used  to  describe  the  nonlinear  response  of  fiber-reinforced  composites  by 
using  the  unified  viscoplastic  theory  of  Bodner  and  Partom  [69]  to  describe  the  matrix 
material.  The  theory  is  referred  to  as  Method  of  Cells,  and  is  quite  accurate,  but  requires 
a  considerable  amount  of  calculation. 

A  combination  type  micromechanical  model  based  on  Mechanics  of  Materials  was 
presented  by  Pecknold  [57].  The  model  is  quite  simple  in  description,  and  the  predictions 
of  elastic  moduli  showed  excellent  agreement  with  results  from  detailed  numerical  analy¬ 
ses.  The  model  can  portray  the  three-dimensional  stress  state  of  the  constituent  phases, 
and  can  be  used  for  general  loading  conditions.  The  three-dimensional  matrix  restraint 
effect  is  also  incorporated  in  this  model. 
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1.3.2  Macromechanical  Theories 


In  macromechanical  analyses  of  laminates,  a  lamina  is  modelled  as  an  anisotropic 
homogeneous  material;  consequently  it  can  not  provide  information  about  the  state  of 
stress  and  strain  in  the  constituent  phases.  Hill’s  anisotropic  plasticity  theory  [58,59], 
originally  intended  to  model  the  orthotropic  behavior  of  rolled  metals,  has  been  modified 
and  extended  by  many  researchers  to  describe  the  behaviour  of  fibrous  composites.  Hill 
used  a  general  quadratic  yield  criterion;  however,  the  yield  criterion  was  not  based  on 
stress  invariants,  and  consequently  the  criterion  is  applicable  only  in  the  principal  axes  of 
the  material.  The  usual  idealizations  that  hydrostatic  stress  states  do  not  influence  the 
yield  criterion  and  that  the  plastic  flow  is  incompressible  were  retained. 

Griffin  et  al.  [60]  developed  a  three-dimensional  finite  element  analysis  program 
for  analysis  of  composite  materials.  Hill’s  anisotropic  yield  criterion  was  employed  along 
with  isotropic  hardening.  To  account  for  nonlinear  deformation  characteristics  of  com¬ 
posites,  Ramberg-Osgood  type  relations  were  assumed  with  no  stress  interaction. 

Petit  and  Whddoups  [61]  proposed  a  method  for  predicting  the  nonlinear  beha¬ 
viour  and  strength  of  laminated  composites.  The  response  of  the  lamina  to  each  stress 
component  was  assumed  to  be  independent  of  stresses  in  other  directions  (such  materials 
are  referred  to  as  quasi-linear  materials).  The  bimodularity  of  the  lamina  in  tension  and 
compression  was  also  modelled.  Maximum  Strain  criteria  were  used  to  determine  failure 
in  a  lamina.  After  failure,  the  lamina  was  gradually  unloaded  in  the  mode  in  which  the 
failure  had  occurred.  The  analysis  was  conducted  by  combining  the  laminae  using  the 
framework  of  Classical  Laminated  Plate  Theory.  Their  analysis  captured  the  nonlinear 
behaviour  of  laminates  fairly  well,  but  the  strength  predictions  made  by  their  procedure 
ranged  from  good  to  poor.  The  authors  emphasized  the  need  for  conducting  more  tests  to 
give  their  procedure  a  statistical  basis  as  well.  The  presentation  of  results  and  other  data 
pertaining  to  the  behaviour  of  laminates  investigated  by  them  was  quite  comprehensive. 


12 


and  their  work  still  remains  as  one  of  the  well  documented  investigations  into  the 
behaviour  of  laminated  composites. 

Hahn  and  Bai  [62,63]  proposed  constitutive  equations  for  a  lamina  in  plane-stress 
from  consideration  of  the  complementary  energy  density  function.  The  complementary 
energy  density  was  expressed  as  a  fourth  order  polynomial  in  the  in-plane  stress  compo¬ 
nents.  Modifications  were  made  to  the  general  expression  so  that  only  shear  nonlinearity 
remained  and  the  response  to  the  other  stress  components  was  linear-elastic.  The  model 
has  been  successfully  used  to  describe  the  nonlinear  behaviour  of  laminates. 

Hashin  et  al.  [64]  modelled  the  nonlinear  response  of  unidirectional  composites 
subjected  to  transverse  and  shear  loadings,  using  Ramberg-Osgood  type  relations  for  the 
transverse  direction  and  for  the  case  of  shear  loading.  The  constitutive  relations  of  the 
lamina  were  incorporated  into  the  Classical  Laminated  Plate  Theory,  to  obtain  responses 
of  laminates  with  various  stacking  sequences.  Their  predictions  of  nonlinear  response  and 
failure  were  in  good  agreement  with  predictions  from  other  methods  and  experimental 
data. 


Dvorak  and  Rao  [65]  developed  an  axisymmetric  plasticity  theory  for  fiber-rein¬ 
forced  composites.  The  yield  function  was  constructed  from  stress  invariants  valid  for 
transversely  isotropic  materials.  The  composite  was  assumed  to  be  plastically  deformable  in 
the  fiber  direction,  and  capable  of  exhibiting  plastic  dilatation.  A  simple  hardening  rule  and 
an  an  associated  flow  rule  were  also  proposed 

Pipkin  and  Rogers  [66]  presented  a  continuum  theory  for  finite  plane  deformations 
of  composites  reinforced  by  inextensible  fibers.  The  composite  was  assumed  to  be 
incompressible  and  the  fibers  were  assumed  to  be  continuously  distributed.  These 
assumptions  are  restrictive  enough  to  render  some  of  the  problems  of  deformation 
kinematically  determinate,  and  if  a  displacement  boundary  value  problem  is  prescribed, 
the  displacement  field  within  the  composite  is  determinable  without  the  consideration  of 
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constitutive  relations  for  the  composite.  However,  the  solution  of  a  traction  boundary 
value  problem  would  require  information  about  the  constitutive  relations  for  the  compos¬ 
ite.  The  theory  was  later  extended  by  Mulhem  et  al.  [67]  to  include  plastic  deformations. 

Mulhera  et  al.  [67]  idealized  the  composite  as  transversely  isotropic,  comprised  of 
rigid-perfectly-plastic  matrix  reinforced  by  inextensible  fibers.  The  composite  material 
was  assumed  to  be  incompressible,  non-hardening  and  rigid-plastic.  The  yield  function 
was  formulated  in  terms  of  the  five  invariants  of  the  stress  tensor.  In  a  later  work  by 
Mulhem  et  al.  [68],  the  theory  was  extended  to  include  elastic-perfectly  plastic  matrix  and 
elastic  fibers. 

1.4  Choice  of  Constitutive  Theory 

A  broad  consensus  on  constitutive  modelling  for  composites  has  not  yet  emerged. 
Here,  the  merits  and  demerits  of  both  types  of  constitutive  theories  are  briefly  discussed, 
and  the  factors  influencing  the  choice  of  the  constitutive  theory  selected  in  this  study  are 
presented. 

The  primary  advantage  of  macromechanical  theories  is  that  they  require  less 
computational  effort  The  macromechanical  theories  are  well-enough  developed  to 
adequately  describe  nonlinear  response  to  mechanical  loadings.  Therefore,  structural 
analysis  of  composite  structures  based  upon  macromechanical  theories  is  a  viable  option. 
However,  much  work  needs  to  be  done  to  extend  these  theories  to  incorporate  the  effects 
of  complex  load  histories,  temperature  changes,  strain  rate  sensitivity  and  creep  etc. 

Micromechanical  theories,  besides  requiring  more  computational  effort,  require 
the  in  situ  properties  of  the  constituent  phases  which  have  to  be  determined  (back-calcu¬ 
lated)  out  from  tests  on  laminates.  The  advantage  associated  with  micromechanical 
theories  is  that  issues  like  complex  loading  histories,  thermal  effects,  and  creep  and  strain 
rate  sensitivity  can  be  addressed  with  less  difficulty.  The  idea  of  predicting  the  response  of 
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composites  starting  from  the  basic  constituents  is  intrinsically  appealing  since  the  effect  of 
using  different  fiber/matrix  combinations  can  be  studied. 

In  the  present  study  a  micromechanical  model  is  used  to  describe  the  behaviour  of 
fiber-reinforced  laminated  composites.  A  newly  proposed  micro-model  by  Pecknold  [57] 
is  investigated  and  further  developed.  Preliminary  results  from  this  model  are  very 
encouraging.  It  is  very  simple  in  description  and  in  the  course  of  this  study  was  found  to  be 
comparable  in  accuracy  to  Aboudi’s  Method  of  Cells. 

A  brief  account  of  previous  studies  employing  micromechanical  models  for  analy¬ 
sis  of  composite  structures  is  presented  in  the  next  section. 

1.5  Previous  Studies  Using  Micromechanics 
in  Structural  Analysis 

The  earliest  study  of  which  the  author  is  aware  which  employs  micromechanics  for 
general  structural  analysis  of  composites  is  by  Bahei-El-Din  et  al.  [43,44,45].  The  Vanish¬ 
ing  Fiber  Diameter  model  was  employed  in  a  three-dimensional  finite  element  analysis 
program  to  model  the  response  of  laminated  plates.  The  matrix  material  was  described  by 
thevonMises  yield  criterion  with  the  Prager-Ziegler  kinematic  hardening  rule.  Among 
the  problems  investigated  were  biaxial  loading  and  unloading  of  laminated  plates,  and 
loading  of  a  plate  with  a  hole;  stress  and  strain  distributions,  and  the  development  of 
plastic  zones  in  the  vicinity  of  the  hole.  The  results  showed  fair  agreement  with  experi¬ 
ments  and  with  detailed  finite  element  solutions.  The  strains  in  the  fiber  direction  were 
well  predicted;  however,  the  strains  in  the  transverse  direction  were  over-predicted, 
which  is  a  characteristic  of  the  Vanishing  Fiber  Diameter  micro-model. 

Adams  and  Crane  [37]  used  the  finite  element  method  to  study  the  response  of  a 
representative  volume  element  of  a  composite  lamina.  Nonlinear  behaviour  of  the  matrix 
was  taken  into  account,  and  special  procedures  were  developed  so  that  a  more  general 
state  of  stress  could  be  applied  to  the  representative  volume  element.  This  micromechan- 
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ical  modelling  procedure  was  employed  in  conjunction  with  Classical  Laminated  Plate 
Theory  to  determine  the  nonlinear  response  of  simple  laminates.  Their  micromechanical 
modelling  procedure  is  of  interest  because  it  is  one  of  the  first  attempts  at  incorporating 
micromechanics  into  structural  analysis  for  composites,  but  it  requires  extensive  computa¬ 
tional  resources,  and  its  implementation  for  structural  analysis  purposes  would  be  prohib¬ 
itively  expensive.  A  similar  procedure  was  used  by  Wu  et  al.  [46]  to  study  the  response  of 
metal  matrix  laminates  subjected  to  thermal  loading.  The  analysis  was  performed  by 
incorporating  the  Periodic  Hexagonal  Array  (PHA)  model  of  Dvorak  and  Teply  [47]  into  a 
general  purpose  finite  element  analysis  program  (ABAQUS  [48]).  The  PHA  model  is 
based  on  the  assumption  that  the  microstructure  of  a  composite  lamina  consists  of 
hexagonal  fibers  arranged  in  a  hexagonal  packing  geometry;  this  microstructure  can  be 
represented  by  a  triangular  representative  volume  element  The  effective  response  of  a 
lamina  was  obtained  by  finite  element  analysis  of  the  representative  volume  element.  The 
laminate  response  was  then  calculated  by  modelling  each  lamina  as  a  solid  element  It  was 
reported  by  Wu  et  al  [46]  that  the  laminate  analysis  required  sizable  computational 
resources,  since  obtaining  the  effective  lamina  response  itself  required  the  solution  of  a 
finite  element  problem. 

Rufin  et  al.  [49]  developed  a  computer  program  (MLAP)  which  employed  the 
Vanishing  Fiber  Diameter  model  of  Dvorak  and  Bahei-El-Din  [42,45]  in  conjunction 
with  Classical  Laminated  Plate  Theory  to  predict  the  response  of  metal  matrix  compos¬ 
ites.  They  observed  that,  in  general,  MLAP  tended  to  underestimate  the  elastic  moduli  of 
the  composites  and  overestimate  the  magnitude  of  plastic  strains.  These  errors  were  not 
significant  in  the  fiber  direction,  but  were  more  pronounced  in  the  transverse  direction. 
They  attributed  this  discrepancy  between  the  experimental  results  and  MLAP  results  to 
two  factors:  the  limitations  of  the  Vanishing  Fiber  Diameter  micro-model  and  the  nonlin¬ 
ear  nature  of  matrix  hardening,  which  could  not  be  captured  by  MLAP  since  it  assumed 
linear  matrix  hardening. 
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Arenburg  [70,71]  examined  the  behaviour  of  fiber-reinforced  metal  matrix  com* 
posites  using  Aboudi’s  Method  of  Cells.  The  nonlinear  matrix  behaviour  was  modelled  by 
the  unified  viscoplasticity  theory  of  Bodner  and  Partom[69].  The  effective  lamina  consti¬ 
tutive  relations  provided  by  the  Method  of  Cells  were  incorporated  into  a  First-order 
Shear  Deformation  Plate  Theory.  Results  were  presented  for  tensile  response  of  notched 
coupons,  plates  in  bending,  and  moment-curvature  responses  of  laminates.  The  numeri¬ 
cal  results  showed  good  agreement  with  experimental  results. 

1.6  Outline  of  the  Analysis  Procedure 

Fig.  1.6  shows  a  laminated  plate  structure  discretized  by  finite  elements.  Atypical 
element  within  the  structure  is  also  shown.  Stiffnesses,  strains  and  stresses  are  tracked  at 
sampling  points  (material  points)  within  each  element.  This  information  is  provided  by 
the  material  model,  which  interfaces  with  the  main  analysis  program.  The  interface  is 
called  a  standard  interface  since  the  laminated  heterogeneous  nature  of  the  structure  is 
hidden  from  the  main  analysis  program,  and  the  information  being  transferred  between 
the  material  model  and  the  main  analysis  program  is  essentially  the  same  as  would  be 
required  for  the  analysis  of  a  structure  made  from  a  homogeneous  material. 

Obviously,  the  standard  procedures  for  sampling  die  stiffness  and  stress  are  not 
applicable  since  the  material  in  the  neighborhood  of  a  material  point  is  heterogeneous.  This 
problem  can  be  solved  if  the  material  in  the  neighborhood  of  the  material  point  is 
represented  by  an  equivalent  homogeneous,  possibly  anisotropic,  material.  In  this  study,  a 
material  modelling  strategy  is  presented  which  is  based  on  this  homogenization  procedure. 

As  shown  in  Fig.  1.6,  the  proposed  material  model  has  two  components: 

(1)  The  Micro-model,  which  describes  the  response  of  a  unidirectional  lamina. 

The  stress  and  strain  fields  within  the  fiber  and  matrix  phases  are  combined 
appropriately  to  yield  effective  stresses  and  strains  in  the  lamina.  Similarly, 
the  effective  lamina  constitutive  relations  are  developed  from  the  constitu- 
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tive  relations  of  fiber  and  matrix.  The  three-dimensional  nature  of  the 
stress-field  in  the  constituent  phases  is  accounted  for,  to  better  represent 
the  fiber  Viatrix  interaction  phenomenon. 

(2)  The  Sublaminate  model,  which  generates  the  response  of  a  group  of  lami¬ 
nae,  by  combining  the  stress  and  strain  fields  within  each  lamina  to  yield 
overall  stresses  and  strains.  The  constitutive  properties  of  a  group  of  lami¬ 
nae  are  combined  in  an  appropriate  manner  to  yield  the  constitutive  prop¬ 
erties  of  an  effective  homogeneous  material,  which  represents  the  group 
of  laminae.  A  three-dimer  sional  lamination  scheme  taking  into  account 
the  in-plane  and  out-of-plane  stresses  has  been  employed  in  the  present 
study,  since  the  out-of-plane  stresses  play  a  decisive  role  in  initiation  of  de¬ 
lamination,  and  should  be  accounted  for. 

After  the  stresses,  strains  and  constitutive  properties  of  the  effective  homogeneous 
material  representing  the  group  of  laminae  has  been  determined,  this  information  is 
passed  to  the  structural  analysis  package.  This  process  of  homogenization  of  fields  within 
a  fiber  and  matrix  to  obtain  the  response  of  a  lamina,  followed  by  homogenization  of  the 
various  fields  within  a  group  of  laminae  to  determine  the  response  of  a  ply  cluster 
(sublaminate),  may  be  referred  to  as  the  composition  process. 

At  the  end  of  the  composition  process,  we  have  the  stresses,  strains  and  effective 
properties  of  an  equivalent  material  representing  a  sublaminate.  This  information  is 
passed  to  the  structural  analysis  program  through  a  standard  interface.  This  is  especially 
advantageous,  since  it  makes  available  a  large  number  of  the  elements  contained  in  the 
element-libraries  of  most  structural  analysis  packages.  This  provides  the  capability  to 
analyze  laminated  structures  without  the  investment  of  the  considerable  effort  required  to 
develop  special  purpose  elements. 
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During  loading  of  the  structure,  changes  in  strains  occur  at  the  material  points. 
This  change  is  communicated  to  the  material  model  through  a  standard  interface.  The 
effect  of  this  change  is  propagated  down,  first  to  the  sublaminate  model,  and  then  to  the 
micro-model,  where  the  stresses,  strains,  and  the  properties  of  the  fiber  and  the  matrix 
may  be  updated.  This  process  of  propagating  the  effect  of  change  at  the  material  point 
down  to  the  basic  constituents  may  be  referred  to  as  the  decomposition  process. 

Within  this  general  scheme,  any  type  of  micro-model  and  sublaminate  model  can 
be  used.  In  this  study  the  micromechanical  model  proposed  by  Pecknold  [57]  is  employed 
because  of  its  simplicity.  As  a  part  of  this  study  the  micro-model  was  investigated  and  was 
found  to  be  comparable  in  accuracy  to  the  refined  micromechanical  models,  such  as  the 
one  presented  by  Aboudi  [7,51,52,53,54]. 

In  the  following  chapter  the  development  of  the  micro-model  is  presented. 
Examples  are  presented  to  demonstrate  its  accuracy.  Chapter  Three  presents  the  subla¬ 
minate  model.  In  the  fourth  chapter,  the  procedure  adopted  for  updating  the  stresses  at  a 
material  point  is  presented.  In  the  fifth  chapter,  numerical  examples  are  given;  results  are 
presented  for  predictions  of  nonlinear  response  of  laminae  and  laminates  subjected  to 
uniform  loadings,  strength  response  of  laminae  and  their  initial  yield  surfaces  are 
presented.  The  results  from  the  proposed  material  model  are  compared  to  predictions 
from  other  micromechanical  models.  Chapter  Six  presents  structural  level  applications, 
i.e.  examples  of  laminated  structures  which  were  solved  by  the  material  modelling  proce¬ 
dure  developed  in  this  study.  Chapter  Seven  is  devoted  to  conclusions  from  the  present 
study  and  recommendations  for  future  work. 
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(b)  Fiber  and  matrix  in  series  (Reuss)  arrangement. 


Fig.  1.1  Mechanics  of  Material  models 
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Fig.  1.2  Description  of  the  Two-phase  model 
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Fig.  1.3  Description  of  the  Three-phase  model 
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Fig.  1.4  Composite  Cylinder  Assemblage  (CCA)  model 
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Fig.  1.5  Assumed  microstructure  and  representative  volume 
element  In  Aboudl's  Method  of  Cells. 
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CHAPTER  2 


MICRO-MODEL 


In  this  chapter  the  first  component  of  the  material  model,  i.e.  the  micro-model, 


is 


described.  The  micro-model  describes  the  response  of  a  unidirectional  lamina,  starting 
from  the  fiber  and  the  matrix  constitutive  descriptions.  The  micro-model  employed  here 
is  the  one  proposed  by  Pecknold  [57].  The  assumptions  made  in  order  to  obtain  the 
effective  stiffness  of  the  lamina  are  stated,  and  the  important  features  of  the  model  are 
high-lighted.  Constitutive  relations  used  to  model  the  nonlinear  behaviour  of  the  matrix 
material  are  presented.  Micro-failure  criteria  for  determination  of  various  failure  modes 
and  procedure  for  incorporating  the  effect  of  damage  are  also  presented.  Examples  are 
presented  to  demonstrate  the  accuracy  of  the  micro-model  by  comparing  its  predictions 
of  effective  elastic  moduli  to  results  obtained  from  detailed  numerical  procedures  and 
other  established  methods. 


2.1  Description  of  the  Micro-model 


The  proposed  micromechanical  model  is  based  on  the  assumption  that  die  internal 
microstructure  of  the  fiber-reinforced  lamina  consists  of  square  fibers  arranged  in  a 
doubly  periodic  array,  as  shown  in  Fig.  2.1.  Following  Shaffer’s  [12]  approach,  it  is 
assumed  that  the  stress  acting  parallel  to  the  *3  axis  is  taken  up  by  two  strips  of 
material:  one  is  a  strip  of  pure  matrix  enclosed  by  lines/!  and  mm,  and  the  other  is  a  strip 
of  material  enclosed  by  lines  mm  and  nn.  The  material  in  the  second  strip  consists  of 
fiber  and  matrix  phases  acting  in  series.  Together  the  two  strips  are  assumed  to  be 
acting  in  parallel  to  withstand  the  transverse  loading.  For  axial  loading  ( loading  in  the 
fiber  direction),  both  the  fiber  and  the  matrix  phases  are  assumed  to  be  acting  in 
parallel. 


24 


Also  shown  in  Fig.  2.1  is  the  representative  unit  cell  which  corresponds  to  the 
assumed  geometry  of  the  idealized  composite.  The  unit  cell  is  further  subdivided  into 
three  subcells;  one  fiber  subcell,  denoted  by  f,  and  two  matrix  subcells,  denoted  by  mA 
and  n>B  respectively.  The  dimensions  of  the  unit  cell  are  1  unit  square,  hence  the  name 
unit  cell.  The  dimensions  of  fiber  and  matrix  subcells  are  given  in  terms  of  quanti¬ 
ties  Wt  and  Wm,  which  are  defined  as 

wt  =  Jvt  ,  wm  =  l  -  wt  (2.1) 

where  Vt  =  Fiber  Volume  Fraction. 

In  the  following  section  the  assumptions  made  in  homogenizing  the  stresses  and 
strains  will  be  presented,  however  before  that,  for  the  sake  of  convenience,  the  three 
subcells  are  grouped  to  define  two  elements:  Material  Element  A,  consisting  of  the  fiber 
subcell  f  and  die  series-or-parallel  connected  matrix  subcell  mA,  and  Material  Element 
B,  consisting  of  the  remaining  matrix  subcell  *%• 

2.2  Stress  and  Strain  Homogenization  in  Unit  Ceil 

The  effective  stresses  and  strains  in  the  lamina  are  determined  from  the  subcell 
values  in  two  stages:  first,  fiber  subcell  f  and  matrix  subcell  mA  are  used  to  construct 
Material  Element  A;  then  Material  Element  A  and  Material  Element  B  are  used  to  con¬ 
struct  the  unidirectional  lamina.  The  procedure  that  is  used  is  shown  in  Fig.  2.2  by  means 
of  spring  analogs. 

The  stress  and  strain  homogenization  procedure  is  now  set  down  more  formally  in  the 
next  section. 

2.2.1  Homogenized  Stresses  and  Strains  in  Material  Element  A 

Within  Material  Element  A,  the  axial  (11)  components  are  treated  differently 
from  the  remaining  components.  In  order  to  express  the  equilibrium  and  compatibility 
conditions  compactly,  the  stress  and  strain  vectors  are  partitioned  as  follows: 
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Thus  Ol  and  are  (lxl),  and  Ot  and  €t  are  (5x1)  matrices.  Note  that  the  components 
of  the  vectors  are  not  listed  in  conventional  order.  Although  this  particular  ordering  of 
components  is  not  necessary  for  micro-model  computations,  it  is  most  convenient  for 
subsequent  computations  at  the  sublaminate  model  level. 

It  is  assumed  that  the  axial  (fiber)  direction  strain  (€l)  are  the  same  in  the  fiber  and 
the  matrix,  and  the  remaining  five  “transverse”  stress  components  (07)  are  the  same  in 
the  fiber  and  the  matrix.  These  relations  are  expressed  compactly  as 

ft. ■  ft,  ■  ft 

A 

The  complementary  stress  and  strain  components  (i.e.  oL  and  eT  )  in  Material  Element  A 

are  weighted  averages  of  these  stress  and  strain  components  in  the  fiber  and  the  matrix, 
i.e. 

ft-  "ft  *  '-ft..  “ 

2.2.2  Homogenized  Stresses  and  Strains  in  Unit  Cell 

Material  Elements  A  and  B  are  connected  in  parallel  for  all  component  direc¬ 
tions,  i.e.  their  strains  are  the  same  and  the  unit-cell-average  stresses  are  weighted 
averages  of  the  stresses  in  Elements  A  and  B,  i.e.. 
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H-  (I-  W. 

W  •  *  '-W 


(2.4a) 

(2.4b) 


in  which  the  stresses  and  strains  are  (6x1)  vectors,  the  subscript  C  stands  for  the 
unit-cell  averages. 


Eqs.  (2.3)  and  (2.4)  are  the  homogenization  relationships  that  completely  define 
the  micromechanical  model.  These  relationships  are  also  valid  for  incremental  stresses 
and  strains.  The  development  of  the  tangent  stiffness  of  the  unit  cell  follows  as  a 
consequence  of  assumptions  made  in  these  relations.  The  tangent  stiffness  of  the  unit 
cell  is  also  the  tangent  stiffness  of  the  lamina  since  the  unit  cell  represents  die  idealized 
composite  lamina. 

2.3  Tangent  Stiffness 

The  tangent  stiffness  relates  incremental  (unit-cell-average)  stresses  and  strains. 
The  tangent  stiffness  of  the  unit  cell  is  built  up  from  tangent  stiffness  relations  for  the 
two  constituent  phases;  these  may  be  linear  relationships  or  they  may  reflect  various 
types  of  nonlinearities. 


The  incremental  stresses  and  strains  within  a  material  are  related  to  each  other  by 
relations  of  the  form 
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(2.5) 
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where 


ax,  a2,  a3  =  Coefficients  of  thermal  expansion  in  xv  and  x3 
directions  respectively. 

A7  =  Change  in  Ibmperature. 

The  above  equation  can  be  compactly  written  in  partitioned  form  in  accordance  with  Eq. 
(22)  as 


The  tangent  compliance  matrix  is  now  partially  inverted  to  provide  relations  of  the  form 


(2.7a) 


in  which  P«  is  (lxl),  P u  is  (1x5),  and  P«  is  (5x5),  and  they  are  given  by: 

p it  -  (1/S u) 

Pit  =  -  (l/S//)  •  S„  (2.7b) 

P«  =  S„-(1/Sh)  •  S  ti  S  it 

PaL  and  Par  define  the  vector  of  partially  inverted  thermal  compliances,  and  are  given 
by 


PaL  =  -  (l/S a)  aL 
Pay  *  %  loj,  Stf 


(2.7c) 


2.3.1  Material  Element  A 

With  the  tangent  relations  for  fiber  subcell  f  and  matrix  subcells  mA  expressed  in 
the  form  of  Eq.  (2.7a),  the  homogenization  relation  Eqs.  (2.3)  can  be  applied  to  give  the 
partially-inverted  tangent  compliance  for  Material  Element  A  as 
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>  r*  *  i  ^ 

daL  I  P/<  deL  P°l 

<fcr  -P,H  P„  daT  p«r 

Ja  l  -Ia'*  jA 


(2.8a) 


where 


[P ii  |  Pit  1  f  P«  j  P«  I  f  P«  j  Pu  ‘ 

- 1 - I  =  W,  - 1 - I  +  wm  ---i—- 

-•£!*«  I  *p»l  p«  J  L -p»|  pn _ 

A  i  m* 

(2.8b) 

=  „  Li“4  +  *  1-4 
wA  w,  ■  1  < 

The  partial  inversion  is  now  completed  to  give  the  tangent  stiffness  of  Material 


Element  A  as 


d°L  1  J  Cb  J  c„  deL  Pi 

■ - >  ss  I - 1 -  - - ■+  AT-  ---  ■ 

d°T  C«  j  C  „  d€r  Pt 
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(2.9a) 


in  which 


C «  =  P n 


C u  -  Pit  C 'ft 
C u  -  P u  +  C u  Pjt 

Pt~  -  P  tt  Par 

Pi  =  PaL  +  P It  Pt 

Eq.(2.9a)  can  be  written  compactly  as 


(2.9b) 


M  =  [C].H+ArW. 


(2.10) 


2.3.2  Material  Element  B 

The  tangent  compliance  relations  for  Material  Element  B,  which  consists  of  a 
single  matrix  subcell,  is  inverted  to  give  the  tangent  stiffness  relation 
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H=[c]bH+4tH 


(2.11) 


2.3.3  Unit  Cell 


The  stress  and  strain  homogenization  relations  (Eqs.  (2.4))  are  now  applied  to 

give  the  tangent  stiffness  relation  for  the  unit  cell  as 

Wc=  [c]cMc+*rWc 

where 

(2.12a) 

[  c]  =  Tangent  Stiffness  of  Unit  Cell 

[c]c=  w'[c\+  w-[c\ 

(2.12b) 

{p}  =  wf{p)  +  ^  M 

*  JC  1  JA  1  JB 

(2.12c) 

It  should  be  mentioned  here  that  the  assumed  microstructure  of  the  composite,  as 
shown  in  Fig.  2.1,  corresponds  to  a  material  possessing  tetragonal  symmetry,  the 
stiffness  matrix  of  which  is  described  by  six  independent  moduli.  However,  because  of 
the  simplifying  nature  of  the  homogenization  relations  assumed  here  the  stiffness  matrix 
given  by  Eq.  (2.12a)  corresponds  to  that  of  a  transversely  isotropic  material,  which  is 
characterized  by  five  independent  moduli.  This  is  in  fact  advantageous,  since  it  is  known 
that  a  unidirectionally  reinforced  lamina  tends  to  be  transversely  isotropic.  In  the  next 
section  the  procedure  for  determining  equivalent  effective  moduli  of  the  composite 
lamina  is  illustrated.  The  information  regarding  the  effective  moduli  of  the  lamina  is  not 
needed  for  the  analysis  of  the  lamina  or  the  laminate;  this  information  is  presented  only 
to  facilitate  comparison  of  the  predicted  moduli  of  composite  laminae  with  other  micro¬ 
mechanical  models. 
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2.4  Elastic  Moduli  and  Coefficients 
of  Thermal  Expansion 

With  the  stiffness  matrix  of  the  composite  lamina  determined  from  Eq.  (2.12a), 
the  effective  elastic  moduli  of  the  composite  can  then  be  calculated  by  inversion  of  the 
stiffness  matrix  or  with  the  help  of  following  useful  relations,  which  are  valid  for 
transversely  isotropic  materials.  If  Cn,  C12.  C23  etc.  are  the  terms  of  the  stiffness  matrix 
of  the  composite,  then 


Kjs  “  "2  (  Cn  +  Cu ) 
Gj3  =  —  (  C22  -  CM ) 

Gu  ”  ™  Gi  •*  Cjj 
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(2.13a-f) 
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(2.13g-h) 


K23,  G73,  J?2  and  V23  are  the  plane  strain  bulk  modulus,  shear  modulus,  transverse 
Young’s  modulus  and  Poisson’s  ratio  in  the  *2 -*3  plane,  respectively.  Eiand  G\  are 
the  Young’s  modulus  and  the  shear  modulus,  in  the  direction  parallel  to  the  fibers  and  v\ 
is  the  Poisson’s  ratio  when  a  uniaxial  stress  is  applied  in  the  direction  of  the  fibers. 


Coefficients  of  thermal  expansion  for  the  composite  lamina  can  be  extracted  by 
inverting  Eq.  (2.12a),  yielding 

Wc=  [s]cWc  +  iTMc  (214a) 
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where 


[  s]  =  Tangent  Compliance  of  Composite  =  [  c] 


-l 


'a  ' 


<h 

0 


1-L  - 


=  Coefficients  of  Thermal  expansion 
of  Composite  lamina. 


(2.14b) 


and  a.y,  a2  and  a3  are  coefficients  of  thermal  expansion  in  xy  X2  and  x$  directions  re- 
spectively. 


Having  described  the  procedure  for  obtaining  effective  lamina  stiffnesses  and 
compliances,  attention  is  now  focused  on  the  procedure  for  incorporating  and  modelling 
various  types  of  nonlinearities  in  the  lamina  response.  It  is  often  assumed  that  the  fibers 
within  a  lamina  behave  linearly  up  to  the  point  of  fracture,  and  that  the  observed  nonlin¬ 
earity  in  lamina  response  is  caused  by  the  nonlinear  response  of  the  matrix  material.  With 
this  in  mind,  constitutive  descriptions  for  the  matrix  material  are  presented  in  the  next 
section. 

2.5  Constitutive  Descriptions  of  Matrix  Nonlinearity 

Nonlinear  matrix  behaviour  can  be  an  important  factor,  since  a  degrading  shear 
modulus  may  result  in  premature  buckling  of  composite  shells.  Furthermore,  matrix 
softening  in  shear  may  interact  with  or  initiate  kink-banding.  It  has  also  been  observed 
that  the  nonlinear  response  of  thick  laminates  is  essentially  reversible  over  a  substantial 
loading  range  (Camponeschi  [72,73]).  Keeping  these  considerations  in  view,  a  nonlinear 
elastic  description  for  the  matrix  material  is  proposed  as  a  reasonable  first  step  towards 
describing  the  nonlinear  response  of  composite  laminae.  A  Ramberg-Osgood  type 
stress-strain  relation  is  proposed  here,  and  is  described  in  the  next  section. 
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2.5.1  Ramberg-Osgood  Relations 


A  Ramberg-Osgood  type  stress-strain  relation  which  is  valid  for  general  multi- 
axial  stress  states  may  be  expressed  as 


C‘J  = 


+ 


.  b 
P>“  *  1G 


(2.15) 


The  first  term  on  the  right  hand  side  is  the  elastic  strain,  the  second  term  corresponds  to 
nonlinear  pressure-volume  change  and  the  third  term  corresponds  to  the  nonlinear  shear 
softening  effects.  In  the  above  expression  K  and  G  are  elastic  bulk  and  shear  moduli. 


a,  b,  po,  to,  n  and  m  are  material  constants  to  be  determined  and 


P  —  'y0**  Sfj  —  Oy  -  p  6(j 


(2.16) 


The  quantity p  is  recognizable  as  one-third  of  the  first  invariant  of  the  stress  tensor,  and 
xe  is  the  second  invariant  of  the  deviatoric  stress  tensor. 


The  incremental  form  of  Eq.  (2.15)  is 


Cl  ( D  h  /  r 

'  **  *  W  "  (*J  dp  +  20  (n)  *• 

b  n- 1  [x ,V_1 
+  2  G  2  \to  j  T(t 


?tt  dsu 


where  s,y  are  the  normalized  deviator  stresses  defined  as. 


f  _  £</ 


Eqs.  (2.15)  and  (2.17)  can  be  expressed  in  matrix  form  as 

=  £  ci  [Si]  +  c3  [S2J  J  jcrj 


(2.17) 


(2.18) 


(2.19) 


jdej  =  [  rfi  [S,]  +  d2  [S2]  +  d ,  [Sj]  J  jdaj  (220) 
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The  the  stress  and  strain  vectors  are  ordered  as  in  Eq.  (2.2),  and  the  matrices 
Si,  s2,  Sj  are: 


[Si]  = 


1 

1 

l  !  Symm 

000 

,  [s2]  = 

0  !  Symm 

0  0  2 

1110 

0  0  0  1 

00000 

0  0  0  0  2 

0  0  0  0  0  0_ 

0  0  0  0  0  2m 

[S3]  =  S  s7 


ST  =  *< 


S22 

23T,2 


$33 

2T23 

2Ji3 

V 


(2.21) 


and  the  functions  Ci,  C2  and  d\,  di,  d$  are  given  by 


c.  -  [^  ~  (1  ♦(>,)] 

C  (1*»J 

*' ■  (ltm  “■)  -is  ,1+*‘)] 


di  =  Ci 


(2.22a-g) 


The  standard  Ramberg-Osgood  relations  as  described  above  can  be  used  to  de¬ 
scribe  the  nonlinear  response  of  most  metal  matrix  composites.  However,  in  their  stan¬ 
dard  form  the  Ramberg-Osgood  relations  also  result  in  significant  nonlinearity  in  re¬ 
sponse  to  transverse  loading.  The  limited  experimental  data  available  for  polymer  matrix 
composites  seem  to  suggest  that  even  though  the  lamina  response  is  quite  nonlinear  in 
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shear,  no  significant  nonlinearity  exists  in  response  to  transverse  normal  loading.  The 
above-mentioned  relations  for  matrix  behaviour  cannot  accommodate  both  these  beha¬ 
viours  simultaneously.  Therefore,  modifications  are  proposed  to  the  above  proposed 
Ramberg-Osgood  relations,  to  be  used  in  the  case  when  nonlinearity  exists  only  in  shear 
response.  The  modified  constitutive  relations  for  the  matrix  material  are  presented  in  the 
next  section. 

2.5.2  Modified  Ramberg-Osgood  Relations 

The  suggested  modifications  are  as  follows.  Since  nonlinear  behaviour  only  in 
shear  is  sought,  the  pressure-volume  change  term  can  be  dropped  from  consideration, 
and  the  response  to  normal  loading  in  all  three  directions  can  Vi  assumed  to  be  linear. 
The  detailed  description  of  the  modifications  made  to  the  Ramberg-Osgood  relations  is 
as  follows. 

It  is  assumed  that  the  response  to  On,  On,  and  033  is  essentially  linear,  i.e., 

€a  =  ej  ,  for  i  -  /  (2.23a) 

The  response  in  shear  is  described  as 

•  for  ■"  <2_23b) 

The  definition  of  *e  is  changed  to 

for  H ;  <224> 

The  incremental  strains  are  given  by 

de«  -  de*  ,  for  i  -  j  (2.25a) 

and 

*»  -  *  £  (*)""  *« 

♦  £  r»  J,  *«  ,  for  i  .  i  (2.25b) 
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where  the  normalized  deviator  stresses  are  defined  in  similar  fashion  as  before,  but  with 
respect  to  the  new  definition  of  xe  given  in  Eq.  (2.24). 


Eqs.  (2.23)  and  (225)  can  be  expressed  in  matrix  form  as: 

JeJ  =  [  c,  [SO  +  c2  IS*]  J  JcrJ 


(2.26) 


M=  [ci 


[Sil  +  c2  [S2]  + 


where  the  matrices  Si,  s2,  s3  are  defined  as: 


d3  [S3]  J  jdaj 


(227) 


[Sx]  = 


Symm 
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1110 
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000000 


,  IS*]  = 


l  Symm 

0  2(1  +bi) 

0  0  1 

0  0  0  2(1 +bi) 

0  0  0  0  2(1  +  b{) 


[S3]  =  S  ,  ff  S 


and  the  functions  ct,  c2  and  d3  are  given  by 


-■(w-M 


(2.28) 


,  _  1  n-1 

di-i G~rb' 


(2.29a-d) 


It  should  be  noted  that  because  of  the  modifications  made  in  the  definition  of  xt  the 

modified  Ramberg-Osgood  relations  are  no  longer  based  on  invariants  of  stress,  and  are 
applicable  only  in  lamina  coordinates. 


2.6  Failure  Criteria 


Numerous  macromechanical  criteria  at  the  lamina  level  have  been  proposed, 
which  view  the  lamina  as  an  effective  homogeneous  anisotropic  material,  and  describe  the 
failure  surface  as  a  stress  polynomial,  e.g.  the  “'ftai-Wu  Tfensor  Polynomial”  criterion 
[88],  and  the  “Quadratic  Stress  Polynomial”  criterion  proposed  by  Hashin  [89].  A  com¬ 
prehensive  survey  of  macromechanical  failure  criteria  can  be  found  in  the  works  of  TSai 
[90]  and  Nahas  [91]. 


Aboudi  [52]  employed  the  Method  of  Cells  to  predict  the  off-axis  strength  of 
unidirectional  composite  laminae.  Failure  of  the  composite  lamina  was  predicted  by 
using  micro-failure  criteria  for  the  fiber  and  matrix  phases.  The  lamina  was  subjected  to 
off-axis  loading,  and  stresses  within  the  various  subcells  of  the  unit  cell  representing  the 
lamina  were  monitored.  A  maximum  stress  failure  criterion  was  applied  to  the  constituent 
phases,  to  define  the  failure  of  the  lamina.  Aboudi’s  analytical  results  were  in  good 
agreement  with  experimental  results. 

Here,  an  attempt  is  made  to  employ  micromechanical  failure  criteria  to  detect 
damage  and  failure  in  the  lamina,  since  the  stresses  and  strains  in  the  constituent  fiber  and 
matrix  phases  are  readily  available  through  the  micro-model.  Attention  is  focused  on 
in-plane  failure  mechanisms  only.  A  brief  description  of  these  damage  mechanisms  and 
failure  modes  along  with  the  criteria  defining  their  initiation  are  presented  in  this  section. 
2.6.1  Fiber  Fracture  in  Tension 

A  unidirectional  lamina  subjected  to  tensile  loading  in  the  fiber  direction  may  fail 
by  fiber  fracture.  If  denotes  the  tensile  strength  of  the  fiber,  then  the  fiber  fracture 
mode  of  failure  is  determined  by 


» 

*11 


Fiber  Fracture 
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(230) 


2.6.2  Matrix  Shearing 


If  5(m)  is  the  shear  strength  of  the  matrix,  then  the  matrix  shearing  failure  criterion  is 


max(  \*n\  >|  t^I  )  -  5<m)  =£>  Matrix  Shearing^^ 

2.6.3  Matrix  Cracking  in  Transverse  Tension 


W  (2.31) 


The  criterion  determining  matrix  cracking  in  transverse  tension  is  taken  as 


Matrix  Cracking  in 
Transverse  Tension 


(232) 


where  is  the  tensile  strength  of  the  matrix  material. 


2.6.4  Matrix  Crushing  in  Transverse  Compression 

The  failure  of  a  lamina  in  transverse  compression  is  assumed  to  be  a  result  of 
compressive  failure  of  the  matrix  material.  If  and  are  compressive  stresses  in 

the  matrix  subcells,  and  denotes  the  compressive  strength  of  the  matrix  material, 

then  the  compressive  failure  of  the  lamina  in  transverse  direction  may  be  determined  by 
the  following  criterion 


max 


Matrix  Crushing  in 
Transverse  Compression 


(233) 


2.6.5  Matrix  Cracking  in  Axial  Tension 


One  possible  damage  mechanism  in  fiber-reinforced  composites  is  cracking  of  the 
matrix  due  to  axial  tension.  Although  not  envisioned  as  a  very  serious  damage  mode,  nor  a 
very  commonly  occurring  one,  its  consideration  may  be  warranted  in  composites  with  low 
fiber  volume  fraction.  The  criterion  indicating  axial  direction  matrix  cracking  may  be 
taken  as 
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Matrix  Cracking  in 
Axial  Tension 


(234) 


2.6.6  Lamina  Failure  by  Kink-banding 


In  axial  compression,  die  failure  is  assumed  to  be  triggered  by  fiber  micro-buck¬ 
ling  resulting  in  kink-banding.  At  this  time  (1992)  no  adequate  analytical  formulations 
exist  which  can  predict  with  consistent  accuracy  the  compressive  strength  of  unidirectional 
composites  in  compression.  A  number  of  micromechanics-based  analytical  formulae 
have  been  proposed,  e.g.  Rosen  [92],  Argon  [93],  Budiansky  [94],  and  Hahn  and  Williams 
[95],  but  they  do  not  correlate  well  with  the  experimental  results.  Although  the  various 
analytical  models  for  prediction  of  compressive  strength  have  not  shown  good  quantita¬ 
tive  agreement  with  experimental  results,  they  have  served  to  identify  the  important  pa¬ 
rameters  influencing  the  compressive  strength  of  composites.  Based  on  these  works  the 
important  parameters  influencing  the  compressive  strength  of  unidirectional  composites 
appear  to  be  matrix  shear  modulus,  fiber  misalignment,  fiber  diameter,  fiber  volume 
fraction,  and  fiber/matrix  interface  strength. 


Here,  a  modified  form  of  die  analytical  formula  suggested  by  Hahn  and  Williams 
[95]  is  employed  for  estimating  the  compressive  strength  of  the  lamina.  The  compressive 
strength  of  a  unidirectional  lamina  with  no  defects  or  fiber  misalignment  is  given  by  Hahn 
and  Williams  as: 


“  vt  Gn  (235) 

where,  Vt  is  the  fiber  volume  fraction  and  Gn  is  the  in-plane  shear  modulus  of  die  lamina. 
If  die  lamina  stress-strain  relations  are  nonlinear,  then,  the  tangent  shear  modulus  Gbr 
should  be  used  in  the  above  relation  and  die  critical  compressive  stress  is  given  as 
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(236) 


In  the  present  analysis,  micro-buckling  resulting  in  kink-banding  failure  is  assumed  to 
occur  if  the  axial  compressive  stress  on  the  lamina  exceeds  the  critical  micro-buckling 
stress  given  by  Eq.  (236),  i.e.. 


Micro-buckling  Kink-banding 


(237) 


2.7  Damage  Modelling 

Due  to  the  aforementioned  failure  modes,  some  degree  of  stiffness  degradation  and 
stress  relief  takes  place  within  the  lamina.  The  assumptions  made  about  lamina  stiffness 
degradation  and  stress  relief  are  discussed  here. 

If  tensile  fiber  fracture  mode  is  detected  then,  the  fiber  tangent  compliance  matrix 
is  modified  as  follows: 


t 


where  Dfn  is  the  stiffness  reduction  factor  which  serves  to  reduce  the  fiber  axial  stiffness 

and  is  arbitrarily  assumed.  Since  the  tensile  fiber  fracture  failure  corresponds  to  the 
opening  of  cracks,  the  axial  stress  in  the  fiber  is  also  relieved. 


(239) 
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In  the  case  of  matrix  cracking  in  the  axial  direction,  the  matrix  tangent  compliance 
matrix  is  modified  as  follows 

t 

^1  sn  0  0 

Sn 

Symm 

where  Dmn  is  the  stiffness  reduction  factor  which  serves  to  reduce  the  axial  stiffness  of  the 

matrix.  Since  axial  tensile  matrix  cracking  corresponds  to  opening  of  cracks,  die  axial 
stress  in  the  matrix  subcells  is  relieved. 
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(2.41) 


Transverse  matrix  cracking  is  a  common  phenomenon  in  fiber-reinforced  compos¬ 
ites  and  it  can  cause  significant  stiffness  degradation  or  even  complete  failure  in  lami¬ 
nates.  Transverse  cracking  can  occur  due  to  transverse  tensile  stresses  or  in-plane  shear 
stresses,  or  perhaps  due  to  the  combined  effect  of  the  two  stresses.  If  the  lamina  fails  in  the 
transverse  direction  in  tensile  matrix  cracking  mode,  or  a  matrix  shearing  mode,  then  the 
matrix  tangent  compliance  matrix  is  modified  as  follows: 
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where  and  D-44  are  stiffness  reduction  factors  which  reduce  the  stiffness  of  the 

matrix  material  in  the  transverse  direction  and  in  in-plane  shear  respectively.  If  the 
failure  mode  is  transverse  matrix  cracking  due  to  tension,  then  since  it  corresponds  to 
opening  of  cracks  the  transverse  direction  stress  in  the  matrix  is  also  relieved. 


.<**>  „(■*) 


'22  ’  °22  =  0  (2.43) 

Since  the  stresses  in  matrix  subcell  mx  are  related  to  stresses  in  the  fiber  via  the  constraint 


Eq.  (2.3a)  ,  the  cr®  stress  in  the  fiber  is  also  reduced  to  zero. 


<4°  -  0  (2.44) 

The  mechanics  of  lamina  compressive  failure  in  the  transverse  direction  are  not 
sufficiently  clear  as  yet,  and  there  exists  a  paucity  of  experimental  data  regarding  the 
behaviour  of  unidirectional  composites  loaded  in  transverse  compression.  If  compression 
failure  of  matrix  is  indicated  by  Eq.  (233),  then  the  matrix  tangent  stiffness  is  reduced  in 
accordance  with  Eq.  (2.42);  but,  since  no  opening  of  cracks  is  envisioned  in  this  case,  no 
stresses  are  relieved  in  the  components  of  the  unit  cell. 

Failure  in  a  lamina  can  also  take  place  due  to  compressive  failure  of  the  fiber  or 
kink-banding.  This  is  considered  as  a  catastrophic  failure,  and  the  tangent  stiffnesses  of 
the  constituent  phases  are  reduced  according  to  Eqs.  (238-2.40,2.42);  but  no  stress  reliev¬ 
ing  is  carried  out 

In  the  next  section  numerical  examples  will  be  presented,  and  the  predictions  of  the 
elastic  moduli  of  composites  based  on  the  micro-model  are  compared  to  predictions  from 
other  established  procedures. 

2.8  Results  for  Elastic  Moduli 

The  micro-model  is  the  most  important  component  of  the  material  modelling 
procedure,  and  may  be  considered  the  kernel  upon  which  the  accuracy  of  the  overall 
scheme  depends.  Therefore  it  is  necessary  that  the  accuracy  of  the  micro-model  be  ex- 
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amined  in  detail.  Hie  first  step  is  examination  of  the  accuracy  of  the  elastic  moduli 
predicted  by  the  model.  With  this  in  mind,  the  results  from  the  micro-model  are  com¬ 
pared  with  a  number  of  procedures,  such  as  the  Composite  Cylinder  Assemblage  (CCA) 
model  of  Hashin  and  Rosen  [9,10,25,26],  numerically  obtained  results  of  Chen  and 
Cheng  [33]  and  Pickett  [34],  and  predictions  from  Aboudi’s  Method  of  Cells  [7,54]. 

2.8.1  Comparison  with  CCA  Model 

The  CCA  model  model  of  Hashin  and  Rosen  gives  closed  form  expressions  for 
Ei  >Kzs  ,Gi  and  whereas  closed  form  expressions  for  G23  >£2  and  V23  could  not 
be  obtained;  these  quantities  could  only  be  bounded. 

As  mentioned  in  the  previous  chapter,  the  Three-phase  model  of  Hermans  [22] 
and  Christensen  and  Lo  [23]  is  related  to  the  CCA  model.  Again,  it  was  noticed  by 
Hashin  [9]  that  the  expressions  for  E\  ,Kb  ,Gi  and  v\  obtained  by  Hermans  are 
exactly  the  same  as  given  by  the  CCA  model,  whereas  Hermans’  derivation  for  the 
expression  for  G23  contained  some  errors.  These  errors  were  removed  later  in  the  work 
of  Christensen  and  Lo  [23]  and  the  result  for  G23  was  given  in  the  form  of  a  root  of  a 
quadratic.  Once  the  value  of  G^is  calculated  from  the  Three-phase  model  re¬ 
sults,  E2  and  V23  can  be  determined  by  using  Eqs.  (2.13). 

Figs.  2.3-2.5  show  the  comparison  of  the  effective  elastic  moduli  as  predicted  by 
CCA  model  and  by  the  micro-model.  Results  are  presented  primarily  for  the  transverse 
moduli,  since  these  are  the  ones  which  are  the  most  difficult  to  predict,  and  they 
therefore  serve  as  a  measure  of  accuracy  of  the  method.  In  instances  where  the  CCA 
model  does  not  give  definite  expressions  for  the  effective  moduli  of  the  composite, 
Christensen  and  Lo’s  results  are  employed  along  with  Eqs.  (2.13).  The  assumed 
properties  of  the  fiber  and  matrix  are  the  same  as  used  by  Hashin  and  Rosen  [25]  and 
are  given  in  Table  2.1  for  convenience. 
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Fig.  2.3  shows  the  variation  of  the  transverse  Young’s  modulus  with  fiber  volume 
fraction.  Results  based  on  Christensen  and  Lo’s  solution  and  the  micro-model  are 
practically  coincident.  In  fact  the  micro-model  results  have  been  slightly  offset  to 
improve  the  readability  of  the  graph.  Both  results  are  within  the  bounds  given  by  the 
CCA  model.  Also  shown  in  Fig.  2.3  is  the  transverse  modulus  obtained  by  applying  a 
combination  type  model  to  the  unit  cell  shown  in  Fig.  2.1,  i.e.  by  using  the  second  analog 
model  shown  in  Fig.  2.2,  while  relaxing  the  matrix  restraint  effect  induced  by  the 
compatibility  of  fiber  direction  strain  in  the  fiber  and  matrix.  The  difference  between 
the  combination  model  and  the  Reuss  lower  bound  arises  from  the  fact  that,  in  obtaining 
the  Reuss  lower  bound,  all  of  the  matrix  material  is  lumped  in  series  with  the  fiber 
phase,  whereas  the  combination  model  assumes  that  a  portion  of  the  matrix  (Material 
Element  B)  is  acting  in  parallel  with  the  rest  of  the  material  of  the  unit  cell.  The  latter 
arrangement  results  in  increased  stiffness.  The  difference  between  the  micro-model 
values  and  the  combination  model  values  arises  from  the  matrix  restraining  effect 
induced  by  the  fibers. 

Fig.  2.4  shows  the  longitudinal  shear  modulus  G\  and  the  transverse  direction 
shear  modulus  Q&  as  given  by  the  various  analytical  models  and  the  micro-model.  An 
interesting  fact  is  observed  in  this  result:  the  values  of  the  longitudinal  shear  modulus 
G\  and  the  transverse  shear  modulus  673  are  practically  the  same  according  to  both 
analytical  models  for  the  case  when  the  constituent  phases  are  isotropic.  In  fact,  the 
micro-model  predictions  for  the  two  shear  moduli  are  exactly  the  same  in  this  case.  The 
values  of  shear  moduli  given  by  the  analytical  methods  and  the  micro-model  are  in  good 
agreement  over  the  entire  range  of  the  fiber  volume  fraction. 

Fig.  2.5  shows  the  micro-model  results  for  V23 ,  the  bounds  given  by  the  CCA 
model  and  the  values  given  by  Christensen’s  results  in  conjunction  with  Eqs.  (2.13).  The 
micro-model  results  and  Christensen  and  Lo’s  [23]  results  lie  within  the  CCA  bounds. 
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with  Christensen  and  Lo’s  results  tending  to  favour  the  lower  bound  and  the  micro-mo¬ 
del  results  lying  closer  to  the  upper  bound. 

2.8.2  Comparison  with  Detailed  Numerical  Procedures 
and  Aboudi’s  Method  of  Cells 

Chen  and  Cheng  [33]  determined  the  transverse  moduli  of  a  fiber-reinforced 
composite  with  fibers  arranged  in  a  hexagonal  packing  arrangement,  using  a  series-type 
elasticity  solution.  Fig.  2.6  shows  the  transverse  Young’s  modulus  and  axial  shear 
modulus  for  a  composite  lamina  as  a  function  of  fiber  volume  fraction.  The  properties 
of  fiber  and  matrix  constituents  are  given  in  Table  2.2.  The  micro-model  predictions  are 
in  excellent  agreement  with  the  results  of  Chen  and  Cheng  [33]. 

Fig.  2.7  shows  estimates  of  transverse  and  axial  shear  moduli  for  a  graphite/epoxy 
lamina  from  the  micro-model,  the  numerical  procedure  of  Chen  and  Cheng  [33],  and 
Aboudi’s  Method  of  Cells  [7].  The  assumed  properties  of  the  graphite  fibers  and  epoxy 
matrix  are  given  in  Tkble  2.3.  The  results  from  the  three  procedures  are  in  excellent 
agreement 

Pickett  [34]  presented  numerically  obtained  estimates  for  elastic  moduli  of  a  glass/ 
epoxy  lamina.  The  assumed  properties  of  the  glass  fibers  and  epoxy  matrix  are  presented 
in  Thble  2.4.  Fig.  2.8  presents  the  variation  of  elastic  moduli  of  gl ass/ epoxy  with  fiber 
volume  fraction.  Results  are  presented  for  the  micro-model,  Pickett’s  procedure  [34],  and 
Aboudi’s  micromechanical  model  [7].  The  micro-model  results  compare  very  well  with 
the  results  by  Pickett  and  predictions  from  Aboudi’s  model. 

Figs.  2.9-2.11  show  variation  of  elastic  moduli  with  fiber  volume  fraction  for  a 
graphite/epoxy  composite.  Experimental  results  were  presented  by  Kriz  and  Stinchcomb 
[74].  Results  from  the  micro-model  and  Aboudi’s  Method  of  Cells  [7]  are  also  presented, 
and  are  in  excellent  agreement.  The  agreement  between  analytical  predictions  of  elastic 
moduli  and  experimental  results  is  very  good  in  all  the  cases  except  for  the  transverse 
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Poisson’s  ratio,  in  which  case  the  analytical  results  seem  to  slightly  over-predict  the  Pois¬ 
son’s  ratio. 

2.9  Results  for  Thermal  Expansion  Coefficients 

In  this  section,  the  predictions  of  effective  coefficients  of  thermal  expansion  from 
the  micro-model  are  compared  with  results  from  Aboudi  [54].  Figs.  2.12  and  2.13  show 
the  coefficients  of  thermal  expansion  in  the  axial  and  transverse  directions  for  a  glass/e¬ 
poxy  composite.  The  properties  of  the  the  glass  fibers  and  epoxy  matrix  are  the  same  as 
used  by  Aboudi  [54],  and  are  presented  in  Ihble  2.6  for  convenience.  Figs.  2.14  and  2.15 
show  the  coefficients  of  thermal  expansion  of  a  graphite/epoxy  composite.  The  properties 
of  the  fiber  and  the  matrix  are  given  in  Ihble  2.7. 

Examination  of  the  results  shows  that  the  presence  of  even  a  small  amount  of  fiber 
in  the  matrix  results  in  a  sharp  decrease  in  the  axial  coefficient  of  thermal  expansion, 
whereas  the  effect  of  fibers  in  reducing  the  coefficient  of  thermal  expansion  in  the  trans¬ 
verse  direction  is  limited.  In  fact  it  is  seen  that  for  low  fiber  volume  fractions  there  can 
even  be  an  increase  in  the  transverse  coefficient  of  thermal  expansion.  This  seemingly 
anomalous  behaviour  is  a  result  of  the  constraining  effect  of  the  fibers.  When  fibers, 
which  are  axially  stiff  and  have  small  thermal  expansion  coefficients,  are  added  to  a  rela¬ 
tively  more  compliant  matrix  with  a  large  expansion  coefficient,  the  effect  is  a  sharp  drop 
in  the  axial  thermal  expansion  coefficient  even  for  low  fiber  volume  fractions.  This  is  due 
to  high  fiber  stiffness  in  axial  direction,  but  the  tendency  of  the  matrix  to  expand  in  the  the 
axial  direction  is  channeled  into  the  transverse  direction,  causing  an  increase  in  the  effec¬ 
tive  transverse-direction  expansion  coefficient  As  the  fiber  volume  fraction  increases, 
the  transverse  expansion  coefficient  gradually  decreases,  more  or  less  linearly. 

The  predictions  of  thermal  expansion  coefficients  from  the  proposed  micro-model 
and  from  Aboudi’s  refined  micromechanical  model  are  in  excellent  agreement  and  the 
micro-model  properly  portrays  the  anomalous  behaviour  of  unidirectional  composites  in 
the  case  of  transverse  thermal  expansion. 
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I  m  n 

i  i  i 


Fig.  2.1  Idealized  composite  lamina  and  representative  unit  cell 


For  component  1 1  For  components 

22,33,12,23,13 

Fig.  2.2  Analogs  for  stress  and  strain  homogenization  in  the  unit  cell 
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E 

V 

(Mpsi) 

E-Glass 

10.5 

0.2 

Epoxy 

0.4 

0.35 

Table  1.1  Elastic  properties  of  E-glass  fibers  and  epoxy 
matrix  used  in  comparing  Micro-model,  CCA 
model  [26]  and  the  Three-phase  model  results. 


Fiber  Volume  Fraction  Vf 

Fig.  2.3  Comparison  of  transverse  modulus  calculated  from  the  Micro¬ 
model,  CCA  model  [26],  and  the  Three-phase  model[23]  In 
conjunction  with  Eq  (2.13). 
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Transverse  Poisson’s  ratio*  v23  Shear  Moduli,  G1  and  G23  (Mpsi) 


Fiber  Volume  Fraction  Vf 

Fig.  2.4  Shear  moduli  calculated  from  Micro-model,  CCA 
model  [26]  and  the  Three-phase  model[23]. 


Fiber  Volume  Fraction  Vf 


Fig.  2.5  Transverse  Poison’s  ratio  calculated  from  Micro-model, 
CCA  model  [26],  and  the  Three-phase  model [23]  in 
conjunction  with  Bq.  (2.13). 
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Elastic 


E1 

e2 

Gj2 

Vl2 

V23 

(Mpsi) 

(Mpsi) 

(Mpsi) 

Fiber 

2 

10 

0.2 

0.15 

Matrix 

1  : 

0.6 

0.222 

0.35 

0.35 

Table  2.2  Elastic  properties  of  fibers  and  matrix  used  In  comparison  of  Micro-model 
predictions  and  results  of  Chen  and  Cheng [33 J. 


*  Fiber  Volume  Fraction  Vf 


Fig.  2.6  Comparison  of  Micro-Model  results  with  numerically 
obtained  results  of  Chen  and  Cheng[33J. 


SO 


Elastic  Modulus  (Mpsi) 


Table  2.3  Elastic  properties  of  graphite  fibers  and  epoxy  matrix  used  in  comparison 
of  Micro-model  predictions,  results  of  Chen  and  Cheng  [33]  and  Aboudi's 
Method  of  Cells  [7] 


Fiber  Volume  Fraction  Vf 


Fig.  2.7  Comparison  of  Micro-model  results  with  numerically 
obtained  results  of  Chen  and  Cheng [33]  and  Aboudi's 
Method  of  Cells  17] . 
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Elastic 


Table  2.4  Elastic  properties  of  glass  fibers  and  epoxy 
matrix  used  in  comparing  Micro-model 
predictions  and  results  of  Pickett  [34]  and 
Aboudi  [7]. 


I 

I 


0.0  0.2  0.4  0.6  0.8  1.0 


Fiber  Volume  Fraction  Vf 

Fig.  2.8  Comparison  of  Micro-model  results  with  numerically  obtained 
results  of  Pickett [34]  and  Aboudi 's  Method  of  Cells  [7], 
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Elastic  Modulus  (Mpsi) 


Ex 

e2 

Gi2 

i'll 

V23 

(Gpa) 

(Gpa) 

(Gpa) 

Graphite 

232 

15 

24 

0.279 

0.49 

Epoxy 

5.35 

5.35 

1.975 

0.354 

0.354 

Table  2.5  Elastic  properties  of  graphite  fibers  and  epoxy  matrix  used  In 

comparison  of  Micro-model  Predictions  and  results  from  Aboudl’s 
Method  of  Cells  [7J. 


Fiber  Volume  Fraction  Vf 


Fig.  2.9  Comparison  of  transverse  Young's  modulus  and  axial  shear  modulus 
predictions  from  Micro-model  and  Aboudl’s  Method  of  Cells  [7]. 


Transverse  Poisson's  ratio  •  v23  Transverse  shear  Modulus  G23  (Qpa) 


Fiber  Volume  Fraction  Vf 

Fig.  2.10  Comparison  of  transverse  shear  modulus  predic¬ 
tions  from  Micro-model  and  Aboudi’s  Method  of 
Cells  [7]. 


Fiber  Volume  Fraction  Vf 


Fig.  2.11  Comparison  of  transverse  Poisson's  ratio  predic¬ 
tions  from  Micro-model  and  Aboudi's  Method  of 
Cells  [  7] . 
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E 

V 

a 

(GPa) 

1  o-V0c 

Glass 

72.38 

0.2 

5 

Epoxy 

2.75 

0.35 

54 

Table  2.6  Properties  of  glass  fibers  and  epoxy  matrix  used  in 
comparing  predictions  of  coefficients  of  thermal 
expansion  from  Micro-model  and  Aboudl’s  Method 
of  Cells  [54]. 


Fiber  Volume  Fraction  Vf 

no-  2.12  Prediction  of  axial  thermal  expansion  coefficient 
of  glass/epoxy  from  Micro-model  and  Aboudi's 
Method  of  Cells  [54]. 
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Fiber  Volume  Fraction  Vf 

Fig.  2.13  Prediction  of  transverse  thermal  expansion  coeffi¬ 
cient  of  Glass! Epoxy  from  Micro-Model  and 
Aboudi's  Method  of  Cells[54]. 


Ei 

Ez 

Git 

Vl2 

*23 

a2 

(GPa) 

(GPa) 

(GPa) 

KT4/  °c 

10  ~*/°C 

Graphite 

388 

7.8 

6.8 

0.23 

0.45 

-0.68 

9.74 

Epoxy 

3.45 

3.45 

1.27 

0.35 

0.35 

36 

36 

Table  2.7  Properties  of  graphite  fibers  and  epoxy  matrix  used  in  comparison  of 

coefficients  of  thermal  expansion  predictions  from  Micro-model  and  Aboudi's 
Method  of  Cells  [54]. 


(V8-0t) 


Fiber  Volume  Fraction  Vf 


Fig.  2.14  Prediction  of  axial  thermal  expansion  coefficient 
of  graphlte/epoxy  lamina  from  Micro-model  and 
Aboudl's  Method  of  Cells  154]. 


Fiber  Volume  Fraction  Vf 

Fig.  2.15  Prediction  of  transverse  thermal  expansion  coeffi¬ 
cient  of  graphitelepoxy  from  Micro-model  and 
Aboudi's  Method  of  Cells  [54]. 
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CHAPTER  3 


SUBLAMINATE  MODEL 


A  group  or  a  duster  of  laminae,  from  among  the  total  number  of  laminae  consti¬ 
tuting  the  laminate,  is  called  a  sublaminate.  The  micro-model  combines  the  stresses  and 
strains  within  the  fiber  and  matrix  phases  to  yield  the  stresses,  strains  and  stiffness  of  an 
effective  homogeneous  lamina.  This  information  is  then  passed  onto  the  sublaminate 
model,  as  shown  in  Fig.  3.1,  which  schematically  describes  the  working  of  the  sublaminate 
model.  The  sublaminate  model  then  combines  the  stress  and  strain  fields  within  each 
lamina  to  yield  the  stresses,  strains  and  stiffness  of  an  equivalent  homogeneous  anisotropic 
material ,  which  represents  the  sublaminate. 

Quite  often,  within  the  stacking  sequence  of  a  laminate,  a  typical  repeating  pattern 
is  identifiable.  This  is  usually  the  case  for  thick-section  composites,  where  the  symmetry 
resulting  from  using  the  same  pattern  throughout  the  laminate  thickness  is  desirable  in 
order  to  avoid  the  high  curing  stresses  which  would  otherwise  arise.  An  example  of  a 
thick-section  laminate  with  a  repeating  pattern  is  a  [02/90]i&  laminate,  which  consists  of 
O2/9O  pattern,  repeated  32  times  through  the  thickness  of  the  laminate  for  a  total  of  96 
individual  laminae.  Here,  it  can  be  reasonably  assumed  that  the  material  in  the  vicinity  of 
a  typical  material  point  consists  of  3  plies  arranged  in  this  basic  pattern.  If  these  layers  are 
replaced  by  an  equivalent  homogeneous  material  then  we  can  employ  standard  finite  ele¬ 
ments  and  the  usual  sampling  procedures  at  the  material  points  to  analyze  this  particular 
laminate.  The  prospective  advantages  of  this  method  also  become  obvious  now,  since  we 
can  carry  out  the  analysis  by  sampling  at  discrete  locations  through  the  laminate  thickness, 
as  compared  to  ply-by-ply  type  analysis  procedures  which  require  considerably  more 
computational  effort  and  storage,  in  order  to  track  the  nonlinear  response  of  a  laminate. 
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The  method  of  representing  the  cluster  of  plies  forming  the  sublaminate  by  an 
equivalent  homogeneous  continuum  maybe  referred  to  as  “Equivalent  Continuum  Mod* 
elling.”  In  the  next  section  the  general  ideas  behind  the  method  are  presented. 

3.1  Equivalent  Continuum  Modelling 

The  idea  of  equivalent  continuum  modelling  has  been  applied  in  many  different 
areas  of  mechanics  and  to  materials  other  than  laminated  composites.  Here,  discussion  is 
focused  upon  the  method  as  it  has  been  applied  to  laminated  composites. 

The  basic  procedure  for  obtaining  the  properties  of  an  equivalent  continuum  rep¬ 
resenting  a  laminated  material  is  as  follows.  The  actual  laminated  material  is  replaced  by 
an  equivalent  homogeneous  material.  The  properties  of  this  equivalent  material  are  de¬ 
termined  by  requiring  that  the  actual  material  and  the  equivalent  material  behave  in  the 
same  way  when  subjected  to  certain  fundamental  patterns  of  stresses  or  strains.  These 
fundamental  patterns  of  stre&jes  or  strains  may  be  thought  of  as  calibrating  stress  or  strain 
fields.  The  properties  of  the  equivalent  material  so  obtained  can  be  used  to  analyze  the 
structure  under  consideration  and  should  yield  reasonably  accurate  results,  provided  that 
the  stress  and  the  strain  fields  at  the  material  points  of  the  structure  are  not  too  different 
from  the  calibrating  stress  and  strain  fields. 

3.2  Previous  Work  and  Adopted  Procedure 

Pagano  [75]  gave  expressions  for  the  effective  elastic  moduli  of  laminated  compos¬ 
ites,  using  a  three-dimensional  lamination  scheme.  Pagano  considered  all  the  laminae 
within  the  thickness  of  the  laminate.  It  was  suggested  that  moduli  so  determined  could  be 
used  for  the  purpose  of  structural  analysis,  as  opposed  to  laminated  plate  type  analysis 
procedures,  which  operate  on  a  ply-by-ply  basis.  Once  the  displacement  field  is  known, 
the  detailed  ply-by-ply  solution  can  be  determined  from  it.  A  similar  procedure  was 
employed  by  Sun  and  Li  [76]  to  determine  the  effective  elastic  moduli  of  a  typical  repeat¬ 
ing  sublaminate  within  the  thickness  of  the  laminate.  It  was  suggested  that  since  the  lami- 
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nate  was  made  up  from  this  repeating  sublaminate,  the  elastic  moduli  of  the  laminate 
would  be  the  same  as  the  elastic  moduli  of  the  sublaminate. 

A  sublaminate  model  for  determining  equivalent  elastic  properties  of  a  sublami¬ 
nate  from  the  properties  of  constituent  laminae  was  presented  by  Pecknold  [57].  The 
choice  of  the  calibrating  stress  and  strain  fields  in  Pecknold’s  work  is  essentially  the  same 
as  the  one  employed  by  Pagano  [75]  and  Sun  and  Li  [76];  therefore  it  was  suggested  by 
Pecknold  that  the  results  obtained  by  sublaminate  model  [57]  should  be  identical  to  re¬ 
sults  from  procedures  proposed  by  Pagano  and  Sun  and  Li.  However,  the  formulation 
presented  in  Pecknold’s  work  is  quite  simple,  and  this  author  is  in  agreement  with  the  as¬ 
sertion  made  in  Pecknold’s  work  [57]  that  the  formulation  presented  in  the  works  of  Paga¬ 
no  [75]  and  Sun  and  Li  [76]  is  such  that  it  obscures  the  essential  simplicity  of  the  basic 
ideas  which  are  involved.  Therefore,  in  the  present  study  the  sublaminate  model  as  pro¬ 
posed  by  Pecknold  [57]  has  been  employed.  In  this  chapter,  certain  characteristics  of  the 
stress  and  strain  distributions  within  laminated  composites  are  first  discussed,  since  these 
have  a  bearing  on  the  type  of  calibrating  stress  and  strain  fields  which  are  permissible. 
The  procedure  for  obtaining  the  equivalent  continuum  properties  of  the  sublaminate  is 
then  formally  outlined. 

3.3  Characteristics  of  Stress  and 
Strain  Fields  in  Laminates 

Fig.  3.2  shows  a  sublaminate  in  assembled  and  disassembled  forms.  The  coordi¬ 
nate  axes  attached  to  sublaminate  are  referred  to  as  laminate  coordinate  axes',  the  primed 
coordinate  axes  attached  to  the  lamina  are  the  lamina  coordinate  axes.  The  kth  lamina  has 
also  shown  on  its  surface  various  stresses  referred  to  the  laminate  coordinate  system.  The 
11, 22,  and  12  components  of  stress  and  strain  are  referred  to  as  in-plane  stresses  and 
strains.  Correspondingly  the  33,23,  and  13  components  of  stress  and  strain  are  referred  to 
as  out-of-plane  stresses  and  strains.  The  vectors  of  stress  and  strain  can  be  partitioned 
using  this  classification  of  stresses  and  strains  as  follows: 
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Referring  to  Fig.  3.2,  if  there  exists  a  perfect  bond  between  two  laminae  at  their 


interface,  then  the  compatibility  of  displacements  at  the  interface  requires,  not  only  that 


the  displacements  u\>  us  be  continuous,  but  also  that  the  in-plane  derivatives  of 
these  displacements  be  continuous.  This  in  turn  translates  into  the  requirement  that  the 
in-plane  strains  should  be  continuous  at  the  lamina  interface,  i.e.. 


(3.2) 


Interface  Interface 

The  equilibrium  of  tractions  at  the  interface  requires  that  the  out-of-plane  stresses  be 


continuous  across  the  lamina  interface,  i.e.. 


r  1* 

r  i*+J 

H 

=  H 

Interface  Interface 


These  two  requirements  can  be  written  compactly  as 


(3.3) 


(3.4) 


Interface  Interface 

The  nature  of  the  remaining  components  of  stress  and  strain,  i.e.  the  in-plane 


stresses  and  out-of-plane  strains,  is  dictated  by  these  requirements,  and  is  explained  as 


follows. 


Referring  again  to  Fig.  3.2,  it  is  obvious  that  if  the  elastic  moduli  of  various  laminae 
are  referred  to  the  global  coordinates,  then  there  will  often  exist  a  discontinuity  in  elastic 
moduli  at  the  lamina  interface,  since  the  fiber  orientations  in  the  two  layers  forming  the 
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interface  may  be  different.  If  this  is  the  case,  and  the  requirement  of  the  continuity  of 
in-plane  strains  is  imposed,  then  obviously  there  would  be  a  discontinuity  in  the  in-plane 
stresses  at  the  lamina  interface.  If  the  requirement  of  the  continuity  of  the  out-of-plane 
stresses  is  imposed  at  the  interface,  then  the  out-of-plane  strains  may  be  discontinuous. 
The  magnitude  of  the  discontinuities  in  the  in-plane  stresses  and  out-of-plane  strains  de¬ 
pends  upon  the  degree  of  mismatch  between  the  elastic  moduli  of  the  two  layers  forming 
the  interface  and  upon  the  nature  of  the  imposed  loading. 

The  problems  encountered  by  conventional  displacement-based  finite  elements, 
such  as  the  elements  based  on  First-order  Shear  Deformation  Theory  and  Higher-order 
Shear  Deformation  Theories  for  laminates,  can  now  be  better  realized.  These  elements 
employ  smooth  strain  distributions  within  a  single  element.  This  results  in  assumption  of 
smoothly  varying  out-of-plane  strains  through  the  laminate  thickness.  The  computed 
out-of-plane  stresses  may  then  exhibit  discontinuities  at  the  laminae  interfaces,  violating 
the  requirement  of  continuity  of  out-of-plane  stresses.  Accurate  assessment  of  out-of- 
plane  stresses  in  some  cases  may  be  important,  since  these  stress  components  play  an  im¬ 
portant  role  in  causing  delaminations. 

In  the  following  section  the  transformation  of  stresses,  strains,  and  material  prop¬ 
erties  from  lamina  coordinates  to  laminate  coordinates  is  described. 

3.4  Transformation  of  Stresses  and  Strains 

The  stresses,  strains,  and  stiffness  matrix  of  a  lamina  in  lamina  coordinates  are 
denoted  by  o',  e\  andC',  respectively.  The  same  quantities  when  referred  to  the  lami¬ 
nate  coordinate  system  are  denoted  by  a,  6,  andC,  respectively.  The  axes  A3  and  A3' 
are  colinear,  the  fiber  orientation  angle  is  denoted  by  0 ,  and  the  superscript  k  denotes  the 
kth  lamina.  The  stress  and  strain  transformation  laws  are 

M  =  M  W  ■  M  =  [t/]  W  (3-5) 


* 
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where  the  strains  and  stresses,  in  the  vector  of  strains  {e }  and  stresses  [o } ,  are  arranged 
according  to  the  convention  employed  in  Eq.  (2.5) 


The  transformation  laws  describing  the  transformation  of  the  lamina  stiffness  ma¬ 
trix  and  compliance  matrix  from  lamina  coordinates  to  laminate  coordinates  follow  from 
Eq.  (3.5)  and  the  invariance  of  strain  energy  under  coordinate  transformation. 
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The  strain  and  stress  transformation  matrices,  [  T*  ]  and  [  T*  ]  are  defined  as: 
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(3.7b) 


where,  /i,  h,  mi,  m2  are  direction  cosines,  given  by 

/,  ”  cos(  XuXS  )  -  cos(0*)  ,  m,  -  cosiX^Xi'  )  -  sin(0*) 

h  -  cos(XuX2'  )  -  -sin(0*)  ,  m2  -  cos(  X^XJ  )  -  cos(0*) 


(3.8) 
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Some  additional  useful  relations  between  the  strain  and  stress  transformation  ma¬ 
trices  [  T*  ]  and[  T*  ]  are 

M  ■  M  M  =  M  (3  -9) 

3.5  Calibrating  Stress-Strain  Fields 

In  order  to  determine  the  equivalent  stiffness  of  the  sublaminate,  it  must  be  sub¬ 
jected  to  some  specified  calibrating  patterns  of  stresses  and  strains.  Also,  preferably,  the 
equivalent  stiffness  of  the  sublaminate  should  be  determined  in  such  a  manner  that  the 
stress  and  strain  continuity  requirements  specified  in  Eq.  (3.4)  are  satisfied.  If  it  is  as¬ 
sumed  that  stresses  and  strains  in  each  lamina  within  the  sublaminate  are  constant,  then 
the  following  calibrating  pattern  of  stresses  and  strains  can  be  proposed,  which  also  satis¬ 
fies  the  requirements  in  Eq.  (3.4): 

ioxk  —  l, ..  .N  (3.10) 

where  ej-  and  aQ  are  homogenized  in-plane  strains  and  out-of-plane  stresses  in  the 
equivalent  continuum  representing  the  sublaminate,  and  JV  is  the  number  of  laminae  com¬ 
prising  the  sublaminate.  Eq.  (3.10)  suggests  that  the  calibrating  stress  and  strain  fields 
chosen  are  uniform  in-plane  strains  and  uniform  out-of-plane  stresses. 

In  the  next  section,  homogenization  relations  will  be  presented  for  the  remaining 
components  of  stress  and  strain  in  the  equivalent  material. 

3.6  Homogenized  In-piane  Stresses  and  Out-of-plane  Strains 

The  in-plane  stresses  and  out-of-plane  strains  in  the  equivalent  continuum  are 
defined  as 

R}  ■  l 
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where  f*  is  the  thickness  of  the  kth  lamina,  t  is  the  thickness  of  the  sublaminate,  and  N  is 
the  number  of  laminae  in  the  sublaminate. 

Eqs.  (3.10)  and  (3.11)  are  sufficient  to  completely  characterize  the  sublaminate 
model.  The  above  relations  also  apply  to  incremental  stresses  and  strains.  In  the  next 
section  these  relations  are  employed  to  obtain  the  tangent  stiffness  of  the  equivalent  con¬ 
tinuum  representing  the  sublaminate. 


3.7  Tangent  Stiffness  of  Sublaminate 

The  tangent  stiffness  relations  for  a  lamina  are  given  by  Eq.  (2.12).  In  the  lamina 
coordinate  system,  they  take  the  form 

\d(f\  =  [  C'*]  {dej  +  AT  {jff'}  (3-12) 

Eq.  (3.12)  is  transformed  to  the  laminate  coordinate  system  using  Eqs.  (3.5-3 .8)  to  yield 

[dof  =  [c*]  {*•}*+ AT  {fif  (3-13) 


The  vector  of  thermal  properties  {/?'}  transforms  like  the  vector  of  stresses,  i.e., 

i  k  P  1  \  k 


W  ■  [-.*]  w 


(3.14) 


The  lamina  stiffness  relations  in  laminate  coordinates  can  be  written  in  partitioned 
form,  in  accordance  with  Eq.  (3.1)  as 


-c'-l  1*4  +  at  tU 

I dO0)  L  !  C00  J  [d€ol  [po  ] 


(3.15) 


The  above  lamina  tangent  stiffness  relation  is  now  partially  inverted  to  give 
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(3.16) 


w  =  r 
Ul  L-* 


£-]{£}  +Arfe_] 


where 


P OO  -  C  OO 


P io  —  C io  Poo 
Pit  =  Cjj  -  Pio  C oi 
Pfio  ~  ~  p°°  A> 

P^  -  A  +  P&  fio 


(3.17) 


The  partially  inverted  stiffness  relations  for  the  individual  laminae  within  the  sublami¬ 
nate  are  assembled  according  to  Eqs.  (3.10)  and  (3.11)  to  form  partially  inverted  stiffness 
relations  for  the  sublaminate  : 


where 


L-f'«  1  F~J  W  UJ 
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(3.18) 


(3.19) 


The  partially  inverted  stiffness  relations  (Eq.  (3.18))  can  now  be  partially  re-in¬ 
verted  to  obtain  the  tangent  stiffness  relations  for  the  equivalent  material  representing  the 
sublaminate 


ui  _  r_* 
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(3.20) 
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where 


A>  =  ~^oo 

~fii  ~  P/fy  +  P/o  fio 


(3-21) 


If  the  elastic  moduli  of  the  sublaminate  are  of  interest,  they  can  be  obtained  from  the  com¬ 
pliance  matrix  of  the  sublaminate,  which  is  found  by  inversion  of  the  stiffness  matrix  in 
Eq.  (320). 


Micro-Model 


Fig.  3.1  Schematic  representation  of  the  functioning  of 
the  Sublaminate  model 
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Disassembled  Sublaminate 

Fig.  3.2  Stresses  on  laminae  and  the  various  coordinate  systems. 
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CHAPTER  4 


THE  STRESS  UPDATE 


X  he  solution  of  finite  element  models  of  nonlinear  structural  systems  is  usually 
obtained  by  applying  the  loading  in  an  incremental  manner.  In  addition,  after  application 
of  each  load  increment,  stresses  are  corrected  and  equilibrium  iterations  are  normally 
performed  to  ensure  equilibrium  to  within  a  specified  tolerance.  As  the  loading  is  applied 
to  the  structure,  its  geometry  and  the  properties  of  the  material  may  change,  making  the 
response  of  the  structure  nonlinear.  These  changes  in  material  properties  are  monitored 
at  the  material/sampling  points  within  each  element  The  effect  of  material  property 
changes  is  reflected  in  the  element  stiffness  matrices,  which  in  turn  convey  the  effect  of 
these  changes  to  the  structure  equilibrium  equations,  when  element  stiffness  matrices  are 
assembled  to  form  the  structure  stiffness  matrix.  The  segment  of  the  structural  analysis 
program  which  provides  the  constitutive  description  of  the  material  at  a  sampling  point 
and  serves  to  update  the  material  properties  and  the  stresses  at  the  sampling  point  is  called 
the  material  model. 

Specifically  stated,  within  the  framework  of  finite  element  analysis,  the  material 
model  is  called  upon  to  perform  the  following  important  tasks: 

1)  Update  the  stresses  at  the  material  points,  as  the  strains  change. 

2)  Update  the  tangent  stiffness  of  the  material  at  the  material  points,  if  re¬ 
quired. 

3)  Record  and  indicate  other  changes  in  the  state  of  the  material,  such  as  the 
onset  of  yielding  or  failure  etc. 

4.1  Solution  of  Nonlinear  Structural  Systems 

The  solution  of  nonlinear  systems  is  usually  carried  out  by  a  combined  incremen¬ 
tal-iterative  scheme,  nidi  as  the  Newton-Raphson  or  the  Modified  Newton-Raphson 
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method.  The  loads  are  applied  in  increments  P„,  pn+1, ... ,  and  at  each  load  step  equilib¬ 
rium  iterations  are  performed  to  satisfy  equilibrium  to  within  specified  tolerance.  The 
equilibrium  equations  of  the  structure  are  cast  in  an  incremental-iterative  form  as 

kt  Pn+i  -  inf;V  =  (4.i) 

in  which  KT  is  the  tangent  stiffness  matrix  of  the  structure,  IRn  +  .  is  the  internal  resis¬ 
tance  vector  and  R„  +  x  is  the  residual  load.  The  superscripts  i  and  i+1  refer  to  ith  and 

i+lth  iterations.  After  each  iteration  the  previous  estimate  of  nodal  deflections  is  up¬ 
dated  as 


«?+i  -  i  <4-2> 

The  graphical  interpretation  of  the  incremental-iterative  scheme  is  shown  in  Fig.  4.1. 

4.1.1  Role  of  the  Material  Model 

The  tangent  stiffness  matrix  of  the  structure  Kr  and  the  internal  resistance  vector 
I*  are  formed  by  assembling  element  tangent  stiffness  matrices  and  internal  resis¬ 
tances  vectors  I* .  The  element  stiffness  matrices  and  internal  resistance  vectors  are 
evaluated  using  a  numerical  integration  scheme  such  as  Gaussian  quadrature.  The  nu¬ 
merical  integration  scheme  requires  the  material  tangent  stiffnesses  and  stresses  at  the 
material  points,  as  is  indicated  in  the  following  Eqs.  (4.3): 


K?=  B(e)TC?  B WrfFW 

v*> 

=  [  BwVe) 

jX«) 

where 


NP 
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NP  =  Number  of  Material  Points  or  Sampling 
Points  within  the  element 


(43) 


71 


Weighting  Factor  associated  with  the  ith 
Material  Point 


Wi 


J,  =  Jacobian  Matrix  defining  the  mapping  from  element 
coordinates  to  parametric  coordinates. 

The  function  of  the  material  model  becomes  clear  when  Eqs.  (4.3)  are  considered.  Corre¬ 
sponding  to  a  deflected  state,  strain  increments,  previous  stresses  and  strains,  and  history 
parameters  are  passed  to  the  material  model,  which  then  updates  the  stresses,  material 
tangent  properties,  and  history  parameters  at  the  material  point.  Geometric  nonlineari¬ 
ties,  if  present,  are  reflected  in  the  strain-displacement  matrix  B . 

The  procedure  for  obtaining  the  equivalent  tangent  stiffness  of  the  material  at  the 
material  point  has  been  described  in  detail  in  Chapters  1-3,  and  is  quite  straight-forward. 
The  stress  update,  however,  is  more  involved  as  a  consequence  of  the  “two-level  micro¬ 
structure”  employed  in  the  material  modelling  procedure.  The  problems  associated  with 
updating  the  stresses  at  a  material  point  will  be  outlined  in  detail  in  the  next  section. 

4.2  Outline  of  the  Stress  Update  Problem 

Suppose  that  during  the  incremental-iterative  solution  scheme,  a  change  AU 
takes  place  in  the  nodal  displacements.  The  displacement  increment  causes  an  increment 
of  strain  Ae  at  a  material  point  The  state  of  strains  at  the  material  points  within  the 
elements  is  directly  updated,  i.e., 

e  *-  e  +  Ae  (4.4) 

The  stresses,  however,  cannot  be  updated  directly  by  adding  the  predicted  increments  of 
stresses  A o  based  on  tangential  estimates  to  the  previous  state  of  stress,  since  this  would 
result  in  a  cumulative  error  as  the  solution  proceeds.  This  problem  is  a  direct  consequence 
of  nonlinear  material  behaviour,  as  is  illustrated  in  Fig.  4.2,  which  shows  the  difference 
between  the  stress  increment  based  on  tangential  estimate  Ao  and  the  true  stress  incre- 
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ment  A a .  What  is  sought  is  an  accurate  estimate  of  the  true  stress  increment  Act  consistent 
with  the  material  constitutive  description  and  the  imposed  strain  increment  A?. 

The  procedure  for  determining  the  true  increment  of  stress  Act  corresponding  to 
an  increment  of  strain  Ae  is  more  complicated  than  for  metal-plasticity.  This  is  a  result  of 
the  two-level  microstructure  of  the  material  model  employed  here.  The  problem  posed  in 
updating  the  stresses  at  a  material  point  is  described  in  Fig.  4.3,  which  shows  the  equiva¬ 
lent  homogeneous  material  at  a  typical  material  point  undergoing  an  increment  of  strain. 
The  problem  is  to  find  the  corresponding  true  increment  of  stress  in  the  equivalent  materi¬ 
al,  and  update  the  state  of  stress  there.  For  this,  the  effect  of  the  strain  increment  has  to  be 
cascaded  down  to  the  micromechanical  level,  where  the  stresses  in  the  fiber  and  matrix 
phases  are  updated.  Once  the  stresses  in  the  fiber  and  matrix  have  been  updated,  the 
tangent  stiffness  and  stress  in  the  equivalent  material  can  be  evaluated  using  the  stress  and 
strain  homogenization  procedure  in  the  micro-model  and  the  sublaminate  model.  For  the 
homogenization  procedure  to  be  valid,  it  is  essential  that  the  homogenization  relations  are  not 
violated  during  the  process  of  transferring  the  effect  of  the  strain  increment  to  the  fiber  and  the 
matrix  phases.  This  means  that  the  same  homogenization  relations  used  to  define  the  Micro¬ 
model  and  the  Sublaminate  model  now  act  as  constraints  which  must  not  be  violated  during 
the  stress  update  procedure.  Fig.  43  shows  these  constraints/homogenization  relations, 
which  have  to  be  satisfied  at  the  micro-model  level  and  the  sublaminate  model  level.  It  is 
important  to  note  that  the  constraints  are  hybrid  constraints,  i.e.  the  constraints  are  speci¬ 
fied  on  the  strains  as  well  as  the  stresses. 

4.3  Stress  Update  Procedure 

The  procedure  adopted  here  to  solve  the  stress  update  problem  is  a  relaxation  type 
procedure.  The  general  strategy  employed  is  that  some  of  the  constraints  are  imposed 
while  temporarily  relaxing  the  complementary  set  of  constraints.  Residuals  are  then  cal¬ 
culated  by  determining  the  imbalance  in  the  constraints  which  are  violated.  The  out  of 
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balance  constraints  are  satisfied  by  redistributing  the  residuals  to  remove  the  imbalance. 
The  procedure  can  be  repeated  a  number  of  times,  until  the  desired  degree  of  accuracy  is 
achieved. 

The  solution  scheme  is  a  five-step  process.  Description  of  these  steps  and  the 
specific  tasks  performed  during  each  of  these  steps  is  presented  in  the  following  sections. 

4.3.1  Tangential  Update  (Step  1) 

In  this  step  the  effect  of  the  strain  increment  Ae,  prescribed  over  the  equivalent 
material  at  a  material  point,  is  transferred  to  the  fiber  and  matrix  phases  at  the  micro-mo- 
del  level  using  tangential  estimates  based  on  material  properties  at  the  end  of  last  con¬ 
verged  load  step.  The  algorithm  for  the  Ihngent  Update  procedure  is  presented  schemat¬ 
ically  in  Fig.  4.4,  and  the  description  of  the  various  operations  performed  in  this  step 
follows  next. 

Given  the  strain  increment  A?  prescribed  over  the  equivalent  material  at  a  materi¬ 
al  point,  the  in-plane  strains  within  each  lamina  of  the  sublaminate  are  known  via  con¬ 
straint  (SI);  but  the  out-of-plane  strains  within  the  individual  laminae  are  not  known.  At 
the  same  time  the  continuity  of  out-of-plane  stresses  has  to  be  assured.  Keeping  both 
these  considerations  in  mind,  the  tangential  prediction  of  the  increment  of  out-of-plane 
stresses  A (%,  is  accepted  for  the  time  being  as  a  reasonably  good  estimate  of  the  actual 
increment  in  out-of-plane  stresses. 

After  the  in-plane  strains  and  out-of-plane  stresses  have  been  transformed  to  the 
lamina  local  coordinates,  the  complementary  out-of-plane  strains  can  be  determined  by 
forming  the  partially  inverted  tangent  stiffness  of  the  lamina  as  shown  in  Fig.  4.4,  and  in 
the  following  relation 


The  complete  set  of  strains  in  a  lamina  is  now  known.  The  strain  in  the  lamina  is  equiva¬ 
lent  to  the  strain  in  the  unit  cell,  and  by  virtue  of  constraint  (Ml)  is  equal  to  the  increment 
of  strain  in  material  elements  A  and  B.  The  increment  of  stresses  in  material  elements  A 
and  B  can  be  readily  estimated  by  forming  the  respective  tangent  stiffnesses  of  material 
elements  A  and  B, 


(4.6) 


Within  material  element  A  the  state  of  stress  and  strain  within  the  fiber  subcell  f 


and  the  matrix  subcell  iua  can  be  estimated  first  by  invoking  constraint  (M3)  and  then 
determining  the  remaining  stress  and  strain  increments  by  using  the  partially  inverted 
compliance  relations  for  the  fiber  and  the  matrix  subcell  mA . 
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This  completes  the  Thngent  Update  step  of  the  stress  update  procedure. 

4.3.2  Constitutive  Relations  Update  (Step  2) 

After  the  Thngent  Update  step,  the  tangential  estimates  of  stresses  and  strains 
within  the  various  subcells  of  the  unit  cell  are  known.  At  this  point  all  constraints  are 
satisfied.  However,  if  the  material  is  nonlinear,  the  estimated  stresses  and  the  strains 
within  the  subcells  are  not  consistent  with  the  material  constitutive  description.  Some  of 
the  stress  or  strain  components  have  to  be  altered  to  make  the  stresses  and  strains  consis¬ 
tent  with  the  material  constitutive  description.  The  manner  in  which  this  is  done  is  that 
some  of  the  tangential  estimates  of  stress  and  strain  components  are  accepted,  and  the 
remaining  components  of  stress  and  strain  are  determined,  such  that  they  are  consistent 
with  the  constitutive  description  of  the  material. 
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For  matrix  subcell  hib  the  tangential  estimate  of  strain  is  accepted,  and  the  corre¬ 
sponding  set  of  stresses  {a]*,  which  is  consistent  with  the  constitutive  description,  is  deter¬ 
mined.  The  procedure  which  is  used  is  subincrementation,  in  which  the  prescribed  strain 
increment  is  divided  into  a  pre-determined  number  of  smaller  subincrements.  At  the  end 
of  each  subincrement  the  tangent  stiffness  of  the  material  is  updated,  and  a  subincrement 
of  stress  corresponding  to  the  next  subincrement  of  strain  is  calculated.  The  subincre¬ 
ments  of  stress  are  summed  to  give  the  increment  in  stress,  which  is  added  to  the  previous 
state  of  stress  to  give  the  value  of  stresses  at  the  end  of  the  prescribed  strain  increment. 
The  consistent  stresses  within  matrix  subcell  me  are  then  given  by 
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where  M  is  the  number  of  subincrements  within  the  increment.  The  subincrementation 
procedure  is  illustrated  in  Fig.  4.5,  and  is  adopted  here  because  of  its  simplicity.  The 
accuracy  of  the  subincrementation  method  can  be  improved  by  increasing  the  number  of 
subincrements;  alternatively,  more  accurate  procedures  such  as  Runge-Kutta  could  be 
used. 


Within  material  element  A  the  fiber  is  linear-elastic,  and  therefore  the  tangential 
estimates  of  stresses  and  strains  within  the  fiber  are  consistent  with  one  another.  It  should 
be  noted  that  a  nonlinear  constitutive  relation  could  be  used  for  the  fiber,  subincrementa- 
tion  scheme  would  then  be  used  here  also.  Within  matrix  subcell  mu  the  tangential 

estimates  of  longitudinal  strain  tL  and  transverse  stresses  <jT  are  accepted,  and  the  com¬ 
plementary  set  of  stresses  and  strains  is  determined  by  subincrementation: 
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After  determining  die  consistent  set  of  stresses  and  strains  in  the  various  subcells  of  the 
unit  cell,  the  material  properties  and  stresses  within  the  various  subcells  of  the  unit  cell  are 
updated. 

4.3.3  Residual  Calculation  (Step  3) 

The  difference  between  the  tangential  estimates  of  stresses  and  strains,  and  the 
consistent  stresses  and  strains  results  in  violation  of  some  of  the  constraints  at  the  micro¬ 
model  level  and  at  the  sublaminate  model  level.  The  process  of  residual  calculation  entails 
determining  the  amount  of  imbalance  in  those  constraints  which  are  violated,  and  calcu¬ 
lating  the  corrective  quantities  needed  to  remove  the  imbalance.  The  process  of  residual 
calculation  is  shown  in  Fig.  4.6;  it  starts  at  the  micro-model  level  and  progresses  upwards 
to  the  sublaminate  model  level. 


At  the  micro-model  level,  within  element  B  (matrix  subcell  mB  )  the  difference 
between  the  tangential  estimates  of  stresses  and  the  consistent  set  of  stresses  defines  a 
stress  residual  {<5o}mB .  Since  the  in-plane  strains  in  material  element  B  are  in  accordance 

with  constraints  (SI),  only  the  residual  {<5ob}mi  pertaining  to  out-of-plane  stresses  need  be 
considered,  and  is  given  as 

M. s  KL  =  {4  -  M  (410) 

JD  tug  ms  lug 

Stress 

Residual 

Similarly,  within  matrix  subcell  mx ,  the  difference  between  the  consistent  strains  and 
tangential  estimates  define  a  strain  residual  given  by: 


K}„ "  (4m  -  {4m 

mA  mA  mA 

Strain 
Residual 

The  strain  residual  within  matrix  subcell  mx  translates  into  a  strain  residual  within 


(4.11) 


material  element  A  as 
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Strain 

Residual 

The  strain  residual  in  material  element  A  is  the  amount  of  strain  mismatch  between  the 
material  elements  A  and  B,  which  ought  to  be  removed  to  satisfy  constraint  (Ml).  This 
would  require  application  of  a  corrective  stress  to  element  A  given  by 

wA-fe}A=  HfeL  (4-i3> 

Stress 

Correction 

This  stress  correction  would  have  to  be  applied  to  element  A  in  order  to  satisfy  constraint 
(Ml),  but  the  application  of  this  stress  correction  is  deferred  for  the  time  being. 

In  material  element  A  the  stress  correction  terms  {do0}A  pertaining  to  out-of-plane 

stresses  are  calculated  from  {<krr}A  by  selecting  only  the  out-of-plane  stress  terms. 

The  out-of-plane  stresses  within  the  lamina  are  now  in  violation  of  constraint  (SI), 
and  the  amount  of  imbalance  from  the  initial  tangential  estimate  is  given  by 

{<5o.}  ■  {do.}  =  -  ( W,{6o.}  +  !T„{<5o.}  )  <4-14) 

Stress 

Residual 

The  out-of-plane  stress  residual  in  the  lamina  would  result  in  an  out-of-plane  strain 
correction  term,  given  by 

{4— ■ tef =  [  **]  {£}*  (415) 

Strain 

Correction 

where  [p*l  is  the  partially  inverted  lamina  tangent  compliance.  The  effect  of  the  lamina 


strain  correction  term  calculated  in  Eq.  (4.15)  is  to  be  carried  back  to  the  fiber  and  the 
matrix  phases,  but  its  application  is  deferred  until  the  next  step. 


Hie  strain  correction  contributions  from  all  the  laminae  within  the  sublaminate 
are  summed  up  to  yield  a  strain  residual  at  the  sublaminate  level,  given  by 

M  - 1  (t)M‘  (4,6) 

Strain  k  =  1 
Residual 

Corresponding  to  the  sublaminate  strain  residual,  an  out-of-plane  stress  correction  term 
is  calculated  as 

M— fel“l>]fe}  (4-17> 

Stress 

Correction 

This  out-of-plane  stress  correction  term  is  applied  to  the  sublaminate  to  remove  the 
strain  residual  in  the  sublaminate.  This  is  explained  in  the  next  step  of  the  stress  update 
procedure. 

4.3.4  Residual  Application  (Step  4) 

Hie  schematic  algorithm  for  application  of  residuals  to  satisfy  the  out-of-balance 
constraints  is  shown  in  Fig.  4.7.  The  stress  correction  term  determined  in  Eq.  (4.17)  is  applied 
to  the  sublaminate,  and  its  effect  is  transferred  to  the  micro-model  level.  Hie  procedure  is 
essentially  die  same  as  that  followed  in  the  Thngent  Update  (Step  1)  with  a  slight  modifica¬ 
tion:  the  residuals  calculated  in  the  preceding  step  of  die  Residual  Calculation,  and  whose 
application  was  deferred,  are  also  taken  into  account  while  distributing  die  residuals. 

4.3.5  Updated  Stress  Calculation  (Step  5) 

After  all  the  constraints  are  satisfied  to  within  acceptable  tolerances,  the  effective 
stress  at  the  material  point  is  calculated  by  using  die  stress  homogenization  relations  used  in 
die  definition  of  the  micro-model  and  the  sublaminate  model  The  procedure  employed  is 
shown  schematically  in  Fig.  4.8.  After  the  updated  stresses  have  been  calculated,  this  infor¬ 
mation  is  passed  on  to  the  main  analysis  program. 
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4.3.6  Overall  Stress  Update  Algorithm 

The  overall  stress  update  algorithm  consists  of  the  five  previously  described  steps; 
step  5  is  the  last  step  in  the  procedure,  and  is  taken  only  when  the  constraints  are  within 
specified  tolerances.  The  check  for  satisfaction  of  constraints  is  made  after  Step  3  (Resid¬ 
ual  Calculation).  If  the  constraints  are  not  satisfied,  some  of  the  steps  have  to  be  repeated. 
The  flowchart  of  the  stress  update  procedure  is  shown  in  Fig.  4.9. 

As  is  shown  in  Fig.  4.9,  the  prescribed  strain  increment  Ac  at  a  material  point  is 
received  as  input  from  the  main  structural  analysis  program.  Step  1  (Thngent  Update), 
Step  2  (Constitutive  Properties  Update),  and  Step  3  (Residual  Calculation)  are  per¬ 
formed,  following  which  a  decision  is  made  about  the  convergence  of  residuals.  In  order 
to  ensure  that  the  out-of-plane  strain  residual  {<5ca}  is  small,  the  following  convergence 
criterion  is  used, 

^  .  <,  Tolerance  (4-18) 

U  ArQ 

If  the  above  convergence  criterion  is  satisfied,  then  the  updated  stresses  are  calculated 
and  the  results  are  conveyed  to  the  main  program.  If  the  convergence  criterion  is  not 
satisfied,  then  Step  4  (Residual  Redistribution)  is  carried  out,  and  Step  2  and  Step  3  are 
repeated. 

4.4  Comments  on  the  Stress  Update  Procedure 

The  procedure  presented  in  this  chapter  for  obtaining  the  updated  stresses  was  found 
to  be  quite  accurate.  After  the  stresses  are  updated  all  the  constraints  are  satisfied,  with  some 
error  remaining  in  die  constitutive  relations.  This  is  partly  the  result  of  the  subincrementa¬ 
tion  procedure  adopted  here,  and  also  because  of  the  fact  that  the  relaxation  type  scheme 
adopted  here  is  iterative  in  nature.  Subincrementation  is  an  approximate  solution  scheme, 
and  therefore  the  results  obtained  contain  some  error.  More  refined  integration  schemes  can 
be  used  for  integration  of  constitutive  relations;  however,  this  would  result  in  making  the 
stress  updating  scheme  more  complex. 
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Fig.  4.1  Incremental  solution  of  nonlinear  structural  systems 
with  Newton-Raphson  Iterations. 


Fig.  4.2  Schematic  view  of  stress  update  Illustrating  the  difference  between 
true  stress  Increment  bJ5  and  the  tangential  estimate  hJj 
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Fig.  4.3  Schematic  statement  of  the  stress  update  problem. 


82 


Fifl.  4.4  Schematic  algorithm  tor  Tangent  Update (•,  o)  procedure,  Step  (1) 
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Fig.  4.5  Schematic  algorithm  for  Constitutive  Relations  Update! «  o),  Step  (2) 
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Fig.  4.6  Schematic  algorithm  for  Residual  Calculation,  Step  (3). 
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Fig.  4.7  Schematic  algorithm  for  Residual  Application  (•),  Step(4). 
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Fig.  4.9  Flowchart  for  Stress  Update  at  a  material  point 


CHAPTER  5 


RESPONSE  OF  LAMINAE  AND  LAMINATES 
TO  UNIFORM  LOADINGS 


In  this  chapter  some  applications  are  presented  with  the  objective  of  evaluating 
the  accuracy  of  the  modelling  procedure  and  demonstrating  its  potential.  It  is  desirable  to 
separate  the  two  issues  of  structural  modelling  and  modelling  of  the  material  behaviour. 
For  this  reason,  the  performance  of  the  material  model  is  evaluated  in  modelling  and 
predicting  the  response  of  fiber-reinforced  laminae  and  laminates  under  homogeneous 
loading  conditions,  i.e.  when  there  is  no  spatial  variation  in  the  effective  stresses  and 
strains.  Results  are  presented  for: 

•  Response  of  various  types  of  laminae,  comprised  of  different  fiber/ 
matrix  combinations,  subjected  to  off-axis  loading,  and  laminates 
subjected  to  uniaxial  loading.  Only  the  effects  of  matrix  stress- 
strain  nonlinearity  are  taken  into  account,  since  the  applied  loading 
is  not  high  enough  to  cause  damage. 

•  Prediction  of  the  off-axis  strength  of  various  types  of  laminae 

•  Initial  yield  surfaces  of  composite  laminae  subjected  to  different 
stress  fields. 

•  Stress-strain  response  of  laminates  to  the  point  of  ultimate  failure. 

In  this  case  both  the  effects  of  matrix  nonlinearity  and  other  damage 
mechanisms  are  taken  into  account.  Micromechanical  criteria  are 
employed  for  damage  detection  and  incorporating  the  effect  of 
damage. 
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5.1  Nonlinear  Response  of  Laminae  and  Laminates 


Aboudi  [53],  using  the  Method  of  Cells,  modelled  the  nonlinear  response  of  vari¬ 
ous  types  of  laminae  and  laminates,  comprised  of  different  fiber/matrix  combinations.  A 
standard  Ramberg-Osgood  description  was  adopted  for  the  matrix  material.  The  stress- 
strain  response  of  laminae  and  laminates  predicted  by  Aboudi  showed  fairly  good  agree¬ 
ment  with  experimentally  observed  results.  Here,  the  proposed  material  modelling  proce¬ 
dure  is  used  to  predict  the  responses  of  the  laminae  and  laminates  studied  by  Aboudi  [53], 
using  the  same  elastic  properties  and  constitutive  description  for  the  constituent  fiber  and 
matrix  phases  as  used  by  him. 

5.1.1  Response  of  Boron/Epoxy  Laminae  and  Laminates 

Fig.  5.1  shows  the  response  of  a  Boron/epoxy  lamina  in  shear.  It  is  worth  mention¬ 
ing  that  the  shear  strain  along  the  abscissa  is  the  tensorial  shear  strain  and  not  the 

engineering  shear  strain,  commonly  denoted  by  yn .  The  elastic  properties  of  the  Boron 

fibers  and  epoxy  matrix  are  given  in  Ihble  5.1,  and  are  the  same  as  the  ones  used  by 
Aboudi  [53]. 

As  can  be  seen  from  Fig  5.1,  the  results  from  the  present  formulation  and  Aboudi’s 
Method  of  Cells  are  in  excellent  agreement.  There  is  discrepancy  between  the  theoretical 
predictions  and  experimental  data  in  the  vicinity  of  en  =  0.75 % ;  Aboudi  attributed  this 

to  error  in  experimental  data  obtained  from  the  15°  off-axis  coupon. 

Figs.  5.2  and  5.3  show  the  response  of  off-axis  coupons  (0  =  30°,  45°) ,  subjected  to 
uniaxial  tension;  results  are  shown  for  analytical  predictions  from  Aboudi’s  Method  of 
Cells,  and  the  proposed  material  modelling  procedure,  along  with  experimental  data  from 
Cole  and  Pipes  [78].  The  predictions  form  both  of  the  analytical  procedures  agree  well 
with  the  experimental  data.  The  predictions  of  ?n  from  both  analytical  procedures  show 
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excellent  agreement,  but  the  *22  strains  predicted  from  the  present  formulation  are  slight¬ 
ly  less  than  the  ones  predicted  by  Aboudi. 

Figs.  5.4  and  5.5  show  the  response  of  symmetric  angle-ply  laminates  to  uniaxial 
loading.  Results  are  presented  for  [+  30/  -  30js  and  [+  45/  -  45}s  laminates.  In  both  cases 
the  predictions  from  the  present  method  are  in  excellent  agreement  with  results  from 
Aboudi  [53].  However,  for  the  [+  45/  -  45fe  laminate  the  theoretical  predictions  show  a 
deviation  from  the  experimentally  observed  results  of  Pipes,  Kaminski  and  Pagano  [80]. 

5.1.2  Response  of  E-Glass/Epoxy  Laminates 

Fig.  5.6  presents  the  shear  stress-strain  response  of  E-glass/epoxy  lamina.  The 
approximated  shear  stress-strain  response  attributed  to  Lifshitz  was  obtained  by  Aboudi 
[53]  by  integrating  the  following  expression  for  the  axial  tangent  shear  modulus  Grn  of 

the  lamina,  reported  by  Lifshitz  [82]: 

GnT  -  6  -  692?^  +  32672^  -  539248^  (5.1) 

Also  shown  in  Fig.  5.6  is  the  shear  stress-strain  response  of  the  lamina  as  modelled  by 
Aboudi’s  Method  of  Cells  and  by  the  present  material  modelling  procedure.  The  assumed 
properties  of  E-glass  fibers  and  epoxy  matrix  are  given  in  Thble  5.2,  and  are  consistent 
with  the  properties  used  by  Aboudi  [53].  As  can  be  seen  from  Fig.  5.6,  the  results  from  the 
two  analytical  modelling  procedures  are  in  excellent  agreement  with  each  other  and  with 
the  derived  experimental  results. 

Figs.  5.7  and  5.8  show  the  response  of  [+  30/ -30]s  and  [  +  41/  -  41]s  laminates 
when  subjected  to  uniaxial  tension.  The  curves  representing  the  experimental  results  of 
Lifshitz  [82]  are  actually  he  upper  and  lower  bounds  of  the  experimentally  observed 
results.  The  response  of  the  [+  41/  -41fe  laminate  shows  significant  nonlinearity.  The 
agreement  between  the  analytical  predictions  and  the  experimental  results  is  very  good. 
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5.1.3  Response  of  Graphite/Polyimide  Laminae 

Fig.  5.9  shows  the  response  of  graphite/polyimide  off-axis  coupons.  Experimental 
results  were  presented  by  Pindera  and  Herakovich  [83].  Results  from  Aboudi’s  Method 
of  Cells  and  the  present  material  modelling  procedure  are  also  presented,  and  are  seen  to 
be  in  excellent  agreement  with  each  other.  The  agreement  between  analytical  predictions 
and  experimental  data  is  fair.  The  properties  of  graphite  fibers  and  polyimide  matrix  used 
in  the  analytical  predictions  are  given  in  Thble  53. 

5.2  Off-Axis  Strength  of  Composite  Laminae 

Off-axis  strengths  of  various  types  of  fiber-reinforced  composite  laminae  were 
predicted  by  Aboudi  [52],  using  the  Method  of  Cells.  Failure  of  the  composite  lamina  was 
predicted  by  using  micro-failure  criteria  for  the  fiber  and  matrix  phases.  The  lamina  was 
subjected  to  off-axis  loading,  and  stresses  within  the  various  subcells  of  the  unit  cell 
representing  the  lamina  were  monitored.  A  maximum  stress  failure  criterion  was  applied 
to  the  constituent  phases,  to  define  the  failure  of  the  lamina.  Aboudi’s  analytical  results 
were  in  good  agreement  with  experimental  results. 

In  this  section,  predictions  of  off-axis  strength  from  the  present  method  will  be 
compared  with  the  results  from  Aboudi’s  work  [52].  It  was  assumed  by  Aboudi  that  an 
off-axis  lamina  loaded  in  uniaxial  tension  can  foil  in  one  of  the  following  three  modes: 
Fiber  fracture.  Matrix  Shearing,  and  Matrix  Cracking  in  transverse  tension.  The  micro- 
failure  criteria  determining  these  modes  of  failure  have  already  been  presented  in  Chap¬ 
ter  TWo . 

5.2.1  Strength  of  Boron/Epoxy  Lamina 

Fig.  5.10  shows  the  off-axis  strength  of  a  boron/epoxy  lamina  with  variation  in 
fiber  orientation  angle.  Results  are  presented  for  analytical  predictions  from  Aboudi’s 
[52]  Method  of  Cells  and  the  present  modelling  procedure.  Experimental  results  from 
Pipes  and  Cole  [77]  are  also  presented.  The  elastic  properties  and  strengths  of  boron 
fibers  and  epoxy  matrix  are  given  in  foble  5.4.  It  is  seen  that  the  predictions  of  both 
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analytical  procedures  agree  very  well  over  the  full  range  of  fiber  orientation  angle.  The 
present  model  not  only  predicts  the  off-axis  strength  with  good  accuracy,  but  also  success¬ 
fully  delineates  the  transition  from  fiber  fracture  mode  to  matrix  shearing  mode  and 
matrix  tensile  failure  mode.  The  agreement  between  analytical  predictions  and  exper¬ 
imental  results  is  quite  good. 

5.2.2  Additional  Examples  of  Off-Axis  Strength  of  Laminae 

The  off-axis  strength  of  a  graphite/polyimide  lamina  is  shown  in  Fig.  5.11.  The 
experimental  results  were  reported  by  Pindera  and  Herakovich  [83].  The  properties  of 
graphite  fibers  and  epoxy  matrix  used  in  the  analyses  are  listed  in  Ihble  5.5.  Excellent 
agreement  between  die  results  from  Aboudi  [52]  and  the  present  method  is  seen  again. 

Fig.  5. 12  shows  the  off-axis  strength  of  AS/3501  graphite/epoxy  lamina  with  varia¬ 
tion  in  fiber  orientation  angle.  The  experimental  results  were  reported  by  Ifcai  and  Hahn 
[6].  The  properties  of  graphite  fibers  and  epoxy  matrix,  used  in  the  micromechanics  based 
predictions  of  strength  are  listed  in  Ihble  5.6. 

The  off-axis  strength  of  a  Kevlar/epoxy  lamina  is  shown  in  Fig.  5.13.  The  exper¬ 
imental  results  are  from  the  work  of  Pindera  et  al.  [84],  and  the  elastic  properties  of  Kevlar 
fibers  and  epoxy  matrix  used  in  micromechanical  analysis  are  given  in  Ihble  5.7. 

Fig.  5.14  shows  the  off-axis  strength  of  E-glass/epoxy  lamina;  the  properties  of 
E-glass  fibers  and  epoxy  matrix  are  presented  in  Ihble  5.8.  It  is  interesting  to  note  that  in 
this  case  the  zone  in  which  the  fiber  fracture  mode  occurs  is  very  small;  the  transition  from 
fiber  failure  mode  to  the  matrix  shear  failure  mode  occurs  at  about  an  off-axis  angle 
0*1°.  This  is  due  to  the  relatively  low  shear  strength  of  the  epoxy  matrix. 

Fig.  5.15  shows  the  variation  in  off-axis  strength  of  a  boron/aluminum  lamina  as 
the  fiber  orientation  angle  changes.  The  experimental  results  presented  are  from  the  work 
of  Becker  et  al.  [86].  The  elastic  properties  of  the  boron  fibers  and  aluminum  matrix  are 
given  in  Ihble  5.8.  The  three  distinct  modes  of  failure,  and  the  transitions  from  one  mode 
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of  failure  to  another  are  dearly  visible  in  Fig.  5.15.  The  region  in  which  the  fiber  fracture 
mock  governs  is  significantly  large  compared  to  the  previously  presented  examples, 
because  the  shear  strength  of  the  aluminum  matrix  is  relatively  high. 

5.3  Initial  Yield  Surfaces  of 
Metal  Matrix  Composites 

Metal  matrix  composites  are  increasingly  being  used  in  aerospace  applications, 
where  the  use  of  stiff  fibers,  such  as  boron  or  graphite  in  a  light-weight  metallic  matrix 
like  aluminum,  results  in  a  light,  yet  stiff  and  strong  material.  The  fibers  in  a  composite 
can  usually  be  assumed  to  be  elastic;  and  the  yielding  observed  in  metal  matrix  composites 
is  attributable  to  yielding  of  the  metallic  matrix.  It  is  of  interest  to  see  if  the  present 
micro-model  can  predict  the  initiation  of  yielding  in  metal  matrix  composites;  therefore 
in  this  section,  results  for  initial  yield  surfaces  of  metal  matrix  composites  under  different 
stress  states  are  presented. 

Results  for  initial  yield  surfaces  of  fiber-reinforced  laminae  and  laminates  under 
different  stress  states  were  presented  by  Dvorak  et  al.  [40],  using  the  finite  element  meth¬ 
od  to  model  a  representative  volume  element  of  a  composite  lamina.  Pindera  and  Aboudi 
[87]  used  the  Method  of  Cells  to  obtain  the  initial  yield  surfaces  of  the  composite  laminae 
studied  by  Dvorak  et  al.  In  this  section  the  results  of  Dvorak  et  al.  and  Pindera  and  Aboudi 
are  compared  to  results  obtained  from  the  proposed  modelling  procedure. 

To  predict  the  onset  of  yielding  in  matrix  material,  the  von  Mises  yield  criterion  is 
used,  which  is  described  as  follows.  Yielding  is  supposed  to  occur  if  the  following  condi¬ 
tion  is  satisfied. 


F(<7y)  £  0 


where 


(52) 

(53) 
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tg  is  the  square  root  of  the  second  invariant  of  the  deviator  stress  tensor  sy  and  is 

the  yield  stress  of  the  matrix  in  uniaxial  tension.  Eq.  (5.3)  may  be  expressed  in  terms  of 
stresses  as 


F(Oy)  =  |  j  Jjtfu-ofc)2  +  (^n  ~  a33}2  + 

+  3(o^  +  cfc  +  ©§3)  J  _  a£m) 


(5.4) 


The  above  yield  criterion  is  applied  to  both  matrix  subcells,  dia  and  ms ,  and 

yielding  is  assumed  to  occur  if  the  above  yield  criterion  is  satisfied  in  either  of  the  two 
matrix  subcells.  Specifically  stated,  yielding  is  assumed  if  either  of  the  two  following 
criteria  are  satisfied. 


F(4"°)  a  0 

(5.5a-b) 

a  0 

Fig.  5.16  shows  the  initial  yield  surface  of  a  boron/aluminum  lamina  in  [022,  <733) 

stress  space;  the  stresses  are  normalized  with  respect  to  the  yield  stress  of  the  matrix  oJP . 
The  properties  of  the  boron  fibers  and  aluminum  matrix  are  given  in  Tkble  5.9.  The  results 
from  the  present  material  model  show  very  good  agreement  with  results  from  Aboudi’s 
Method  of  Cells  and  the  finite  elements  based  results  of  Dvorak  et  al.  However,  the  yield 
surface  predicted  by  Aboudi’s  Method  of  Cells  contains  comers  which  are  not  predicted  by 
the  detailed  finite  element  approach. 

Fig.  5.17  shows  the  initial  yield  surface  of  the  same  boron/aluminum  lamina  in 
(oil,  <^2)  stress  space.  The  results  from  the  present  material  model  agree  very  well  with 
the  results  from  Aboudi’s  Method  of  Cells  and  the  finite  element  based  results. 
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Fig.  5.18  shows  the  initial  yield  surface  of  the  boron/aluminum  lamina  in 
(oil,  ^2  =  #33)  stress  space.  It  is  noticed  that  in  the  vicinity  of  the  hydrostatic  stress  state, 

both  the  Method  of  Cells  and  the  present  micro-model  perform  poorly  and  predict  yield 
stress  states  which  are  significantly  higher  than  the  finite  element  based  results.  However, 
the  results  from  both  models  improve  as  the  imposed  stress  states  move  away  from  the 
hydrostatic  stress  state. 

5.4  Response  and  Ultimate  Strength  of 
Boron/Epoxy  Laminates 

Petit  and  Whddoups  [61]  presented  a  method  for  modelling  the  nonlinear  response 
of  laminated  composites  in  which  the  nonlinear  stress-strain  behaviour  was  taken  into 
account,  but  no  interaction  between  stresses  was  recognized.  Furthermore  the  bimodu¬ 
larity  of  the  lamina  in  tension  and  compression  was  taken  into  account.  The  laminate 
response  was  derived  by  assembling  the  laminae  comprising  the  laminate  within  die 
framework  of  Classical  Laminated  Plate  Theory.  A  maximum  strain  criterion  was 
employed  to  detect  the  onset  of  failure  in  a  lamina,  and  the  lamina  was  then  gradually 
unloaded  by  giving  it  a  negative  tangent  stiffness  in  the  loading  direction  being  considered 
(longitudinal,  transverse  or  shear). 

5.4.1  Micromechanical  Modelling  of  Boron/Epoxy  Lamina 

The  response  of  the  boron/epoxy  lamina  as  reported  by  Petit  and  Waddoups  [61]  is 
now  modelled  using  the  proposed  material  model.  Fig.  5.19  shows  the  shear  response  of 
the  boron/epoxy  lamina;  also  shown  is  the  response  predicted  by  the  material  model.  The 
properties  of  the  constituent  boron  fibers  and  epoxy  matrix  were  inferred  from  lamina 
responses  under  various  loading  conditions  (Fig.  5.19-5.23),  and  are  listed  in  Thble  5.10. 
It  should  be  mentioned  that  the  properties  of  the  fibers  and  the  epoxy  matrix  are  to  be 
used  in  the  “Modified  Ramberg-Osgood  Relations”  (Eqs.  (2.23-2.29),  Chapter  2).  Modi¬ 
fied  Ramberg-Osgood  relations  are  used  in  order  to  eliminate  excessive  nonlinearity  in 
transverse  direction  loading  which  is  not  observed  in  experimental  results. 
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5.4.2  Failure  Criteria  and  Damage  Modelling 

Hie  employed  micro-failure  criteria  for  the  various  failure  modes  have  already 

been  presented  in  chapter  2  (Section  2.6).  The  procedure  for  modelling  the  effects  of 
damage  have  also  been  described  in  chapter  2  (Section  2.7).  The  stiffness  reduction 
factors  used  in  the  analysis  are  listed  in  Thble  5.10. 

Results  for  the  response  of  boron/epoxy  laminates  are  presented  in  the  next  section. 

5.4.3  Response  of  Boron/Epoxy  Laminates  in  Tension 

Figs.  5.24-5.27  show  the  responses  of  boron/epoxy  laminates  subjected  to  uniaxial 
tension.  Results  are  presented  for  the  proposed  material  model,  analytical  predictions, 
and  the  experimental  results  of  Petit  and  Waddoups  [61]. 

Fig  5.24  shows  the  response  of  a  [0/90]  s  laminate  in  tension.  The  response  pre¬ 
dicted  by  the  proposed  material  model  agrees  well  with  the  experimental  results  and  the 
analytical  results  of  Petit  and  Widdoups.  According  to  the  present  analysis  procedure, 
tensile  matrix  cracking  in  the  transverse  direction  takes  place  in  the  90°  plies  at  the  stress 
level  of  70  ksi,  resulting  in  some  loss  of  stiffness  and  load  redistribution.  This  is  consistent 
with  the  analytical  results  of  Petit  and  Waddoups.  The  ultimate  failure  of  the  laminate 
occurred  by  fiber  fracture  in  the  0°  plies  at  92  ksi;  the  ultimate  strain  being  0.6  %. 
Experimental  results  indicate  failure  at  100  ksi,  with  a  strain  at  failure  of  0.71%. 

The  response  of  a  [+  30/  -  30]s  laminate  is  shown  in  Fig.  5.25.  The  present  proce¬ 
dure  indicates  some  nonlinearity  in  response,  but  not  as  much  as  is  indicated  by  Petit  and 
Whddoups.  Also,  shear  failure  of  the  matrix  in  the  plies  is  predicted,  whereas  the  analysis 
conducted  by  Petit  and  Waddoups  indicated  the  mode  of  failure  to  be  compressive  failure 
of  the  matrix  in  the  transverse  direction. 

Fig.  5.26  presents  the  response  of  a  [  +  60/  -  60]  s  laminate.  The  indicated  mode  of 
failure  from  the  present  analysis  is  tensile  matrix  cracking  in  the  transverse  direction  at  15 
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ksi,  the  ultimate  strain  being  0.6%.  The  experimental  results  indicate  an  ultimate  stress  of 
17.5  ksi,  and  strain  at  failure  of  about  0.8%. 

Fig.  5.27  shows  the  behaviour  of  a  [+  45/-45/903]s  laminate.  The  indicated  mode 
of  failure  is  tensile  matrix  cracking  of  the  90°  plies  in  the  transverse  direction,  followed 
by  matrix  shearing  failure  in  the  45°  plies.  The  predicted  ultimate  strength  is  22  ksi  at  the 
ultimate  strain  of  0.44%.  The  experimental  results  indicated  the  mode  of  failure  to  be 
matrix  cracking  of  the  90°  plies  in  the  transverse  direction  at  20  ksi. 

5.4.4  Response  of  Boron/Epoxy  Laminates  in  Compression 

The  response  of  a  [0/90]s  laminate  in  compression  is  shown  in  Fig.  5.28.  The 

presentanalysispredictedfailureofthe  0°  plies  by  compressive  failure  of  the  matrix  in  the 
transverse  direction  at  about  230  ksi.  This  resulted  in  a  slight  loss  of  stiffness,  but  did  not 
cause  complete  loss  of  load  carrying  capacity.  The  laminate  finally  failed  at  275  ksi  by 
micro-buckling  of  the  0°  plies.  Experimental  results  indicate  that  the  laminate  failed  at 
255  ksi  by  failure  of  the  0°  plies. 

Compression  response  of  a  [  +  20/  -  20]s  laminate  is  shown  in  Fig.  5.29.  No  nonlin¬ 
earity  is  observed  in  the  laminate  response,  and  the  analytical  procedure  predicts  laminate 
failure  by  tensile  failure  of  the  plies  in  the  transverse  direction  at  117  ksi.  The  [+30/-30]s 

laminate  fails  in  a  similar  fashion  (Fig.  530).  The  shapes  of  the  analytical  response  curves 
and  the  experimental  results  show  very  good  agreement  Experimental  results  indicated  the 
failure  stress  to  be  44  ksi,  as  compared  to  50  ksi  predicted  by  Petit  and  Whddoups,  and  60  ksi 
predicted  by  the  present  material  model  and  micromechanical  failure  criteria. 

The  predicted  response  and  experimental  results  for  a  [+  60/  -  60]  s  laminate  are 

given  in  Fig.  531.  A  significant  amount  of  nonlinearity  is  observable  in  the  response.  The 
present  analysis  indicates  failure  at  41  ksi  due  to  compressive  failure  of  the  matrix  in  the 
transverse  direction,  which  is  consistent  with  the  experimental  results  and  the  analytical 
predictions  of  Petit  and  Waddoups.  The  present  analysis,  however,  predicts  an  elastic- 
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plastic  type  response,  whereas  the  analytical  results  of  Petit  and  Waddoups  show  a  more 
gradual  transition  to  the  point  when  the  laminate  loses  all  stiffness. 

5.4.5  Comments  on  Modelling  of  Laminate 

Stress-Strain  Response 

Modelling  the  response  of  fiber-reinforced  laminated  composites  up  to  the  point 
of  ultimate  strength  is  a  challenging  problem,  made  difficult  by  a  number  of  factors.  The 
experimentally  observed  ultimate  strengths  of  unidirectional  composites  show  a  consider¬ 
able  amount  of  scatter,  indicating  that  statistically  reliable  properties  should  be  used  in 
the  analysis  to  obtain  better  statistical  correlation  between  the  predicted  laminate 
strength  and  the  experimentally  observed  values.  Keeping  in  view  the  observed  scatter  in 
the  experimentally  obtained  ultimate  strength  values  of  laminated  composites,  the  true 
measure  of  the  performance  of  an  analysis  procedure  for  predicting  the  strength  of  com¬ 
posites  would  be  its  better  statistical  correlation  with  experimental  results. 

The  choice  of  the  failure  criteria  adopted  in  the  analysis  procedure  can  significantly 
affect  the  quality  and  accuracy  of  the  predicted  results.  The  failure  mechanisms  exhibited  by 
composites  are  many  and  quite  complex..  Although  numerous  failure  criteria  have  been 
proposed  for  strength  prediction  of  composites,  much  work  needs  to  be  done  in  this  area  to 
develop  criteria  which  are  accurate  and  take  into  account  the  complex  mechanisms  of  failure 
associated  with  composites. 

Finally,  the  manner  in  which  the  effects  of  damage  are  accounted  for  in  the  lami¬ 
nate  analysis  also  significantly  affects  the  predicted  response  of  the  laminate.  Damage 
mechanisms  tend  to  dissipate  the  stored  strain  energy  in  the  laminate  and  result  in  stif¬ 
fness  degradation.  The  manner  in  which  the  stiffness  degradation  and  strain  energy  dissi¬ 
pation  is  modelled  would  influence  the  calculated  response  of  the  laminate.  The  field  of 
Damage  Mechanics  aims  to  provide  damage  evolution  laws  and  framework  for  incorporat¬ 
ing  the  effects  of  damage.  It  is  felt  that  incorporating  a  suitable  damage  theory  would 
result  in  improved  analytical  predictions  of  laminate  response. 
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E 

V 

to 

a 

b 

n 

(GPa) 

(MPa) 

Boron 

413.6 

0.15 

Epoxy 

5.2 

0.35 

157.6 

0 

5.74 

4 

Table  5.1  Elastic  properties  of  boron  fibers  and  epoxy  matrix,  along  with  Ramberg- 
Osgood  parameters  for  epoxy  matrix.  Fiber  volume  fraction  =  0.5 


Fig.  5.1  Shear  stress-strain  response  of  boronfepoxy  lamina  as  predicted  by 

Aboudi's  Method  of  Cells  and  the  proposed  material  model,  along  with 
experimental  results  of  Pipes  and  Cole. 
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i.  5.2  Uniaxial  stress-strain  response  of  30°  off-axis  boron! epoxy  lamina  as 
predicted  by  Aboudi's  Method  of  Cells  and  the  proposed  material 
model,  along  with  experimental  results  of  Cole  and  Pipes. 


Fig.  5.3  Uniaxial  stress-strain  response  of  45 0  off-axis  boron/epoxy  lamina 
as  predicted  by  Aboudi’s  Method  of  Cells  and  the  proposed  material 
model,  along  with  experimental  results  of  Cole  and  Pipes. 
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5.4  Uniaxial  Stress-Strain  Response  of  [ +301-30] s  Boron! epoxy  laminate, 
as  predicted  by  Aboudi's  method  of  cells  and  the  proposed  material 
model,  along  with  experimental  results  of  Petit. 


Fig.  5.5  Uniaxial  Stress-Strain  Response  of  [+451-45]$  Boron/epoxy  laminate, 
as  predicted  by  Aboudi's  method  of  cells  and  the  proposed  material 
model,  along  with  experimental  results  of  Pipes,  Kaminski  and  Pagano. 
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(MPa) 


E 

V 

T0 

a 

b 

n 

(GPa) 

(MPa) 

E-Glass 

73 

0.22 

Epoxy 

6 

0.35 

106 

0 

3.33 

3 

Table  S.2  Elastic  properties  of  E-glass  fibers  and  epoxy  matrix,  along  with  Ramberg- 
Osgood  parameters  for  epoxy  matrix.  Fiber  volume  fraction  =  0.64 


Tensorial  Shear  Strain  €j2  (*) 


FiQ.  5.6  Shear  stress-strain  response  of  E-glasslepoxy  lamina  as  predicted  by 
Aboudl’s  Method  of  Cells  and  the  proposed  material  model,  along  with 
derived  experimental  results  of  Ufshltz  and  Gilat. 
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Fig.  5.7  Uniaxial  stress-strain  response  of  [+301-30]  s  E-glasslepoxy 

laminate,  as  predicted  by  Aboudi's  Method  of  Cells  and  the  proposed 
material  model,  along  with  experimental  results  of  Ufshitz. 
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Fig.  5.8  Uniaxial  stress-strain  response  of  [+451-45] s  E-glasslepoxy 

laminate,  as  predicted  by  Aboudi’s  Method  of  Cells  and  the  proposed 
material  model,  along  with  experimental  results  of  Ufshitz. 
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Ei 

Vi 

Ez 

Vzs 

G\z 

To  a 

b 

n 

(GPa) 

(GPa) 

(GPa) 

(MPa) 

Graphite 

222 

0.33 

29.5 

0.73 

24.1 

-  - 

— 

— 

Polyimide 

3.1 

0.39 

3.1 

0.39 

1.1 

186.9  0 

5.61 

4 

Table  5.3  Elastic  properties  of  graphite  fibers  and  polylmide  matrix,  along  with  Ramberg- 
Osgood  parameters  for  polyimide  matrix.  Fiber  volume  fraction  ■  0.61 


°n 


Fig.  5.9  Uniaxial  stress-strain  response  of  10°, 15°  and  30 0  off-axis  graphite i polyimide 
laminae,  as  predicted  by  Aboudl’s  Method  of  Cells  and  the  proposed  material 
model,  along  with  experimental  results  of  Pindera  and  Herakovich. 
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E 

(GPa) 

V 

(MPa) 

(MPa) 

s(m) 

(MPa) 

Boron 

400 

0.2 

2566 

— 

Epoxy 

5.2 

0.35 

— 

93 

87 

Table  5.4  Elastic  properties  and  strengths  of  boron  fibers  and  epoxy 
matrix.  Fiber  volume  fraction  =  0.5 


$  (Degrees) 

Fig.  5.10  Off-axis  strength  of  boron /epoxy  lamina  as  predicted  by  Aboudi's  Method  of  Cells 
and  the  proposed  material  model,  along  with  the  experimental  results  of  Pipes  and 
Cole. 
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Vi 

e2 

V23 

Gi  2 

*r 

if0 

5(m) 

( GPa ) 

(GPa) 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

Graphite 

222 

0.33 

29.5 

0.73 

24.1 

2517 

— 

— 

Polyimide 

3.1 

0.39 

3.1 

0.39 

1.1 

— 

69 

71 

Table  S.S  Elastic  properties  and  strengths  of  graphite  fibers  and  polyimlde  matrix. 
Fiber  volume  fraction  =  0.61 


0  (Degrees) 

FIq.  5.11  Off-axis  strength  of  graphitelpolylmlde  lamina  as  predicted  by  Aboudl's  Method  of 
Cells  and  the  proposed  material  model,  along  with  the  experimental  results  of 
Pindera  and  Herakovlch. 


107 


Ei 

Vi 

e2 

V23 

G\2 

X[f) 

if0 

s(m) 

(GPa) 

(GPa) 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

AS/3501 

Graphite 

213.7 

0.2 

13.8 

0.25 

13.8 

2250 

— 

— 

Epoxy 

3.45 

0.35 

3.45 

0.35 

1.3 

—  ■ 

62.9 

108 

Table  5.6  Elastic  properties  and  strengths  of  graphite  fibers  and  epoxy  matrix. 
Fiber  volume  fraction  =  0.66 


9  (Degrees) 


Fig.  5.12  Off-axis  strength  of  AS/3501  graphlte/epoxy  lamina  as  predicted  by  Aboudi's 
Method  of  Cells  and  the  proposed  material  model,  along  with  the  experimental 
results  of  Tsai  and  Hahn. 
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E , 

Vi 

e2 

V23 

Gi2 

Xt(f) 

y(n>) 

•** 

s(m) 

(GPa) 

(GPa) 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

Kevlar 

124.1 

0.35 

4.1 

0.35 

2.9 

2031 

— 

— 

Epoxy 

3.45 

0.35 

3.45 

0.35 

1.3 

— 

62.9 

108 

Table  5.7  Elastic  properties  and  strengths  of  Kevlar  fibers  and  epoxy  matrix. 
Fiber  volume  fraction  =  0.55 


$  (Degrees) 

HO*  5.13  Off-axis  strength  of  Kevlar/ epoxy  lamina  as  predicted  by  Aboudl's  Method  of  Cells 

and  the  proposed  material  model,  along  with  the  experimental  results  of  Pindera  et  at. 


109 


E 

(GPa) 

V 

jf 

(MPa) 

if0 

(MPa) 

s(m) 

(MPa) 

E-Glas$ 

73 

0.22 

1990 

— 

— 

Epoxy 

3.45 

0.35 

— 

40.4 

45.5 

Table  5.8  Elastic  properties  and  strengths  of  E-glass  fibers  and  epoxy 
matrix.  Fiber  volume  fraction  =  0.60 


0  (Degrees) 

Fig.  5.14  Off-axis  strength  of  E-glass/epoxy  lamina  as  predicted  by  Aboudi’s  Method  of 
Cells  and  the  proposed  material  model,  along  with  the  experimental  results  of 
Hashln  and  Rotem. 
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E 

V 

Af 

if” 

s(m) 

(GPa) 

(MPa) 

(MPa) 

(MPa) 

Boron 

400 

0.2 

2313.8 

— 

72.5 

0.33 

— 

155.6 

167.4 

Table  5.8  Elastic  properties  and  strengths  of  boron  fibers  and  6061-0 
aluminum  matrix.  Fiber  volume  fraction  =  0.46 


$  (Degrees) 

Fig.  5.15  Off-axis  strength  of  boron/6061 -O  aluminum  lamina  as  predicted  by  Aboudi's 
Method  of  Cells  and  the  proposed  material  model,  along  with  the  experimental 
results  of  Becker  et  at. 
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V 

a  (m) 

°Y 

(GPa) 

(MPa) 

Boron 

413.7 

0.21 

— 

Aluminum 

72.5 

0.33 

262 

Table  5.9  Elastic  properties  of  boron  fibers  and  aluminum  matrix  for 
calculation  of  Initial  yield  surfaces  of  composite  lamina. 
Fiber  volume  fraction  =  0.3 


Boron/Aluminum 

- F.  E.  Results.  Dvorak  et  ai.  [39] 

- Plndera  and  Aboudi  [87] 

-  Material  Model 


Fig.  5.16  Initial  yield  surface  of  a  boron! aluminum  composite  as  predicted 
by  Aboudi’ s  Method  of  Cells  and  the  proposed  material  model, 
along  with  finite  element  based  results  of  Dvorak  et  al. 
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Boron/ Aluminum 


022/0™ 

Fig.  5.17  Initial  yield  surface  of  a  boron! aluminum  composite  as  predicted 
by  Aboudl’s  Method  of  Cells  and  the  proposed  material  model, 
along  with  finite  element  based  results  of  Dvorak  et  al. 
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Boron/Aluminum 

F.  E.  Resuits,  Dvorak  et  al.  [40] 
Pindera  and  AboucU  [87] 

Material  Model 


a22/^ym)= 


Fig.  5.18  Initial  yield  surface  of  a  boron/aluminum  composite  as  predicted 
by  Aboudl’s  Method  of  Cells  and  the  proposed  material  model, 
along  with  finite  element  based  results  of  Dvorak  et  al. 
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to 
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b 

n 

x® 

Jm) 

*  C 

s(m) 

4 

(Ksi) 

(Ksi) 

(Ksi) 

(Ksi) 

(Ksi) 

(Ksi) 

Boron 

59985 

0.25 

372 

Epoxy 

1000 

0.35 

13.5 

0 

1.0 

4 

15.5 

50.0 

21.0 

Table  5.10  Elastic  properties  of  boron  fibers  and  epoxy  matrix,  along  with  Ramberg-Osgood  parameters 
for  epoxy  matrix  used  in  modelling  response  of  boron/epoxy  laminates  tested  by  Petit  and 
Waddoups  [61 J.  Fiber  volume  fraction  =  0.5 
Stiffness  reduction  factors  Dfn,Dmn  ,  Dm22 ,  =  0.001 


Fig.  5.19  Shear  stress-strain  response  of  boron/epoxy  lamina  as  modelled  by 
the  proposed  material  model,  and  by  Petit  and  Waddoups. 
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Tensile  Stress,  (Ksi)  “  Tensile  Stress,  °n  (Ksi) 
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Tensile  Strain,  ?n  (%) 

1.20  Response  of  boron  /epoxy  lamina  In  axial  tension,  as  modelled  by  the 
material  model,  and  by  Petit  and  Waddoups. 
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Tensile  Strain,  €22  (%) 

Fig.  5.21  Response  of  boron/epoxy  lamina  In  transverse  tension,  as  modelled 
by  the  material  model,  and  by  Petit  and  Waddoups. 
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Compressive  Strain,  eu  (%) 


Fig.  5.22  Response  of  boron/ epoxy  lamina  in  axial  compression,  as  modelled 
by  the  material  model,  and  by  Petit  and  Waddoups. 


Compressive  Strain,  622  (%) 

Fig.  5.23  Response  of  boron  /epoxy  lamina  In  transverse  compression,  as  modelled 
by  the  material  model,  and  by  Petit  and  Waddoups. 
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Fig.  5.24  Response  of  a  [0/90]s  boron/epoxy  laminate  In  uniaxial  tension 
as  predicted  by  the  material  model,  along  with  prediction  and 
experimental  results  from  Petit  and  Waddoups. 


Tensile  Strain,  eii  (%) 


Fig.  5.25  Response  of  a  [+301-30]$  boron/epoxy  laminate  In  uniaxial  tension 
as  predicted  by  the  material  model,  along  with  prediction  and 
experimental  results  from  Petit  and  Waddoups. 
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Tensile  Strain,  Su  (%) 


Fig.  5.26  Response  of  a  [+60l-60]s  boronlepoxy  laminate  in  uniaxial  tension 
as  predicted  by  the  material  model,  along  with  prediction  and 
experimental  results  from  Petit  and  Waddoups. 


Tensile  Strain,  ?n  (%) 

Fig.  5.27  Response  of  a  [+45I-45I903Js  boronlepoxy  laminate  in  uniaxial 

tension  as  predicted  by  the  material  model,  along  with  prediction  and 
experimental  results  from  Petit  and  Waddoups. 
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Compressive  Stress,  5,,  (Ksl)  *  Compressive  Stress,  °n  (Ksl) 


Compressive  Strain,  ?u  (%) 

i.28  Response  of  a  [0190] s  boron/epoxy  laminate  In  uniaxial  compression  as 
predicted  by  the  material  model,  along  with  prediction  and  experimental 
results  from  Petit  and  Waddoups. 
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Compressive  Strain,  ?n  (%) 

Fig.  8.29  Response  of  a  [+201-20]$  boron f epoxy  laminate  In  uniaxial  compression 
as  predicted  by  the  material  model,  along  with  prediction  and  experimental 
results  from  Petit  and  Waddoups. 
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Compressive  Strain,  ?u  (%) 


Fig.  5.30  Response  of  a  [+30/-30]s  boron/epoxy  laminate  In  uniaxial  compression 
as  predicted  by  the  material  model,  along  with  prediction  and  experimental 
results  from  Petit  and  Waddoups. 


Compressive  Strain,  ?n  (%) 


Fig.  5.31  Response  of  a  [+601-60]$  boron/epoxy  laminate  In  uniaxial  compression  as 
predicted  by  the  material  model,  along  with  prediction  and  experimental 
results  from  Petit  and  Waddoups. 
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CHAPTER  6 


ANALYSIS  OF  LAMINATED  STRUCTURES 


The  material  model  was  incorporated  in  the  finite  element  analysis  package  PO¬ 
LO-FINITE  [96]  developed  at  die  University  of  Illinois.  Since  the  material  model  is 
linked  to  the  structural  analysis  program  through  a  standard  interface,  a  wide  variety  of 
elements  in  POLO-HNITE’s  element  library  are  readily  accessible.  However,  in  this 
study  only  the  2-dimensional  and  3-dimensional  isoparametric  elements  were  employed. 
A  number  of  example  problems  were  solved,  and  results  are  presented  for: 

•  A  thick-walled  laminated  cylinder  in  compression  and  in  circumfer¬ 
ential  bending. 

•  Notched  metal  matrix  laminates  in  tension. 

•  Bending  of  a  simply  supported  plate. 

•  Compression  of  a  notched  laminate  up  to  the  point  of  failure,  using 
micromechanics-based  criteria  for  damage  detection  and  loss  of 
stiffness. 

The  predictions  from  the  proposed  material  model  are  compared  with  experimen¬ 
tal  data  and  analytical  predictions  from  sources  in  the  literature. 

6.1  Response  of  Thick-walled  Cylinder 
in  Compression  and  Bending 

The  geometry  of  the  thick  cylinder  under  consideration  is  shown  in  Fig.  6.1.  It  is 
similar  to  the  cylinders  tested  at  the  David  Thylor  Research  Center  (DTRC)  (Garala  [97]) 
in  an  ongoing  program  to  determine  the  compression  behaviour  of  laminated  cylinders. 
The  loading  cases  which  are  considered  are  shown  in  Fig.  6.2:  axisymmetric  pressure  pro¬ 
ducing  compressive  response  in  the  cylinder  wall;  and  an  ovaling  pressure  with  pzcos29 
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variation  which  produces  primarily  circumferential  flexure.  A  state  of  plane  strain  is  im¬ 
posed  in  the  cylinder  axis  direction  in  both  cases.  Exact  elasticity  solution  for  a  cylinder 
composed  of  anisotropic  laminae  subjected  to  axisymmetric  pressure  is  given  by  Lekhnits- 
kii  [98].  For  cross-ply  laminated  cylinders  subjected  to  harmonic  loading  around  the  cir¬ 
cumference  exact  elasticity  solutions  were  presented  by  Ren  [99].  The  results  from  the 
material  modelling  procedure  will  be  compared  with  these  exact  elasticity  solutions. 

The  cylinder  consists  of  96  plies  of  AS4/3501-6  prepreg  tape  arranged  in  the  stack¬ 
ing  sequence  [02/90]16S .  The  experimentally  observed  properties  of  the  prepreg  tape  were 
reported  by  Camponeschi  [72,73]  and  are  listed  in  Thble  6.2.  Also  listed  in  Thble  62  are 
the  elastic  properties  of  the  prepreg  tape  as  predicted  by  the  micro-model  portion  of  the 
material  model.  The  assumed  properties  of  the  graphite  fibers  and  the  epoxy  matrix  are 
given  in  Thble  6.1.  The  micro-model  predicts  the  lamina  properties  quite  well;  better  fits 
with  experimental  results  can  be  obtained  by  further  tuning  the  fiber  and  the  matrix  prop¬ 
erties. 


Thble  6.3  lists  die  equivalent  elastic  properties  of  a  typical  repeating  sublaminate 
(  [O2/90] )  identifiable  within  the  laminate  stacking  sequence.  The  results  labelled  “Sub¬ 
laminate  Model”  are  obtained  by  using  the  experimentally  observed  lamina  properties  as 
an  input  to  the  sublaminate  model.  The  results  labelled  “Material  Model”  are  from  the 
complete  material  model,  consisting  of  both  micro-model  and  sublaminate  model,  with 
the  constituent  properties  of  Thble  6.1  as  input. 

The  two  loading  cases  described  above  are  analyzed  using  the  finite  element  meth¬ 
od  in  conjunction  with  the  material  modelling  procedure.  The  results  from  the  material 
modelling  procedure  are  referred  to  as  “Homogenized  responses”  since  the  material 
modelling  procedure  is  based  on  homogenization  of  stresses,  strains  and  the  properties  of 
the  material  comprising  the  cylinder.  Exact  homogenized  responses  are  obtained  by  em¬ 
ploying  the  material  modelling  procedure  in  an  exact  elasticity  formulation.  The  accuracy 
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of  the  finite  element  solution  is  assessed  by  comparing  the  finite  element  results  to  the 
exact  homogenized  responses.  To  examine  the  implications  of  the  homogenization  pro¬ 
cess,  the  homogenized  responses  are  compared  to  the  results  obtained  from  the  exact  ply- 
by-ply  elasticity  solution. 

The  finite  element  solution  was  carried  out  by  providing  as  input  the  properties  of 
the  graphite  fibers  and  the  epoxy  matrix  along  with  the  stacking  sequence  of  the  typical 
repeating  sublaminate  ([02/90]) .  The  material  model  generated  the  response  of  the  equiv¬ 
alent  material  representing  the  typical  repeating  sublaminate  via  the  micro-model  and 
the  sublaminate  model.  The  exact  homogenized  responses  were  obtained  by  assuming 
that  the  cylinder  is  composed  of  a  homogeneous  material  with  the  properties  same  as 
those  of  equivalent  material  representing  the  ([O2/9O))  sublaminate. 

6.1.1  Through-thickness  Stress  and  Strain  Distributions 

Figs.  63  and  6.4  show  the  through-thickness  stress  and  strain  distributions  in  the 
case  of  axisymmetric  compressive  loading.  In  this  case  only  a  small  circumferential  seg¬ 
ment  of  the  cylinder  was  modelled,  since  it  is  an  axisymmetric  problem.  For  the  cylinder  in 
bending ,  one-quarter  of  the  cylinder  was  modelled.  The  finite  element  mesh  employed  is 
shown  in  Fig  6.5:  two  elements  are  provided  through  the  thickness  in  order  to  accurately 
capture  the  expected  parabolic  distribution  of  the  shear  stresses  and  strains.  The  stress 
and  strain  distributions  for  the  case  of  cylinder  in  bending  are  shown  in  Figs.  6.6  and  6.7. 

For  the  stress  and  strain  quantities  that  are  continuous  through  the  laminate  thick¬ 
ness  (eg ,  o>and  Tf$ ),  the  homogenized  responses  are  almost  identical  to  the  response  cal¬ 
culated  by  exact  ply-by-ply  solution.  Furthermore  these  continuous  stress  and  strain 
quantities  are  not  only  continuous  at  the  laminae  interfaces  but  also  smoothly  varying 
throughout  the  laminate  thickness. 

For  the  discontinuous  stress  and  strain  quantities  (o$  ,€r  and  €#),  it  is  of  interest  to 
note  that  although  they  may  vary  markedly  from  ply  to  ply,  the  distribution  for  each  family 
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of  plies  with  the  same  orientation  is  smooth,  forming  two  envelopes  of  response.  The  ho¬ 
mogenized  response  is  a  weighted  mean  of  the  exact  ply-by-ply  response.  The  homogeni¬ 
zation  procedure  provides  excellent  estimates  of  the  exact  ply-by-ply  responses  of  continuous 
stress  and  strain  quantities.  The  complementary  discontinuous  stresses  and  strains  can  then 
be  recovered  using  the  individual  ply  properties. 

6.2  Response  of  Notched  Metal  Matrix 
Laminates  in  Tension 

The  material  model  is  used  to  predict  the  response  of  notched  metal  matrix  lami¬ 
nates.  Experimental  results  for  notched  laminates  were  presented  by  Shukow  [100].  The 
geometry  of  the  test  specimens  and  the  strain  gauge  locations  are  shown  in  Fig.  6.8.  Hie 
specimens  were  formed  from  br  ron/aluminum  laminae  arranged  in  stacking  sequences  of 
[+  45/  — 45]  ^  and  [0/  +  <*5/-45/90]s  respectively.  The  response  of  these  notched  lami¬ 
nates  was  investigated  analytically  by  Arenburg  [70],  who  used  Aboudi’s  Method  of  Cells 
[7]  in  conjunction  with  First  Order  Shear  Deformation  Plate  Theory  for  laminates  to  pre¬ 
dict  the  laminate  response.  He  used  Bodner-Partom  unified  viscoplasticity  theory  [69]  to 
characterize  die  aluminum  matrix.  In  the  present  study  the  matrix  material  is  character¬ 
ized  using  Ramberg-Osgood  relations  (Chapter  2,  Section  2.6.1),  as  illustrated  in  the  next 
section. 

6.2.1  Characterization  of  Boron/Aluminum  Lamina 

Since  micromechanics  is  employed  here  to  predict  the  lamina  behaviour,  material 
characterization  in  this  context  entails  determination  of  the  appropriate  properties  of  the 
constituent  fiber  and  matrix  phases  such  that  the  lamina  behaviour  is  accurately  predicted. 

The  elastic  properties  of  boron  fibers  and  aluminum  matrix  provided  by  Arenburg 
[70]  are  listed  in  Thble  6.4.  The  task  of  determining  appropriate  Ramberg-Osgood  pa¬ 
rameters  for  the  aluminum  matrix  such  that  lamina  nonlinearity  is  accurately  portrayed 
was  accomplished  in  the  following  way.  Fig.  6.10  shows  the  response  of  a  [+  45/  -  45]  ^ 
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unnotched  laminate  in  tension,  Appropriate  Ramberg-Osgood  parameters  were  as¬ 
sumed  for  the  aluminum  matrix  to  match  the  predicted  response  with  the  experimental 
data.  The  employed  Ramberg-Osgood  parameters  are  listed  in  Ihble  6.4;  Fig.  6.9  shows 
the  uniaxial  stress-strain  response  of  the  aluminum  matrix.  Before  proceeding  with  the 
analysis  of  the  notched  laminates  some  additional  examples  are  presented  for  unnotched 
laminates  to  show  that  the  material  characterization  is  not  limited  to  the  case  examined  in 
Fig.  6.10. 

Fig.  6.11  shows  the  predicted  response  of  a  [0/  +  45/-45/90]s  laminate.  The  ex¬ 
perimental  results  suggest  a  change  in  stiffness  at  about  10  ksi.  This  is  accurately  predicted 
by  the  material  model  and  is  due  to  the  softening  of  the  aluminum  matrix  in  the  90°  plies. 
Fig.  6.12  shows  the  response  of  a  [0/90]  ^  laminate.  A  definite  knee  is  again  observed  in 
the  predicted  response  due  to  the  yielding  of  the  90°  ply.  The  agreement  between  pre¬ 
dicted  and  the  experimentally  observed  responses  is  very  good  up  to  80  ksi,  after  which  the 
laminate  failed.  Thble  6.6  presents  the  ‘'predicted”  initial  elastic  moduli  of  the 
[+  45/ -45]  ^  and  [0/  +  45/-45/90]s  laminates.  The  predictions  from  the  material  mod¬ 
el  are  in  very  good  agreement  with  the  analytical  predictions  from  Arenburg  [70]  and  the 
experimental  data  from  Shukow  [100]. 

6.2.2  Response  of  [+45/-45]2*  Notched  Laminate 

The  geometry  of  the  test  coupon  is  shown  in  Fig.  6.8.  The  symmetry  of  the  test 
specimen  is  exploited  so  that  only  a  quarter  of  the  specimen  is  modelled  (Fig.  6.13)  using 
20  node  solid  elements  with  2x2x3  Gauss  rule 

Fig.  6.14  shows  the  axial  strain  en  at  the  edge  of  the  hole  and  at  the  far-field  strain 
gauge  locations,  as  a  function  of  applied  stress.  Experimental  results  from  Shukow  [100] 
are  also  shown.  The  analytical  results  for  far-field  strain  are  in  excellent  agreement  with 
the  experimental  results.  For  the  strain  at  the  edge  of  the  hole,  the  analytical  results  from 
the  present  study  are  in  veiy  good  agreement  with  the  analytical  results  of  Arenburg  [70]; 
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however,  both  sets  of  analytical  results  start  diverging  from  the  experimental  results  in  the 
vicinity  of  0.6%  strain.  Beyond  this  point  both  analytical  models  start  underestimating  the 
axial  strain.  Arenburg  attributed  this  to  geometric  nonlinearity  effects  and  possible  dela¬ 
minations,  which  are  not  accounted  for  in  the  analyses.  Another  source  of  the  discrepancy 
could  be  fiber/matrix  debonding,  which  has  also  not  been  accounted  for  in  the  analysis. 

6.2.3  Response  of  [0/+45/-45/90]s  Notched  Laminate 

The  finite  element  mesh  employed  to  analyze  this  laminate  is  shown  in  Fig.  6.15.  It 
is  more  refined  in  the  region  near  the  hole  because  preliminary  investigations  indicated 
that  the  stress  gradients  near  the  hole  were  much  greater  here  than  in  the  [+  45/  -45]^ 
notched  laminate. 

Fig.  6.16  shows  the  computed  strains  at  the  edge  of  the  hole  and  at  the  far-field 
strain  gauge.  The  far-field  strain  is  in  good  agreement  with  the  analytical  predictions 
from  Arenburg  [70],  and  is  essentially  the  same  as  the  stress-strain  response  of  the 
[0/  +  45/ — 45/90]  s  unnotched  laminate  (Fig.  6.11).  The  axial  strain  at  the  edge  of  the  hole 
is  overestimated  compared  to  analytical  results  of  Arenburg  and  the  experimental  results 
of  Shukow.  A  possible  reason  for  the  discrepancy  between  the  analytical  predictions  is 
that  different  constitutive  descriptions  were  employed  in  the  two  studies  to  describe  the 
matrix  material  nonlinearity.  Arenburg  used  Bodner-Partom  unified  viscoplasticity 
theory  to  describe  the  matrix  material,  whereas  here  Ramberg-Osgood  relations  are  used 
for  characterizing  the  matrix  nonlinearity. 

6.3  Bending  of  Simply  Supported 
Square  Plate 

Analytical  results  for  the  response  of  a  simply  supported  square  plate  subjected  to 
a  uniformly  distributed  transverse  load  (Fig.  6.17)  were  presented  by  Arenburg  [70].  The 
plate  is  10.0  inches  square  and  its  thickness  is  0.1  inch,  resulting  in  a  span  to  thickness  ratio 
of  100,  so  that  it  can  be  characterized  as  a  thin  plate.  The  plate  is  made  of  boron/alumi- 
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num  laminae;  the  fibers  in  all  the  laminae  are  orientated  in  a  single  direction  ( i.e.  a  [0°^ 

laminate).  The  properties  of  the  boron/aluminum  lamina  (Thble  6.4)  are  same  as  the  ones 
used  to  make  the  previously  investigated  notched  and  unnotched  laminates.  Only  a  quad¬ 
rant  of  the  plate  is  modelled  for  finite  element  analysis,  using  Q3DISOP  solid  elements. 
The  finite  element  mesh  and  the  imposed  boundary  conditions  are  also  shown  in  Fig.  6.17. 
The  effects  of  geometric  nonlinearity  were  not  by  accounted  Arenburg,  and  are  neglected 
in  this  study  as  well. 

Fig.  6.18  shows  the  center  deflection  of  the  plate  as  a  function  of  the  applied  trans¬ 
verse  load.  Results  are  presented  from  the  study  conducted  by  Arenburg,  and  from  the 
present  study  in  a  number  of  cases  with  different  meshes  and  Gauss  integration  rules.  In 
these  cases  the  parameters  varied  were,  the  number  of  elements  through  the  thickness, 
and  the  order  of  Gauss  integration  rule;  the  cases  considered  can  be  classified  as  follows: 

•  1  Element  through  the  thickness  with  2x2x2  Gauss  rule  (case  (1)). 

•  1  Element  through  the  thickness  with  2x2x3  Gauss  rule  (case  (2)). 

•  2  Elements  through  the  thickness  with  2x2x3  Gauss  rule  (case  (3)) 

The  maximum  number  of  Gauss  points  that  can  be  employed  in  the 
thickness  direction  with  Q3DISOP  element  is  restricted  to  3,  so  in 
order  to  have  additional  Gauss  points  in  the  thickness  direction,  2 
elements  are  employed  through  the  thickness  with  2x2x3  Gauss 
rule. 

An  interesting  phenomenon  is  observed  in  the  results  presented  in  Fig.  6.18;  i.e. 
the  results  from  the  different  combinations  of  meshes  and  integration  rules  do  not  seem 
to  be  converging  monotonically  to  an  exact  solution,  rather  the  convergence  seems  to  be 
non-monotonic.  This  is  a  result  of  material  nonlinearity  and  the  sensitivity  of  the  com- 
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puted  response  to  the  location  of  the  integration  points  associated  with  the  different 
Gauss  rules. 

Fig.  6.19  shows  the  location  of  the  Gauss  sampling  points  throughout  the  thickness 
in  the  various  cases  considered.  If  the  cross-section  is  undergoing  flexure,  the  flexural 
stresses  will  be  higher  near  the  top  and  bottom  surfaces  of  the  plate.  Case  (1)  will  yield  a 
lower  stiffness  value  for  the  element  than  case  (2),  because  in  case  (2)  the  sampling  points 
are  nearer  to  the  top  and  bottom  surfaces  of  the  plate  where  they  detect  the  onset  of  non¬ 
linearity  earlier  than  in  case  (1).  The  integration/sampling  points  near  the  plate  mid-sur¬ 
face  do  not  detect  significant  straining  and  therefore  do  not  influence  the  results  much. 
Hence  we  see  more  softened  response  for  case  (2)  as  compared  to  case  (1).  In  case  (3) 
there  is  an  additional  layer  of  Guass  points  between  the  Gauss  points  near  the  plate  mid¬ 
surface  and  the  top  and  bottom  surfaces  of  the  plate.  The  Gauss  points  in  the  additional 
layer  are  sufficiently  far  from  the  top  and  bottom  of  the  plate  and  detect  the  onset  of  non¬ 
linearity  later  than  the  extremal  Gauss  points  in  case  (2).  This  explains  the  slightly  stiffer 
response  predicted  in  case  (3)  (Fig.  6.18).  The  results  from  case  (3)  are  accepted  as  more 
accurate  than  cases  (1)  and  (2),  since  the  number  of  Gauss  points  in  the  thickness  direction 
is  greater  and  their  distribution  through  the  thickness  is  more  reasonable. 

The  maximum  load  intensity  is  12  ksi.  At  this  load  level,  the  analytical  results  of 
Arenburg  predict  the  deflection  at  center  of  the  plate  to  be  0.41 1  in.,  whereas  the  accepted 
solution  in  the  present  study  (2  elements  through-thickness  with  2x2x3  Gauss  rule)  pre¬ 
dict  the  center  deflection  to  be  0.446  in. 

Fig.  6.20  shows  the  variation  in  the  maximum  stresses  at  the  plate  top  surface  as  the 
applied  load  increases.  The  maximum  crn  and  &&  stresses  occur  at  the  plate  center, 

whereas  the  maximum  shear  stress  occurs  at  the  comers  of  the  plate.  The  agreement 
with  the  results  presented  by  Arenburg  is  fairly  good.  The  present  analysis  shows  more 
softening  in  the  022  response  than  predicted  by  Arenburg.  This  is  accompanied  by  an 
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elevation  of  the  on  stress  in  order  to  balance  the  applied  load.  The  calculated  shear 
stresses  xn  from  the  two  studies  are  in  excellent  agreement. 

6.4  Compressive  Strength  of  Notched  Laminate 

Chang  and  Lessard  [102,103]  conducted  an  experimental  and  analytical  investigation 
into  the  compressive  strength  of  notched  laminates  composed  of  T300/BP976  graphite/e¬ 
poxy  laminae.  The  lamina  was  modelled  as  an  orthotropic  material  exhibiting  nonlinearity 
only  in  the  in-plane  shear  response,  using  the  relation  suggested  by  Hahn  and  Tfcai  [62,63]. 
Failure  criteria  for  predicting  die  in-plane  failure  modes  ofa  lamina  were  also  given.  After  a 
lamina  Med,  its  stiffness  was  reduced  to  reflect  the  damage.  The  laminate  response  was 
obtained  by  assembling  the  individual  laminae  within  the  frame  work  of  Classical  laminated 
Plate  Theory. 

The  material  modelling  procedure  developed  in  the  present  study  is  used  to  predict 
the  compression  response  of  one  of  the  notched  laminates  studied  by  Chang  and  Lessard 
[102,103].  The  geometry  of  the  selected  notched  laminate  is  shown  in  Fig.  6.21.  The  test 
specimen  has  an  extensometer  attached  to  it  which  measures  the  shortening  over  the  gau¬ 
ge-length  of  1.0  inch. 

The  characterization  of  the  T300/BP976  lamina  is  described  in  the  next  section. 

6.4.1  Characterization  of  T300/BP976  Lamina 

The  elastic  moduli  of  the  T300/BP976  lamina  are  given  by  Chang  and  Lessard 
[102].  From  these  lamina  properties  and  the  information  available  in  the  literature,  the  in 
situ  elastic  moduli  of  the  constituent  graphite  T300  fibers  and  the  BP976  epoxy  matrix 
were  estimated,  and  are  given  in  Ihble  6.7.  The  properties  of  the  T300/BP976  lamina  used 
by  Chang  and  Lessard  are  given  in  Ihble  6.8.  The  predictions  of  lamina  properties  by  the 
material  model  are  also  presented  in  Ihble  6.8.  Fig.  6.22  shows  the  shear  response  of  the 
lamina  according  to  Chang  and  Lessard.  Modified  Rambe r g-Osgood  relations  (Chapter 
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2,  Eqs.  (2.24-2  JO))  are  used  to  model  this  response;  the  parameters  are  presented  in 
Thble  6.7. 

The  micromechanical  failure  criteria  employed  in  die  analysis  along  with  the  proce¬ 
dure  for  damage  modelling  has  been  described  in  (Chapter  2,  Sections  2.7-2.8).  The  micro¬ 
failure  criteria  predict  the  strengths  of  the  lamina  very  well  (Table  6.8).  The  only  exception  is 
the  axial  compressive  strength,  for  which  the  fiber  micro-budding  criterion  grossly  over-pre- 
dicts  the  actual  compressive  strength  of  the  lamina.  This  discrepancy  may  be  a  result  of  fiber 
misalignment  in  the  lamina;  it  may  also  be  that  a  failure  mechanism  other  than  the  fiber 
micro-buckling  such  as  fiber  compressive  failure  precipitates  the  failure. 

6.4.2  Response  of  [(+45/-45)6]s  Notched  Laminate 

The  notched  laminate  shown  in  Fig  6.21  is  analyzed  using  the  the  mesh  shown  in  Fig. 
6.23.  The  properties  of  the  of  fiber  and  the  matrix  along  with  die  stacking  sequence  of  the 
typical  repeating  sublaminate  (+4S/-45)  is  provided  as  input  for  the  material  model. 

Fig.  624  shows  the  load-shortening  response  of  the  notched  laminate.  The  response 
is  quite  nonlinear  since  the  +45/-45  stacking  sequence  produces  primarily  shear  stresses 
within  the  laminae.  The  analysis  does  not  predict  any  damage  up  to  2400  lbs,  at  which  point 
fiber  micro-budding  is  predicted  around  the  periphery  of  the  hole.  Application  of  additional 
load  results  in  matrix  shearing  failure  in  elements  lying  in  die  vicinity  of  die  45°  direction 
starting  from  the  region  of  maximum  stress  concentration.  The  effect  of  matrix  shearing  fail¬ 
ure  is  modelled  by  reducing  the  lamina  in-plane  shear  stiffness.  The  loss  of  lamina  shear 
stiffness  results  in  prediction  of  fiber  micro-buckling.  Figs.  625-621  show  the  predicted 
state  of  damage  at  various  load-levels.  The  predicted  ultimate  load  is  3320  lbs,  which  is 
about  11%  more  than  the  collapse  load  of  3000  lbs  predicted  by  Chang  and  Lessard  [102]. 
The  present  investigation  employed  maximum  stress  failure  criteria  at  the  micromechanical 
level:  It  is  speculated  that  application  of  interaction  type  failure  criteria  at  the  microme¬ 
chanical  level  would  result  in  further  improvement  in  the  predicted  results. 
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(a)  Compression  (&)  Bending 

Flfl*  6-2  (a)  Cylinder  In  compression,  (b)  Cylinder  In  circumferential  bending. 


E\  E2 
(Mpsl)  (MpSi) 


G 12 
(Mpsl) 


Vl 


V23 


Graphite 

27.0 

2.5 

5.0 

0.30 

0.25 

Epoxy 

0.728 

0.728 

0.26 

0.40 

0.40 

Table  6.1  Elastic  properties  of  graphite  fibers  and  epoxy  matrix. 


Source 

Ei 

(Mpsl) 

Ei 

(Mpsl) 

Gj2 

(Mpsl) 

G23 

(Mpsl) 

Vi 

V23 

Experiment, 

Camponeschl 

16.48 

1.40 

0.87  1 

0.55  2 

0.334 

0.54 

[72,73] 

Micro-model 

16.50 

1.56 

0.87 

0.53 

0.337 

0.468 

1  From 

[±  45]2S 

tension  test 

2  From  literature 


Table  6.2  Properties  of  AS4/3501-6  prepreg  tape. 
Fiber  volume  fraction  =  0.60 


Source 

Ei 

Ei 

e3 

Gn 

C?i3 

G23 

V12 

V13 

V23 

(Mpsl) 

(Mpsl) 

(Mpsl) 

(Mpsl) 

(Mpsl) 

(Mpsl) 

Experiment, 

Camponeschl 

[72,73] 

11.63 

— 

— 

— 

— 

— 

0.069 

0.469 

— 

Sublaminate 

Model 

11.53 

6.47 

1.80 

0.87 

0.73 

0.63 

0.073 

0.448 

0.519 

Material 

Model 

11.60 

6.58 

1.85 

0.82 

0.69 

0.60 

0.08 

0.455 

0.467 

Table  6.3  Equivalent  elastic  properties  of  [Oa/90]  repeating  sublaminate. 
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10,000  psi 


-15  -10  -5  0  5  10  15 

Radial  Stress  aT  (ksi) 

Fig.  6.3  Response  of  laminated  cylinder  In  compression, 
continuous  quantities  (e$  ,or ) 


[02/90]^,  laminate 

Finite  Element  1  ,,  .  . 

Exact  f  Homogenized 

Ply-by-Ply 


Typical  Element. 
Q3DIS0P  (20  Node  Brick) 


36  Elements 
355  Nodes 


Fig.  6.5  Finite  element  mesh  employed  In  calculating  the  response  of 
laminated  cylinder  In  bending 
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Fig.  6.7  Response  of  laminated  cylinder  In  bending, 
discontinuous  quantities  (o0  ,€r  ,yu) 


Strain  Guage  at  Edge  of  Hole 

/—  Far-Fleld  Strain  Guage 


E 

V 

To 

a 

b 

n 

(Mpsi) 

(Ksi) 

Boron 

58.0 

0.2 

— 

— 

— 

— 

Aluminum 

8.69 

0.25 

6.0 

0 

1.0 

5 

Table  6.4  Elastic  properties  of  boron  fibers  and  aluminum  matrix,  along  with  Ramberg- 
Osgood  parameters  for  aluminum  matrix.  Fiber  volume  fraction  =  0.44 


En 

E22 

g12 

*12 

(Mpsi) 

(Mpsi) 

(Mpsi) 

Material 

Model 

30.40 

16.70 

6.51 

0.223 

Arenburg  [70] 

30.39 

16.03 

6.51 

0.229 

Experiment 
Shukow  [100] 

30.59 

19.20 

6.50 

0.218 

Table  6.5  Elastic  properties  of  boron! aluminum  lamina. 


Tensile  Strain  €u  (%) 

Fig.  6.9  Back-calculated  uniaxial  stress-strain  response  of 
In  situ  aluminum  matrix. 
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Strain, €n  (%) 


Fig.  6.10  Response  of  1*451-45] &  boron! aluminum  laminate  as  predicted 
by  the  material  model,  along  with  predictions  from  Arenburg,  and 
experimental  data  from  Shukow. 


Strain,  €11  (%) 


Fig.  6.11  Response  of  [0I*45I-45I90]S  boron! aluminum  laminate  as  predicted 
by  the  Material  Model,  along  with  predictions  from  Arenburg,  and 
experimental  data  from  Shukow. 
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Fig.  6.12  Response  of  [0190] 2s  boronlalumlnum  laminate  as  predicted  by  the 

material  model,  along  with  predictions  from  Aboudi,  and  experimental 
data  from  Bigelow  et  al. 


En  (Mpsi) 

^12 

Laminate 

Material 

Arenburg  [70] 

Experiment, 

Material  Arenburg  [70]  Experiment, 

Model 

Shukow  [100] 

Model 

Shukow  [100] 

[*  ^2S 

17.77 

17.67 

17.36 

0.366 

0.358 

0.373 

[0/±45/90]s 

20.92 

20.70 

21.25 

0.253 

0.249 

0.267 

Table  6.6  Elastic  properties  of  boronlalumlnum  laminates. 
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Fig.  8.13  Finite  element  mesh  employed  In  the  analysis  of  [+451-45^  notched  coupon. 


Stress  5  (Ksi) 


Strain,  f a  (%) 

Fig.  6.14  Response  of  [+451-45]  23  notched  boron! aluminum  laminate. 


L44 


Fig.  8.15  Finite  element  mesh  employed  In  the  enelysls  of  [01+451-45190] s  notched  coupon. 


Stress  o  (Ksi) 


X 


0.0  0.2  0.4  0.6  0.8  1.0  1.2 

Strain,  en  (%) 


Fig.  6.16  Response  of  [0+45l-45/90]s  notched  boron! aluminum  laminate. 
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Transverse  Load,  p  (psi) 


Center  Deflection,  d  (in.) 


Fig.  6.18  Center  deflection  of  simply  supported  square  plate. 
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Case  (1)  Case  (2)  Case(3) 

1  Element  1  Element  2  Elements 

2x2x2  Gauss  rule  2x2x3  Gauss  rule  2x2x3  Gauss  rule 


Fig.  6.18  Location  of  Gauss  Integration  points  In  the  thickness  direction, 
for  the  various  cases  considered. 
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m**T12here 


W=  1.0 
Ail  dimensions  are  in  inches 


Fig.  6.21  Geometry  of  the  notched  laminate  tested  by  Chang  and  Lessard  (102,1031, 
stacking  sequence  [(+  45/- 45)6]$ 
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Ei 

Vl 

E2 

v23 

Gn 

To 

b 

n 

if> 

s(m) 

(Mpsi) 

(Mpsi) 

(Mpsi) 

(Ksi) 

(Ksi) 

(Ksi) 

(Ksi) 

33.8 

0.2 

3.0 

0.42 

1.7 

— 

— 

— 

347.0 

— 

— 

0.925 

0.25 

0.925 

0.25 

0.37 

9.25 

1.0 

4 

7.5 

42.0 

17.0 

Table  6.7  Elastic  properties  of  7300  graphite  fibers  and  BP976  epoxy  matrix,  along  with  Ramberg- 
Osgood  parameters  for  the  matrix  used  In  modelling  the  response  of  T300IBP976  lamina. 
Fiber  volume  fraction  -  0.66 

Stiffness  reduction  factors  Dfn ,  Dmn  =  g.01  >  ^44  =  0.25 


Source 

Ei 

V| 

Ei 

V23 

Gn 

Y, 

Ye 

5 

(Mpsi) 

(Mpsi) 

(Mpsi) 

(Ksi) 

(Ksi) 

(Ksi) 

(Ksi) 

(KSI) 

22.7 

0.23 

1.88 

— 

1.01 

220.0 

231.0 

6.46 

36.7 

15.5 

22.63 

0.213 

1.94 

0.367 

0.89 

220.0 

590.0 

6.5 

37.0 

15.5 

*  The  predicted  compressive  strength  as  calculated  from  micro-buckling  criteria  is  much  higher, 
probably  because  of  some  bther  failure-mechanism,  such  as  fiber  compressive  failure. 

Table  6.8  Elastic  properties  and  strengths  of  graphite  T300/BP976  lamina  as  given  by  Chang  and 
Lessard  [102],  and  predictions  from  the  material  model. 
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Fig.  6.22  Shear  stress-strain  response  of  T300/BP976  graphite! epoxy  lamina 
according  to  the  material  model,  and  Chang  and  Lessard  [102], 
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Applied  Load,  P  (lb) 


Extensometer  Measurement,  AE  (in) 


Fig.  6.24  Load-shortening  response  of  [(+  45/ -45)^  notched  laminate. 
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Applied  Load-level  P  =  3320  lbs  (Collapse-load) 


Fiber/Matrix 

Shearing 


Fiber  Buckling 


Fig.  6.27  Sfafe  of  damage  In  [(+  457-45)615  notched  laminate 
at  the  load-level  of  P  *  3320  lbs  (Collapse-load) 
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CHAPTER  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  the  present  study  was  to  develop  general  and  accurate  proce¬ 
dures  for  analysis  of  fiber-reinforced  laminated  structures  which  also  allow  the  use  of 
standard  displacement-based  finite  elements.  The  stress  and  strain  fields  in  laminated 
composites  can  exhibit  discontinuities;  procedures  were  sought  which  take  this  phenome¬ 
non  into  account 

A  micromechanical  constitutive  theory  was  used  to  characterize  the  composite 
lamina.  The  material  model  has  a  two-level  hierarchical  structure.  At  the  bottom  level, 
the  Micro-model  synthesizes  the  constitutive  information  about  fiber  and  matrix  phases  to 
yield  the  effective  stress-strain  response  of  a  unidirectional  lamina.  At  the  top  level  the 
Sublaminate  model,  using  a  3-D  lamination  scheme,  assembles  the  laminae  forming  the 
smallest  repeating  stack  within  the  laminate  to  deliver  the  effective  stress-strain  response 
of  the  laminate.  Local  stresses  and  strains  in  a  lamina  or  in  fiber  and  matrix  phases  can  be 
recovered  from  the  effective  values  at  any  stage.  The  material  modelling  procedure  en¬ 
ables  the  use  of  standard  displacement  based  finite  elements  for  the  analysis  of  fiber-rein¬ 
forced  laminated  composites. 

In  this  study  the  micromechanical  model  presented  by  Pecknold  [57]  is  employed 
and  enhanced  to  model  the  behaviour  of  unidirectional  lamina.  It  was  found  to  yield  very 
accurate  results  in  spite  of  its  simplicity. 

Several  conclusions  based  on  the  results  of  this  study  are  discussed  next 
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7.1  Conclusions 


General  conclusions  about  the  overall  modelling  procedure  are  presented  first. 
Next,  conclusions  regarding  the  micro-model  are  presented.  Finally,  some  observations 
on  micromechanical  modelling  of  fiber-reinforced  composites  are  presented. 

7.1.1  General 

The  material  modelling  procedure  successfully  enables  the  use  of  standard  displa¬ 
cement-based  finite  elements  for  the  analysis  of  laminated  composites.  The  general 
scheme  has  sufficient  flexibility  in  it  so  that  different  types  of  constitutive  descriptions  for 
the  fiber  and  matrix  constituents  can  be  incorporated  easily.  The  modelling  procedure 
satisfies  the  equilibrium  requirements  on  the  tractions  at  the  laminae  interfaces. 

7.1.2  Micro-model 

The  micro-model  is  the  kernel  on  which  the  accuracy  of  the  overall  modelling 
scheme  depends.  Therefore,  the  micro-model  (Pecknold  [57])  was  investigated  in  detail; 
conclusions  are  presented  below. 

The  micro-model  is  quite  simple  in  its  description,  but  accurately  accounts  for  the 
three-dimensional  matrix  restraining  effects  of  the  fiber  reinforcement.  The  micro-mo- 
del  is  valid  for  general  multi-axial  stress  states.  The  elastic  moduli  of  composite  laminae 
predicted  by  the  micro-model  show  excellent  agreement  with  results  from  detailed  nu¬ 
merical  schemes  and  experimental  data;  it  was  found  to  be  comparable  in  accuracy  to 
Aboudi’s  Method  of  Cells  [7],  but  is  much  simpler.  The  predicted  effective  coefficients  of 
thermal  expansion  were  also  in  very  good  agreement  with  results  from  Aboudi’s  Method 
of  Cells. 

The  micro-model  gives  accurate  results  for  nonlinear  response  of  composite  lami¬ 
nae.  The  results  obtained  compared  very  well  with  experimental  data. 

The  results  for  off-axis  strengths  of  composite  laminae  were  in  very  good  agree¬ 
ment  with  experimental  data.  The  micro-model  predictions  of  initial  yield  surfaces  of 
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metal  matrix  composites  compared  very  well  with  finite  element  based  results  and  with  the 
Method  of  Cells;  the  only  exception  being  the  case  when  the  composite  lamina  is  sub¬ 
jected  to  a  hydrostatic  stress  state,  in  which  case  both  the  micro-model  and  the  Method  of 
Cells  performed  poorly.  However,  the  results  improved  as  the  applied  stress  state  moved 
away  from  the  hydrostatic  stress  state. 

7.1.3  Comments  on  Micromechanical  Modelling 

Micromechanics  based  material  models  are  a  viable  option  for  characterizing  the 
behaviour  of  composite  materials  and  can  be  successfully  used  for  analyzing  structures 
made  from  such  materials.  The  advantage  of  micromechanics  based  material  models  is 
that  they  predict  the  behaviour  of  composites  from  the  fundamental  properties  of  its  basic 
constituents  and  also  represent  the  microstructure  of  the  composite  to  some  extent 
Therefore,  they  can  accommodate  relatively  complex  aspects  of  composite  material  beha¬ 
viour,  such  as  debonding  at  the  fiber/matrix  interface  and  matrix  cracking.  The  complex 
constitutive  description  of  composites  can  be  obtained  by  synthesizing  comparatively  sim¬ 
ple  constitutive  descriptions  of  its  constituents  and  their  interaction. 

The  problems  associated  with  the  material  modelling  procedure  are  those  usually 
associated  with  the  use  of  micromechanics.  For  example  it  is  known  that  the  in  situ  prop¬ 
erties  of  the  matrix  material  can  be  quite  different  from  the  properties  of  matrix  material 
in  bulk.  Direct  measurement  of  the  in  situ  properties  of  the  matrix  and  fiber  are  at  this 
time  (1992)  not  possible,  and  these  properties  have  to  be  back-calculated  from  the  lamina 
properties  determined  from  simple  tests. 

In  this  study  micromechanical  failure  criteria  were  employed  for  determining  fail¬ 
ure  in  the  lamina.  Only  the  in-plane  failure  mechanisms  were  considered  and  the  failure 
mechanisms  such  as  delaminatinns  were  not  considered.  Maximum  stress  failure  criteria 
were  successfully  employed  to  determine  failure  modes  such  as  fiber  fracture  in  tension, 
matrix  cracking  in  transverse  direction,  and  matrix  crushing  in  compression.  For  predic- 


161 


tion  of  axial  compressive  strength  of  lamina  a  micromechanics  based  criterion  suggested 
by  Hahn  and  Williams  [95]  was  investigated,  but  was  found  to  overestimate  the  axial  com* 
pressive  strength. 

In  the  author’s  opinion,  micromechanics  based  material  modelling  procedures 
have  great  potential,  and  as  die  problem  of  scarcity  of  experimental  studies  is  solved  and 
better  experimental  procedures  are  developed,  the  area  of  micromechanics  should  see 
significant  improvement.  Therefore,  micromechanics  based  structural  analysis  proce¬ 
dures  deserve  continued  investigation. 

7.2  Recommendations 

Here,  some  possible  changes  and  enhancements  that  can  be  made  to  die  material 
modelling  procedure  are  presented.  The  suggested  recommendations  concern  mostly  die 
micromechanical  modelling. 

•  In  order  to  investigate  problems  involving  thermal  loading,  strain 
rate  effects  and  creep,  incorporation  of  a  unified  viscoplastidty 
theory,  such  as  the  Bodner-Partom  theory  [69]  would  indeed  be  very 
useful. 

•  In  this  study  relatively  simple  failure  criteria  were  employed  at  the 
micromechanical  level  to  determine  the  onset  of  failure.  A  possible 
improvement  would  be  to  employ  interaction  type  criteria,  which 
take  into  account  the  interaction  of  stresses,  and  at  the  same  time 
provide  information  about  the  mode  of  failure. 

•  The  manner  in  which  the  effect  of  localized  damage  is  incorporated 
is  an  important  consideration.  Accurate  modelling  of  the  effect  of 
damage  is  needed  for  accurate  prediction  of  the  laminate  response 
and  the  ultimate  load  carrying  capacity  of  composite  structures.  The 
field  of  Damage  Mechanics  provides  damage  evolution  laws  and  pro* 
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cedures  for  incorporating  the  resulting  deterioration  and  loss  of  stif¬ 
fness  in  a  material.  In  the  present  study  simple  procedures  were 
employed  to  simulate  the  effect  of  damage.  Incorporation  of  con¬ 
cepts  and  procedures  from  the  area  of  damage  mechanics  within  the 
micro-model  may  result  in  further  improvements  in  the  predictions 
from  the  material  modelling  procedure. 

Deterioration  in  the  bond  between  the  fiber  and  the  matrix  at  the 
interface  can  significantly  affect  the  overall  properties  of  the  com¬ 
posite  lamina.  This  has  been  demonstrated  by  Aboudi’s  Method  of 
Cells.  The  capability  to  model  fiber/matrix  interface  failure  should 
also  be  incorporated  in  the  micro-model. 
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